APPENDIX  C 

Final  Technical  Report 
Grant  F49610-92-J-0110 


Gp 

C 


A 


RADUATE  /AERONAUTICAL  LABORATORIES 


ALIFORNIA  INSTITUTE  OF  lECHNOLOGY 


Te 


Shock  wave  interaction 
in  hypervelocity  flow 

Simon  R.  Sanderson 
1995 


Firestone  Flight  Sciences  Laboratory 

Guggenheim  Aeronautical  Laboratory 

Karman  Laboratory  of  Fluid  Mechanics  and  Jet  Propulsion 


19960408  135 


.1  - 


Pasadena 


report  documentation  page 


form  Approved 
0MB  No.  0704^788 


»  JJL  »» to  rrvf»9«  I  iw  por  womo.  WKJudtn^  IhP  wm  Iv  ftywww  irnt^ioin.  t»»tiog  tf»u  tourm. 

C0me**ti«0  »nd  rtfviowmq  th«  cortt«K>«  <H  m»ormoi*on.  rioorrtiM  burd^^tUnatt  Of  any  other  MOt^  INt 

fn  - «Of  f*duei«  th*s  bwrt*#«.  to  w«h*rtgto«  «o«hioanta  WryktC  Owwrttt  for  mionwatJOh  Operttiom  aod  «*ooftu  1215  iofttnon 

Of1«  o«  •*««  Otfm^Mdwct>ei*fhoK<t(oy^  PC  20503. _ 


1.  AGENCY  USf  OfiLY  (U#vf  bisnk)  I  2. 


REPORT  DATE 


3.  REPORT  TYPE  AND  OATES  COVERED 

Final  Technical  Kept.  1  l/15/9r.--8/31/95 


5.  FUNDING  NUMBERS 

PE  - 
PR  - 


SA  - 

&.  authoaU) 

H,  G.  Hornung  and  B.  Sturtevant 

G  -  F4961O-02-J-O11O 

7,  PERFORMMG  ORGANIZATION  NAME(S)  AND  ADDRESSEES) 

Graduate  Aeronautical  Laboratories 

California  Institute  of  Technology 

1201  E.  California  Blvd,  y 

Pasadena,  CA  91125  ^ 

AFOSR-TR-96 

9.  SPONSOWNe/MONtTOWNG  AGENCY  NAME(S)  AND  AOOAESS(£S)  ^  j 

AFOSR/NA  '  ■ 

no  Duncan  Avenue,  Suite  B115 

Bolling  AFB  DC  20332-0001 

to.  SPONSORING/ MONITORING 

AGENCY  REPORT  NUMBER 

n.  SUPPUMENTARY  NOTES 

121.  DlSTRaUTION/AVAlLABlUTY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is 
unlimited 

13.  ABSTRAa  (M^KimwniOOwOfxh) 

The  following  important  results  were  obtained  from  investigations  of  chemically  reacting  blunt-body  flows  and 
of  shock  wave  interactions  in  hypervelocity  flow: 

•  A  new  view  of  the  binary  scaling  concept  resulted  from  experimental  and  theoretical  work  on  flow  over 


spheres,  in  terms  of  two  dimensionless  parameters. 

•  The  evolution  of  the  vorticity  generated  by  a  curved  shock  wave  could  be  visualized  with  resonant  enhance¬ 
ment  of  streak  lines. 

•  The  effect  of  recombination  in  flow  along  a  streamline  after  a  curved  shock  was  explained. 

•  Previous  predictions  of  greatly  enhanced  heat  flux  by  real-gas  effects  in  shock-on-shock  interactions  were 


not  borne  out. 

•  The  study  provided  extensive  corroboration  of  the  computed  free-stream  conditions  in  the  free-piston  shock 
tunnel  T5. 

•  A  number  of  improvements  of  shock  tunnel  technology  were  made,  especially  regarding  diaphragm  manu¬ 
facture',  nozzle  throat  material,  and  computation  of  the  processes  during  operation. 

•  Improved  diagnostic  techniques  were  generated,  including  low-noise  coaxial  thermocouple  gauges  and  inex¬ 
pensive  high-qtiality  holographic,  interferometry. 


i"l4.  SUBJEa  TERMS 

hypervelocity 
heat  transfer. 


U  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclass  if ied 


flow,  shockwave,  shock  interaction, 
vorticity ,  dissociation 

15.  NUMBER  OF  PAGES 

485 

16.  PRICE  CODE 

IB.  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT 

OF  THIS  PAGE 

OF  ABSTRACT 

Unclass  if ied 

Unclassif ied 

UL 

NSN  7 540-01 -2 80-5 500 


Standard  Form  298  {890104  Draft) 

by  akV  Std  JM-U 

29»-Ot 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE 
COPY  FURNISHED  TO  DTIC 
CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO 
NOT  REPRODUCE  LEGIBLY. 


Shock  wave  interaction  in  hypervelocity  flow 


Thesis  by 

Simon  R.  Sanderson 

In  Partial  Fulfillment  of  the  Requirements 
for  the  Degree  of 
Doctor  of  Philosophy 


California  Institute  of  Technology 
Pasadena,  California 


1995 

(Submitted  May  5,  1995) 


11 


©  1995 

Simon  R.  Sanderson 
All  Rights  Reserved 


Acknowledgements 


iii 


Galcit  is  an  institution  that  has  prospered  over  the  past  60  years  because  of  the  collective  dedication 
of  the  diverse  group  of  people  that  comprise  it.  The  author  wishes  to  express  his  gratitude  to 
this  community  for  the  role  it  has  played  in  the  completion  of  this  work.  Through  their  insight, 
encouragement  and  wise  direction,  Professors  Brad  Sturtevant  and  Hans  Hornung  have  contributed 
immeasurably  to  this  endeavour.  Equally  appreciated  are  the  efforts  of  Jacques  Belanger,  Eric 
Cummings,  Patrick  Germain,  Bahrain  Valiferdowsi  and  Chih-yung  Wen  who  helped  to  build  T5  and 
formed  part  of  the  much  larger  group  whose  cooperation  made  involvement  with  this  laboratory 
rewarding  and  enjoyable. 

Personal  thanks  are  due  to  the  author’s  family  who  fostered  his  penchant  for  science  from  an 
early  age.  Their  support  has  been  unfailing  ever  since.  Christina  has  been  a  kindhearted  encourager 
and  the  truest  friend  one  could  wish  for. 

None  of  this  would  have  been  possible  without  the  beneficence  of  contributors  to  the  Insti¬ 
tute.  The  author  is  grateful  for  assistance  received  from  the  Darryl  G.Greenamyer  Fellowship  and 
C.L.  Powell  Fellowship  funds.  This  work  was  supported  by  AFOSR  Grant  Nos.  F49620-92-J-0110 
and  F49620-93- 1-0338. 


IV 


Abstract 

The  interaction  of  a  weak  oblique  shock  with  the  strong  bow  shock  ahead  of  a  blunt  body  in 
supersonic  flow  produces  extreme  heat  transfer  rates  and  surface  pressures.  Although  the  problem 
has  been  studied  extensively  in  low  enthalpy  flows,  the  influences  of  high  enthalpy  real  gas  effects 
are  poorly  understood.  Existing  perfect  gas  models  predict  greatly  increased  heating  with  increasing 
Mach  number  and  decreasing  ratio  of  specific  heats. 

Experiments  are  conducted  in  a  free  piston  shock  tunnel  to  determine  the  effects  of  thermochem¬ 
istry  on  the  problem  at  high  enthalpy.  The  flow  topology  is  simplified  by  studying  the  nominally 
two-dimensional  flow  about  a  cylinder  with  a  coplanar  impinging  shock  wave.  High  resolution  holo¬ 
graphic  interferometry  is  used  to  investigate  changes  in  the  flow  structure  as  the  location  of  the 
impinging  shock  wave  is  varied.  Fast  response  heat  transfer  gauges  provide  time  resolved  measure¬ 
ments  of  the  model  surface  temperature.  The  data  that  are  obtained  do  not  support  the  existing 
predictions  of  greatly  increased  heat  transfer  at  high  enthalpy. 

A  model  is  developed  to  study  the  thermochemical  processes  occurring  in  the  interaction  region. 
The  phenomenon  arises  because  the  stagnation  streamline  is  forced  to  pass  through  a  system  of 
oblique  shock  waves  that  produce  less  entropy  than  the  undisturbed  bow  shock.  Peak  heating  is 
shown  to  result  from  a  balancing  of  the  strengths  of  the  oblique  shock  waves.  This  condition  is 
demonstrated  to  simultaneously  minimize  the  influence  of  thermochemistry  on  the  flow.  Real  gas 
effects  are  shown  to  become  important  at  lower  Mach  numbers  (<  7.5)  and  for  shock  angles  weaker 
or  stronger  than  that  which  produces  maximum  heating.  The  model  accurately  reproduces  the 
experimental  observations. 

A  nonequilibrium  approximation  is  introduced  that  applies  when  the  oblique  waves  are  weak 
with  respect  to  the  undisturbed  bow  shock.  Within  the  scope  of  the  approximation,  non-monotonic 
behavior  with  the  reaction  rate  is  predicted.  The  reaction  rate  is  not  varied  as  an  independent 
parameter  in  the  current  experiments. 
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Chapter  1  Introduction 

1.1  Motivation 

The  impingement  of  a  weak  oblique  shock  wave  on  the  bow  shock  ahead  of  a  blunt  body  in  steady 
supersonic  flow  is  known  to  cause  extremely  high  local  heat  transfer  rates  and  surface  pressures. 
This  phenomenon  represents  a  severe  constraint  on  the  design  of  hypervelocity  vehicles.  Korkegi 
[26]  has  conducted  an  extensive  review  of  the  literature  and  discusses  the  available  data  from  flight 
test  vehicles.  A  well  documented  example  of  the  problem  occurred  during  the  X-15  flight  research 
program.  Damage  to  the  vehicle  resulted  from  the  impingement  of  an  oblique  shock  wave  on  the 
leading  edge  of  the  lower  vertical  stabilizer  fin  (Burcham  &  Nugent  [2]).  Shock  interactions  that 
occurred  in  the  flow  about  a  model  of  the  X-15  vehicle  are  illustrated  in  figure  1.1.  The  damage 
sustained  by  the  full  scale  vehicle  (Figure  1.2)  demonstrates  the  need  for  an  improved  understanding 
of  the  problem.  The  phenomenon  is  also  observed  in  vehicles  with  a  delta  planform  where  the  bow 
shock  intersects  the  swept  leading  edge  of  the  wing.  Complex  interaction  flow  fields  also  arise 
between  vehicles  and  externally  mounted  tanks  or  boosters. 

1.2  Description  of  the  phenomenon 

The  most  complete  description  of  the  shock  impingement  phenomenon  is  that  due  to  Edney  [8]  [9]. 
Edney  investigated  the  impingement  problem  for  a  variety  of  model  geometries  in  a  blowdown  tunnel 
at  Sweden’s  Flygtekniska  Forsbksanstalten  (FFA).  The  FFA  tunnel  provided  Mach  4.6  and  7.0  flows 
of  air  with  maximum  stagnation  conditions  of  625  K  and  1.7  MPa  through  a  200  mm  diameter 
test  section.  This  wind  tunnel  was  equipped  with  a  model  injection  system  that  allowed  the  model 
to  be  removed  from  the  test  section  during  the  startup  process.  Edney’s  investigation  arose  out 
of  the  need  to  calibrate  thin  film  heat  transfer  gauges.  Anomalous,  elevated  heat  transfer  rates 
were  observed  as  the  calibration  model  was  traversed  through  a  spurious  shock  wave  generated  at 
the  tunnel  wall.  Subsequent  high  quality  schlieren  cinematography  of  the  controlled  injection  of 
models  through  shocks  of  known  strength,  combined  with  novel  heat  transfer  and  surface  pressure 
measurement  techniques,  allowed  Edney  to  effectively  classify  the  spectrum  of  interaction  regimes. 

Consider  the  impingement  of  a  shock  wave  on  the  leading  edge  of  a  blunt  body  as  illustrated 
schematically  in  figure  1.3.  For  the  compressible  flow  of  a  perfect  gas  a  sufficient  set  of  dimensionless 


Figure  1.1:  Schlieren  visualization  of  the  flow  field  about  a  model  of  the  X-15-2  vehicle  with  inter¬ 
ference  from  a  scram-jet  mounted  on  the  vertical  stabilizer  (Burcham  &  Nugent  [2]). 


Figure  1.2:  Damage  to  X-15  vertical  stabilizer  caused  by  impingement  of  scram-jet  inlet  shock  wave 
(Edney  [8]). 
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Figure  1.3:  Schematic  representation  of  the  impingement  of  an  oblique  shock  wave  on  the  leading 
edge  of  a  blunt  body. 

parameters  to  describe  any  quantity,  in  the  inviscid  flow  away  from  the  vehicle  surface  is, 


(j)  = 

where  Moo  is  the  free-stream  Mach  number,  7  is  the  ratio  of  specific  heats,  Pi  is  the  angle  of  the 
impinging  shock  wave,  g  describes  the  position  of  the  impinging  shock  wave  with  respect  to  the  body 
and  r  is  a  set  of  parameters  defining  the  body  geometry.  Although  viscous  effects  are  important  at 
the  surface  of  the  body  where  the  heating  occurs,  the  fundamental  mechanisms  that  control  the  flow 
will  be  seen  to  be  of  an  inviscid  nature.  In  the  case  of  fixed  free-stream  conditions,  impinging  shock 
strength  and  body  geometry  the  only  remaining  dependence  is  the  location  of  the  impingement  point 
relative  to  the  body;  i.e.  (j)  =  (j)[g].  As  the  incident  shock  wave  is  translated  relative  to  the  body  it 
potentially  intersects  with  all  possible  shock  strengths  of  both  positive  and  negative  slope.  On  the 
basis  of  an  experiment  such  as  this,  with  spherical  and  other  axisymmetric  blunt  bodies,  Edney  [8] 
[9]  observed  and  categorized  six  distinct  interaction  regimes.  The  interactions,  that  are  known  as 
types  I- VI,  occur  for  different  locations  of  the  impinging  shock  wave  with  respect  to  the  body  as 
indicated  in  figure  1.3.  A  review  of  this  method  of  classification  is  contained  in  appendix  A. 

Severe  heating  occurs  when  the  incident  shock  wave  impinges  in  the  vicinity  of  the  geometrical 
stagnation  point  and  a  type  IV  flow  results  (Figures  1.4  &  1.5).  Our  discussion  will  concentrate 
initially  on  the  type  IV  flow  since  the  heating  rates  for  the  remaining  regimes  are  less  intense.  A  three 
shock  A-pattern  is  observed  at  the  point  where  the  impinging  shock  wave  interacts  with  the  bow 
shock.  The  interaction  produces  a  strong  vortex  sheet  that  separates  the  upper  region  of  subsonic 
flow  from  the  lower  region  of  supersonic  flow.  A  second  inverted  A-pattern  is  observed  to  appear 
where  the  strong  bow  shock  continues  below  the  interaction  region.  Embedding  of  the  supersonic 
region  between  the  two  A-shocks  in  the  surrounding  subsonic  shock  layer  causes  the  formation  of  a 


Figure  1.4:  Holographic  shadowgraph  of  the  impingement  of  an  oblique  shock  wave  near  the  geomet¬ 
rical  stagnation  point  of  a  cylinder  in  hypersonic  flow.  Following  Edney  [8],  this  flow  is  categorized 
as  type  IV. 


supersonic  jet.  Greatly  elevated  heat  transfer  rates  and  surface  pressures  are  observed  where  the  jet 
impinges  on  the  body. 

Figure  1.5  shows  the  model  of  the  type  IV  flow  field  that  was  originally  proposed  by  Edney  [8]. 
Whilst  the  global  flow  field  is  quite  complex,  some  insight  may  be  obtained  by  mapping  the  flow 
field  in  the  vicinity  of  the  shock  wave  intersection  points  into  the  pressure-flow  deflection  angle  {p-6) 
plane  (Figure  1.6).  In  this  method  the  shock  interaction  solution  is  represented  graphically  in  the 
p-S  plane  by  the  intersection  of  the  loci  of  possible  downstream  states  for  each  of  the  interacting 
waves.  Further  details  of  the  technique  may  be  found  in  appendix  A.  Here  it  is  sufficient  to  recognize 
that  given  only  the  free  stream  conditions  and  the  incident  shock  angle,  the  intersection  of  the  loci 
representing  states  2  and  3  in  figure  1.6  completely  determines  the  wave  angles  and  flow  properties 
at  the  shock  impingement  point.  Since  this  fixes  the  strength  of  the  wave  that  connects  the  two 
A-points,  the  solution  for  states  4  and  5  at  the  second  inverted  A-point  is  also  determined.  Note  that 
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Figure  1.5:  Schematic  of  the  type  IV  shock  interference  flow  field. 


the  formation  of  the  jet  and  the  state  of  the  gas  in  it  are  therefore  independent  of  the  interaction 
of  the  flow  with  the  body.  This  represents  the  fundamental  inviscid  mechanism  that  controls  the 
shock  impingement  heating  phenomenon. 

Consider  next  the  problem  of  estimating  the  pressure  loads  and  heat  fluxes  at  the  jet  impingement 
point.  In  order  to  estimate  these  quantities  a  local  model  is  formulated  under  the  assumption  that 
the  oblique  waves  in  the  supersonic  jet  remain  straight.  This  is  justified  on  dimensional  grounds. 
The  Mach  number  of  the  jet  is  high  so  that  the  inertia  of  the  supersonic  flow  dominates  the  pressure 
forces  of  the  subsonic  shock  layer.  Similarly,  the  pressures  in  the  subsonic  shock  layer  above  and 
below  the  supersonic  jet  are  assumed  to  remain  constant.  Propagation  of  the  oblique  wave,  5, 
along  the  supersonic  jet  proceeds  by  multiple  reflections  from  the  shear  layers  that  are  treated  as 
free  streamline  boundaries.  The  oblique  wave,  5,  is  initially  reflected  from  the  shear  layer  as  an 
expansion  fan  so  that  the  constant  pressure  boundary  condition  is  satisfied.  This  expansion  fan,  and 
all  subsequent  reflected  waves,  are  assumed  to  be  isentropic  in  order  to  compute  the  state  upstream 
of  the  terminating  strong  shock.  Specification  of  the  pressure  and  entropy  completely  defines  this 
state  however  the  limitations  of  local  modeling  prevent  exact  determination  of  the  pressure.  Since 
the  solution  for  the  type  IV  flow  is  based  on  local  models  that  apply  only  in  the  vicinity  of  the 
shock  interactions,  the  scale  of  the  impingement  flow  field,  the  location  of  the  terminating  shock 
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Figure  1.6:  Representation  of  the  type  IV  shock  interference  flow  field  in  the  p-6  plane. 

and  hence  the  exact  state  of  the  gas  upstream  of  it  remain  indeterminate.  The  unknown  length 
scales  are  determined  by  the  interaction  of  the  jet  with  the  body.  Without  further  knowledge  of  the 
number  of  oblique  wave  reflections  in  the  supersonic  jet,  it  is  only  possible  to  bound  the  pressure 
upstream  of  the  terminating  shock  by  the  pressures  in  the  subsonic  regions  above  and  below  the  jet. 
For  simplicity,  the  pressure  corresponding  to  state  5  is  assumed  for  all  calculations  in  the  current 
study. 

A  portion  of  the  supersonic  jet  passes  both  above  and  below  the  cylinder  and  so  the  stagnation 
streamline  must  pass  through  the  supersonic  jet  (Figure  1.5).  Computation  of  the  pressure  at  the 
jet  impingement  point  is  based  on  the  further  assumption  that  the  terminating  shock  is  normal  to 
the  distorted  stagnation  streamline.  Whilst  stagnation  enthalpy  is  conserved  along  all  streamlines 
in  steady  flow,  the  stagnation  density  is  higher  for  streamlines  that  pass  through  the  weak  jet  shock 
system  than  for  streamlines  that  experience  a  larger  entropy  rise  across  the  adjacent  strong  bow 
shock.  The  elevated  density  and  the  strong  velocity  gradients  produced  by  the  impingement  of 
the  supersonic  jet  provide  the  mechanism  for  locally  increased  heat  transfer.  Edney  [8]  used  the 
stagnation  point  boundary  layer  solution  of  Cohen  &  Reshotko  [5]  to  correlate  the  heat  transfer  rate 
at  the  impingement  point  with  the  peak  pressure  and  body  diameter  to  jet  width  ratio,  D/Ax; 


Here  the  peak  heat  transfer  rate  and  pressure  are  normalized  with  respect  to  the  undisturbed  stag- 
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Figure  1.7:  Ratio  of  jet  impingement  pressure  to  undisturbed  stagnation  point  pres¬ 
sure  AS  A  function  of  INCIDENT  SHOCK  FLOW  DEFLECTION  ANGLE,  6i.  —  •  —  CUrves  are  for 

the  perfect  gas  model  with  Moo  =  10  (after  Edney  [8];  Figure  7-9). - curve  is  for  the  IDG  model 

in  the  frozen  limit  with  Poo  =  0.0075  and  aoo  =  0. 

nation  point  values.  The  pressures  that  appear  in  equation  1.1  are  determined  from  the  local  type 
IV  solution  in  the  p-6  plane. 

Curvature  of  the  supersonic  jet  is  observed  in  experimental  investigations  of  the  type  IV  flow. 
We  have  noted  that  the  strength  of  the  oblique  wave,  5,  is  determined  entirely  by  the  free  stream 
conditions  and  the  strength  of  the  incident  shock.  The  strength  of  this  wave  controls  the  pressure 
difference  that  exists  between  the  subsonic  regions  above  and  below  the  interaction.  It  follows  that 
the  pressure  gradient  across  the  jet,  and  hence  its  curvature,  vary  inversely  with  the  width  of  the 
jet. 

The  preceding  discussion  summarizes  the  mechanisms  that  control  the  shock  impingement  heat¬ 
ing  phenomenon.  The  key  conclusion  is  that  a  non-linear  interaction  of  the  shock  waves  at  the 
impingement  point  leads  to  the  formation  of  a  supersonic  jet  that  is  embedded  in  the  surrounding 
subsonic  shock  layer.  This  inviscid  jet  formation  process  is  local  in  nature  and  is  independent  of  the 
interaction  of  the  jet  with  the  body.  Elevated  pressure  loads  and  heat  transfer  rates  are  produced 
by  the  impingement  of  the  low  entropy  jet  fluid  on  the  surface  of  the  body. 


■  /  /> 
/  //' 

.  / 

/  ^ 


'y  7-1.2 

\ 

\ 

\ 

\ 

iDG  frozen;  7=1.33 

7=1 .4  ^ 

7=1.67 


8 


1.3  Real  gas  effects 


Dimensional  analysis  of  the  inviscid  reacting  flow  away  from  the  surface  of  a  hypervelocity  vehicle 
leads  to  three  parameters  involving  the  velocity; 


M  = 


«/ 


j.  mute 
k6d 


A  = 


(1.2) 


where  «oo  is  the  velocity  of  the  vehicle,  a;  is  the  frozen  speed  of  sound,  k  is  Boltzmann’s  constant,  6d 
is  a  temperature  characterizing  the  energy  of  the  dissociation  reaction,  ^  characterizes  the  reaction 
rate  and  D  is  some  characteristic  dimension  of  the  vehicle.  For  the  case  of  a  vehicle  traveling  at 
near  orbital  speeds  the  Mach  number,  M,  is  large  and  strong  shock  waves  will  be  dominant  features 
of  the  flow.  Oxygen  and  nitrogen  dissociation  occur  at  specific  kinetic  energies  of  approximately  17 
MJ /kg  and  34  MJ /kg  respectively  and  so  strong  dissociation  will  occur  in  the  flow  about  the  vehicle. 
This  influence  on  the  equilibrium  equation  of  state  is  expressed  by  the  dimensionless  specific  kinetic 
energy,  K.  The  dissociation  rate  parameter.  A,  which  is  equivalent  to  the  Damkohler  number  in 
related  chemical  kinetics  problems,  expresses  the  finite  length  over  which  the  reactions  occur  relative 
to  the  size  of  the  vehicle.  The  reaction  rate  parameter  will  increase  from  zero  for  a  vehicle  in  orbit  to 
very  large  values  at  low  altitudes.  Since  the  rates  of  dissociation  reactions  scale  with  the  density  of 
the  flow,  finite  rate  chemical  reactions  become  coupled  with  the  gas  dynamics  at  some  intermediate 
condition  where  A  1. 

By  mapping  the  type  IV  flow  field  into  the  p-8  plane,  as  discussed  in  section  1.2,  Edney  [8]  was 
able  to  compute  the  variation  of  the  pressure  at  the  jet  impingement  point  for  various  free-stream 
conditions.  Figure  1.7  shows  the  ratio  of  jet  impingement  pressure  to  undisturbed  stagnation  point 
pressure  as  a  function  of  the  flow  deflection  angle  across  the  incident  shock,  (5i,  and  the  ratio  of 
specific  heats,  7.  The  peak  pressure  initially  increases  with  impinging  wave  strength  before  reaching 
a  maximum  for  moderate  strength  waves.  Significant  real  gas  effects  are  suggested  by  the  increase 
in  the  peak  pressure  with  decreasing  ratio  of  specific  heats.  Edney  [8]  showed  that  the  peak  pressure 
also  increases  strongly  with  increasing  free-stream  Mach  number. 

Although  the  variable  7  model  provides  only  a  crude  estimate  of  the  influence  of  equilibrium  real 
gas  effects,  figure  1.7  provides  strong  motivation  for  further  study  of  the  problem  in  high  enthalpy 
test  facilities.  The  primary  effect  of  dissociation  reactions  on  blunt  body  flows  at  large  values  of  the 
flow  specific  kinetic  energy,  K,  is  to  increase  the  density  of  the  gas  in  the  shock  layer.  Since  shock 
standoff  distances  vary  inversely  with  the  shock  layer  density,  equilibrium  real  gas  effects  must  be 
expected  to  influence  the  length  scales  of  the  type  IV  interaction  flow  field. 

In  the  nonequilibrium  regime,  where  A  ~  1,  a  thermochemical  length  scale  is  introduced  into  the 
problem.  Nonequilibrium  phenomena  are  most  pronounced  in  regions  of  a  flow  field  where  discontin- 
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uous  changes  in  the  thermodynamic  variables  produce  large  deviations  from  chemical  equilibrium. 
A  familiar  example  is  the  relaxation  that  occurs  downstream  of  the  bow  shock  ahead  of  a  blunt  body 
(Hornung  [19]).  The  disparate  shock  strengths  produced  at  the  two  mutually  inverted  A-shocks  in 
the  type  IV  flow  must  therefore  be  expected  to  produce  strong  nonequilibrium  effects.  Figure  1.8 
indicates  the  possible  real  gas  influences  on  a  A-shock  intersection  point.  The  three  shock  waves  pro¬ 
duce  widely  differing  dissociation  levels  and  reaction  lengths.  Variations  in  the  relative  magnitudes 
of  the  thermochemical  and  fluid  mechanical  length  scales  introduce  effects  that  are  not  encompassed 
by  existing  models. 

Figure  1.8  illustrates  an  additional  effect  on  the  shear  layers  that  are  generated  at  the  A-points. 
Energy  that  is  released  by  mixing  of  the  dissociated  low  speed  fluid  with  the  lower  temperature 
supersonic  fluid  influences  the  density  of  the  shear  layer  mixture.  Although  this  effect  is  not  addressed 
by  the  current  study,  turbulent  diffusion  must  become  significant  when  the  jet  width  is  small  relative 
to  the  shock  standoff  distance. 

Despite  the  inviscid  nature  of  the  mechanisms  that  control  the  flow,  real  gas  effects  also  influence 
the  viscous  flow  in  the  boundary  layer  that  forms  at  the  jet  impingement  point.  Recombination 
occurs  in  the  boundary  layer  because  the  wall  temperature  is  low  and,  in  the  immediate  vicinity  of 
the  jet  impingement  point,  the  behavior  parallels  that  of  the  classical  blunt  body  problem.  Fay  & 
Riddell  [11]  showed  that  when  the  recombination  rate  in  the  boundary  layer  is  large, 


Here  q  is  the  heat  transfer  per  unit  area  and  time,  pe^  p>e  and  are  the  density,  viscosity  and 

transverse  velocity  gradient  respectively  at  the  outer  edge  of  the  boundary  layer.  Jet  impingement 
does  not  influence  the  total  enthalpy  at  the  stagnation  point,  /iq,  since  stagnation  enthalpy  is  con¬ 
served  at  all  points  in  steady  flow.  This  also  implies  a  limited  influence  on  the  viscosity,  ^e,  since 
for  diatomic  gases  pe  ^  ^  The  most  important  real  gas  influences  on  the  stagnation  point 

flow  are  therefore  the  increased  density  of  the  outer  flow,  pej  Q^nd  the  coupled  effect  on  the  velocity 
gradient.  Continuity  implicitly  couples  the  density  and  velocity  gradient  at  the  edge  of  the  boundary 
layer.  Finally,  note  that  shock  impingement  greatly  increases  the  recombination  rates  in  the  bound¬ 
ary  layer  since  these  scale  with  the  square  of  the  density  intensification  at  the  jet  impingement 
point. 

We  conclude  that  both  equilibrium  and  finite  rate  real  gas  effects  may  be  expected  to  exert  a 
significant  influence  on  the  type  IV  interaction  flow  field.  Shock  impingement  was  shown  to  produce 
elevated  pressures,  and  hence  densities,  at  the  jet  impingement  point.  A  broad  range  of  possible 
influences  arise  as  a  consequence  of  the  strong  interdependence  of  real  gas  effects  and  density. 
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Figure  1.8:  Real  gas  effects  at  a  A-shock  intersection  point. 

1.4  Literature  review 

Before  proceeding  we  wish  to  note  additional  important  contributions  to  the  literature.  The  need  for 
an  improved  understanding  of  the  shock  impingement  phenomenon  motivated  some  of  the  earliest 
investigations  that  were  summarized  by  Keyes  &  Hains  [24].  Early  investigations  provided  flow 
visualization  and  surface  measurements  for  a  variety  of  configurations,  however  these  works  failed  to 
elucidate  the  fundamental  mechanisms  of  the  problem.  The  utility  of  the  data  that  were  obtained  was 
limited  by  the  large  number  of  parameters  involved  and  restricted  spatial  and  temporal  resolution  of 
the  sensors.  Edney  [8]  discussed  early  attempts  at  modeling  such  as  the  acoustic  disturbance  model 
of  Fontenot  [12],  the  vorticity  model  of  Francis  [13]  and  the  boundary  layer  model  due  to  Hiers  & 
Loubsky  [18].  These  were  based  on  limited  observations  and  the  physical  models  that  were  utilized 
now  appear  to  be  inadequate.  Keyes  &  Hains  [24]  reported  experiments  that  further  advanced  the 
work  of  Edney  [8]  and  conducted  a  parametric  study  of  the  heat  transfer  rates  predicted  by  Edney’s 
models. 

Hiers  &  Loubsky  [18]  conducted  experiments  in  a  shock  tunnel  at  enthalpies  sufficient  to  produce 
vibrational  excitation  and  limited  oxygen  dissociation  in  air  however  the  results  did  not  demonstrate 
any  conclusive  trends  with  enthalpy.  Wieting  &  Holden  [48]  have  conducted  a  recent  experimental 
study  and  obtained  data  on  the  influence  of  Reynolds  number,  Mach  number  and  shock  strength 
variations  at  modest  enthalpies.  Kortz  [27]  and  Kortz,  McIntyre  k  Eitelberg  [28]  report  on  a  recent 
experimental  investigation  conducted  at  enthalpies  sufficient  to  cause  significant  nitrogen  dissociation 
in  the  DLR  HEG  free  piston  shock  tunnel.  This  work  included  the  first  quantitative  interferometric 
visualization  of  shock  impingement  flows.  Elevated  heat  transfer  rates  were  not  observed  for  the 
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range  of  parameters  that  were  investigated  and  this  appears  to  have  resulted  from  a  relatively  coarse 
variation  of  the  shock  impingement  location. 

Progress  in  the  development  of  monotonicity  preserving  shock  capturing  algorithms  is  evident  in 
recently  reported  numerical  studies  of  the  shock  impingement  problem.  Shock  interference  flow  fields 
are  simulated  by  the  provision  of  an  inflow  boundary  condition  ahead  of  the  blunt  body  that  accounts 
for  the  discontinuity  across  the  impinging  shock.  Tannehill  &  Holst  [43]  produced  one  of  the  earliest 
attempts  at  such  a  simulation  however  the  results  exhibited  excessive  numerical  diffusivity.  Recent 
applications  of  monotonicity  preserving  algorithms  such  as  those  of  Klopfer  fc  Yee  [25]  succeeded  in 
capturing  the  gross  features  of  the  mean  flow  field.  The  unstructured  adaptive  meshing  techniques 
employed  by  Vemaganti  [44]  provided  remarkable  resolution  of  the  internal  structure  of  the  impinging 
jet  for  the  type  IV  flow  field.  Vemaganti  [44]  used  turbulence  models  to  simulate  the  effect  of  the 
shear  layers  generated  at  the  shock  intersection  points  on  the  surface  heat  flux.  The  effect  of  finite 
rate  chemical  reactions  was  considered  by  Hannemann  et  al.  [16]  who  simulated  the  experimental 
results  reported  by  Kortz  et  al.  [27]  [28].  These  computations  indicated  the  persistence  of  the 
type  IV  behavior  for  shock  impingement  locations  intermediate  to  those  discussed  by  Kortz  et  al. . 
Kinetic  theory  based  simulations  of  the  shock  impingement  problem  that  were  reported  by  Carlson 
k  Wilmoth  [3]  represent  an  unusual  application  of  DSMC  methods.  Results  were  presented  for  the 
variation  of  species  concentrations  across  the  shock  layer  and  the  impinging  jet  for  a  non-catalytic 
wall.  An  additional  recent  contribution  is  that  of  Yeung  k  Hall  [49]  who  discussed  differences  in  the 
computed  flow  fields  about  spherical  and  cylindrical  models. 

Numerical  simulations  of  the  mean  flow  fail  to  reproduce  the  unsteady  flow  features  that  are 
observed  in  the  experimental  data  (Figure  1.4).  The  potential  instabilities  include  those  of  the  shear 
layers  generated  at  shock  intersection  points  and  oscillations  of  the  jet  structure  itself.  Edney  [8] 
noted  oscillations  that  were  observed  in  the  jet  impingement  study  of  Henderson  [17]  but  did  not 
observe  any  analogous  fluctuations  of  the  type  IV  jet.  Jet  oscillations  appeared  in  the  computations 
of  Prabhu  et  al.  [36]  and  Gaitonde  [15]  where  limit  cycle  oscillations  were  observed  in  the  com¬ 
putational  residual  as  the  mesh  was  adaptively  refined.  Gaitonde  [15]  addressed  the  issue  of  time 
accuracy  of  the  algorithms  that  were  used  and  concluded  that  the  oscillations  were  attributable  to 
the  physics  of  the  problem.  Similar  conclusions  may  be  drawn  from  the  results  of  Zhong  [50]  who 
also  obtained  time  accurate  solutions  for  the  jet  fluctuations  using  essential  non-oscillatory  (ENO) 
schemes.  The  possibility  of  coherent  fluctuations  of  the  impinging  jet  restricts  the  gains  attainable 
by  adaptive  meshing  if  extremely  high  resolution  is  required  since  the  regions  of  grid  refinement  be¬ 
come  dynamic.  Despite  these  difficulties  the  numerical  solutions  are  of  great  pragmatic  value  since 
they  appear  capable  of  providing  the  mean  flow  data  necessary  for  the  application  of  correlations 
for  the  sub-grid  behavior. 

The  paucity  of  the  existing  data  in  the  hyperveiocity  regime  and  the  importance  of  the  phe- 
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nomenon  provide  strong  motivation  for  further  study  of  the  problem. 

1.5  Scope  of  the  current  study 

The  objective  of  the  present  work  is  to  determine  the  quantitative  effects  of  thermochemistry  on  the 
shock  impingement  phenomenon.  Models  that  were  advanced  by  Edney  [8]  suggest  that  increases 
of  the  Mach  number  and  reductions  in  the  ratio  of  specific  heats  should  produce  strong  increases  in 
the  peak  heat  transfer  rates.  Systematic  experimental  or  parametric  numerical  studies  to  investigate 
these  predictions  have  yet  to  appear.  The  current  work  consists  primarily  of  an  experimental  inves¬ 
tigation  of  the  shock  interaction  problem  that  utilizes  the  capabilities  the  GALCIT  T5  hypervelocity 
shock  tunnel. 

The  variable  7  model  used  by  Edney  [8]  represents  only  a  crude  estimate  of  the  influence  of  flow 
thermochemistry.  The  foregoing  survey  identifies  a  number  of  potentially  significant  mechanisms  by 
which  equilibrium  and  finite  rate  real  gas  effects  may  influence  the  shock  interaction  phenomenon. 
This  work  aims  to  develop  models  of  these  mechanisms  and  provide  quantitative  predictions  of  the 
influences  of  the  controlling  parameters. 

In  that  the  work  of  Edney  [8]  was  made  possible  by  effective  instrumentation,  a  significant 
preliminary  goal  of  the  current  study  was  to  develop  effective  diagnostic  techniques  for  use  in  the 
extreme  conditions  that  occur  in  high  enthalpy  test  facilities.  This  included  the  development  of 
surface  heat  transfer  gauges,  differential  and  holographic  interferometry,  high  speed  data  acquisition 
electronics  and  novel  flow  models.  These  topics  are  discussed  further  in  the  appendices. 
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Chapter  2  Experimental  program 


2.1  Description  of  the  T5  hypervelocity  shock  tunnel 


The  experimental  program  was  based  upon  the  capabilities  of  the  T5  free  piston  driver  shock  tunnel 
that  has  recently  been  constructed  at  GALCIT.  Dimensional  analysis  of  the  inviscid  flow  away  from 
the  surface  of  a  hypervelocity  vehicle  leads  to  three  parameters  (Equation  1.2)  involving  the  velocity; 


M=  — 
«/ 


_  mugo 

k0d 


ioL  n 


(2.1) 


These  dimensionless  parameters  express  the  requirements  for  similitude  in  a  ground  testing  facility. 
With  the  same  test  gas  the  velocities  of  the  full-scale  and  sub-scale  flows  must  be  matched  in  order 
to  achieve  the  correct  scaling,  K,  of  the  flow  specific  kinetic  energy  with  respect  to  the  dissociation 
energy.  This  requirement  for  high  velocities  implies  that  the  facility  stagnation  temperature  must  be 
large;  typically  10000  K  to  achieve  enthalpies  equivalent  to  low  earth  orbital  velocities.  High  densities 
are  required  for  similarity  of  the  reaction  rate  parameter,  A,  with  sub-scale  models  since  the  rates 
of  dissociation  reactions  are  proportional  to  the  density.  For  typical  vehicle  flight  trajectories  this 
implies  a  nozzle  reservoir  pressure  approaching  100  MPa. 

The  development  of  free  piston  driver  shock  tunnels,  by  Stalker  [41]  and  coworkers,  has  enabled 
the  simulation  of  these  dimensionless  parameters  in  the  laboratory.  Shock  tunnels  of  this  type  are 
unique  in  their  ability  to  produce  the  high  densities  necessary  for  simulation  of  both  the  kinetic 
energy,  K,  and  the  reaction  rate  parameter,  A.  The  T5  facility  is  illustrated  schematically  in 
figure  2.1.  High  enthalpy  and  high  pressure  stagnation  conditions  are  achieved  in  the  nozzle  reservoir 
by  shock  heating  the  test  gas  as  in  a  conventional  reflected  mode  shock  tunnel.  In  order  to  achieve 
a  sufficiently  high  shock  speed  at  high  pressures,  the  speed  of  sound  of  the  driver  gas  is  raised 
by  adiabatically  compressing  helium  to  100  MPa  and  3000  K  with  a  single  stroke  of  a  reusable 
compressed  air  driven  piston.  Further  details  of  the  performance  and  construction  of  this  complex 
facility  are  given  by  Hornung,  Sturtevant,  Belanger,  Sanderson,  Brouillette  &  Jenkins  [23]. 


2.2  Test  conditions 

Three  different  test  conditions  were  used  and  these  are  summarized  in  table  2.1.  The  first  case,  A, 
represents  a  low  enthalpy  control  condition  where  the  only  real  gas  effects  were  a  partial  excitation 
of  the  vibrational  modes  in  the  stagnation  regions  of  the  flow.  Case  C  represents  a  high  enthalpy 
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Figure  2.1:  Schematic  representation  of  the  T5  free  piston  driver  shock  tunnel. 


condition  that  was  chosen  to  highlight  the  influence  of  nonequilibrium  thermochemistry  on  the 
problem.  Case  B  is  an  intermediate  enthalpy  condition  that  was  designed  to  aid  interpretation 
of  the  dilferences  in  the  flow  mechanisms  for  cases  A  and  C.  Nitrogen  test  gas  was  used  for  all 
experiments.  Air  was  considered  too  aggressive  for  prolonged  model  life  with  shock  impingement. 

Appendix  B  contains  an  assessment  of  the  flow  quality  and  the  accuracy  of  the  test  condition 
predictions.  This  assessment  comprises  two-dimensional  nonequilibrium  nozzle  flow  computations, 
interferometric  stagnation  density  measurements  using  Abel  transform  methods  to  eliminate  end 
effect  errors,  reacting  flow  shock  standoff  distance  correlations,  stagnation  point  heat  flux  measure¬ 
ments  and  frozen  oblique  shock  density  comparisons.  The  conclusion  of  these  comparisons  was  that 
the  Mach  number  is  over-predicted  for  test  condition  A.  This  was  attributed  to  nozzle  flow  nonuni¬ 
formity  at  low  enthalpy  off  design  conditions.  Condition  B  was  very  accurately  confirmed  by  all 
the  experimental  measurements.  Condition  C  produced  blunt  body  shock  layer  densities  that  were 
somewhat  less  than  the  predictions.  Contamination  of  the  test  gas  with  the  helium-argon  driver 
gas  mixture  is  known  to  become  limiting  in  shock  tunnels  at  high  enthalpy  conditions.  It  was  not 
possible  to  measure  this  contamination  directly  with  the  apparatus  used  in  the  present  experiments. 
Although  indirect  measurements  that  are  discussed  in  appendix  B  demonstrate  some  anomalies, 
the  data  reported  here  remained  temporally  constant  and  self  consistent  early  during  the  test  time. 
These  data  agree  tolerably  well  with  established  theories.  Flow  contamination  that  was  detectable 
with  the  current  instrumentation  did  not  become  apparent  until  times  beyond  the  period  of  test 
flow  that  was  used  to  obtain  the  data. 

Differing  helium-argon  driver  gas  mixtures  were  used  at  the  three  enthalpies  to  achieve  tailored 
interface  operation.  In  this  condition  the  nozzle  reservoir  pressure  and  stagnation  point  heat  flux 
remained  temporally  constant  to  within  ±10%  during  the  test  time.  Uniform  conditions  were  main¬ 
tained  for  1.25  ms  at  condition  C  and  in  excess  of  5  ms  at  condition  A. 
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Test  condition 

A 

B 

C 

Shock  tube 
conditions 

Fill  pressure 
(kPa  ±0.5kPa) 

75 

32.5 

20 

Fill  temperature 
(K  ±4K) 

297 

297 

297 

Shock  speed 

1890 

3500 

4360 

(m/s) 

±100 

±245 

±255 

Nozzle 

Pressure 

14.6 

22.3 

28.3 

reservoir 

(MPa) 

±0.7 

il.3 

±1.9 

conditions 

Temperature 

(K) 

3210 

7550 

8960 

Enthalpy 

(MJ/kg) 

3.88 

12.0 

19.1 

Test  section 

Velocity  (m/s) 

2540 

4450 

5350 

conditions 

Density  (kg/m^) 

0.0218 

0.0155 

0.0157 

Pressure  (kPa) 

1.03 

5.48 

11.4 

N  concentration  (mole/kg) 

9.9  X  10“® 

6.9  X  10-1 

3.65  X  10° 

Test  section 

IDG  Pco 

7.32  X  10-^ 

1.79  X  10-2 

2.54  X  10-2 

dimensionless 

IDG  Ho^ 

0.102 

0.346 

0.560 

parameters 

IDG  Oioo 

1.4  X  10"® 

9.7  X  10-2 

0.051 

IDG  pd 

6.0  X  10® 

8.4  X  10® 

8.3  X  10® 

IDG  A  per  mm 

4.  X  10-12 

3.4  X  10-2 

0.525 

Re  per  mm 

5540 

1630 

1350 

Table  2.1:  SUMMARY  OF  THE  FREE-STREAM  CONDITIONS.  Test  conditions  were  computed  on  the 
basis  of  measurements  of  the  initial  shock  tube  fill  pressure  and  temperature,  incident  shock  speed 
prior  to  reflection,  and  nozzle  reservoir  pressure  after  shock  reflection.  The  dimensionless  parameters 
for  the  ideal  dissociating  gas  (IDG)  model  are  defined  in  section  4.2.  Error  estimates  given  for  the  fill 
conditions  represent  the  accuracy  of  the  pressure  gauge  and  the  variation  of  the  ambient  temperature. 
Error  estimates  for  the  measured  shock  speed  and  reservoir  pressure  are  the  standard  deviation  of 
the  quantities  sampled  over  the  entire  sequence  of  shots.  The  error  in  the  computed  quantities  may 
be  inferred  from  these  estimates. 
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Figure  2.2:  Abstracted  model  geometries  for  the  shock  impingement  problem. 

2.3  Description  of  experiment 

2.3.1  Abstraction  of  the  problem 

If  a  blunt  body  is  characterized  by  two  principal  radii  of  curvature  at  the  stagnation  point,  Ri  and 
i?2,  then  it  is  reasonable  to  consider  the  following  simplified  cases  (Figure  2.2)  as  abstractions  of  the 
complex  shock  impingement  flows  that  occur  about  flight  vehicles; 

1.  Spherical  body  with  =  i22- 

2.  Cylindrical  body,  with  Ri-^oo,  oriented  normal  to  the  free-stream  flow  direction  and  parallel 
to  the  plane  of  the  incident  shock  wave.  This  is  a  two-dimensional  flow  field. 

3.  Swept  cylindrical  fin,  with  i?2  —  oo,  lying  in  the  plane  defined  by  the  free-stream  flow  direction 
and  the  normal  to  the  incident  shock  surface. 

Each  of  these  cases  exhibits  a  single  geometrical  degree  of  freedom  in  the  location  of  the  impinging 
shock  wave.  In  cases  1  and  2  this  is  the  vertical  displacement  of  the  incident  wave  with  respect  to 
the  model.  If  the  fin  of  case  3  is  so  long  that  the  flow  is  cylindrically  symmetric  far  away  from  the 
ends,  the  sole  geometric  dependence  of  the  flow  is  the  sweep  angle  of  the  fin. 

The  experiments  described  in  this  report  were  intended  to  test  models  of  the  influence  of  thermo¬ 
chemistry  on  the  local  behavior  near  the  impingement  point  and  the  manner  in  which  this  influences 
the  global  flow  about  the  body.  The  local  models  are  rigorous  in  nature  whereas  the  global  flow 
models  represent  interpretations  of  the  flow  structure.  Whilst  these  topological  considerations  are 
qualitative  in  nature,  the  streamline  patterns  that  are  used  to  describe  the  shock  interaction  flow 
fields  must  be  kinematically  admissible.  This  aspect  has  been  somewhat  neglected  in  previous  studies 
of  the  phenomenon. 
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2.3.2  Discussion  of  the  flow  topology 

Consider  the  model  of  the  type  III  flow  field  that  is  shown  in  figure  A. 5.  On  the  basis  of  an  analogy 
with  the  flow  over  a  backward  facing  step,  heating  was  presumed  to  occur  at  the  point  where  the 
shear  layer  attaches  to  the  body  (Edney  [8]).  The  two  flows  are  topologically  equivalent  only  in 
the  extraordinary  case  where  the  stagnation  streamline  coincides  with  the  shear  layer.  For  all  other 
cases  the  analogy  violates  the  continuity  equation  in  the  absence  of  the  upstream  separation  event 
that  occurs  in  the  flow  over  a  step.  A  large  mass  flux  enters  the  shock  layer  across  the  bow  wave 
whereas  the  mass  contained  in  the  recirculation  region  downstream  of  a  backward  facing  step  is 
constant. 

Similar  difficulties  arise  in  the  interpretation  of  the  type  IV  model  as  the  incident  shock  wave 
moves  upward.  This  causes  the  width  of  the  supersonic  jet  to  decrease  and  its  curvature  to  increase. 
Impingement  of  the  type  IV  jet  requires  that  the  stagnation  streamline  pass  through  it.  Since  the 
width  of  the  jet  is  small,  this  can  occur  for  only  a  limited  range  of  locations  of  the  impinging  shock 
wave.  Again,  continuity  prohibits  the  impingement  of  the  jet  if  a  stagnation  point  is  presumed  to 
exist  at  some  other  point  on  the  body. 

Despite  the  different  characteristics  of  the  type  III  and  IV  flows,  the  topological  concerns  are 
identical  in  both  cases.  We  conclude  that  attachment  of  the  disturbance  generated  by  shock  impinge¬ 
ment  is  possible  only  in  the  seemingly  extraordinary  case  when  the  type  III  shear  layer  coincides  with 
the  stagnation  streamline  and  for  the  limited  cases  when  the  stagnation  streamline  passes  through 
the  type  IV  jet. 

2.3.3  Critical  points  in  the  streamline  patterns 

The  anomalies  that  are  noted  above  appear  to  arise  from  a  two-dimensional  interpretation  of  line 
of  sight  integrated  images  of  a  three-dimensional  flow  field.  From  these  images  the  shear  layer  or 
type  IV  jet  away  from  the  axis  of  symmetry  appears  to  meet  the  surface  at  a  finite  angle  as  it  passes 
around  the  three-dimensional  model.  There  exist  alternate  three-dimensional  surface  flow  patterns 
that  can  reconcile  the  difficulty. 

All  three  abstracted  configurations  that  are  considered  here  exhibit  a  plane  of  symmetry.  Stream¬ 
lines  are  necessarily  contained  in  this  plane  in  which  the  out-of-plane  velocity  is  zero.  Streamlines 
are  also  locally  contained  in  the  plane  tangential  to  surface  at  any  given  point  on  the  body.  In  the 
vicinity  of  a  critical  point,  where  the  in-plane  velocities  vanish  and  the  streamline  direction  is  inde¬ 
terminate,  the  kinematical  problem  in  planes  that  contain  streamlines  reduces  to  a  pair  of  coupled 
autonomous  ordinary  differential  equations; 


7-)  I  2\ 

Ui  =  —  =  BijXj  +  0{xj) 


i=l,2  i=l,2. 


(2.2) 
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Figure  2.3:  Surface  streamlines  for  shear  layer /jet  attachment  in  the  flow  over  a  spherical  body. 

The  coefficients,  Bij,  are  determined  by  substituting  the  expansion  into  the  equations  of  motion. 
Here  we  consider  only  inviscid  streamline  patterns  so  that  higher  order  terms  of  the  expansion  may 
be  neglected.  For  the  purposes  of  the  immediate  discussion  it  is  sufficient  to  recognize  that  all 
possible  critical  points  in  the  streamline  patterns  are  described  by  the  node,  saddle,  focus  and  center 
type  solutions  of  equation  2.2.  The  form  of  each  critical  point  is  determined  by  the  eigenvalues  of 
Bij .  The  terminology  that  is  used  in  this  description  is  based  on  the  approach  of  Perry  &  Chong  [35]. 

Consider  first  the  impingement  of  an  oblique  shock  wave  onto  a  spherical  model.  The  apparent 
contradiction  of  continuity  is  resolved  if  a  saddle  point  is  assumed  to  form  in  the  surface  streamline 
pattern  between  the  nodes  at  the  stagnation  point  and  shear  layer  /  jet  attachment  point  (Figure 
2.3).  This  interpretation  is  consistent  with  the  surface  oil  flow  visualization  of  Keyes  &  Hains  [24] 
and  the  computations  of  Yeung  &;  Hall  [49]  who  observed  that  the  velocity  vanished  at  three  points 
along  the  centerline  for  a  spherical  body.  In  the  plane  of  symmetry,  the  stagnation  point  and  shear 
layer  /  jet  attachment  point  appear  as  saddle  points  (Figure  2.4).  The  intermediate  critical  point 
may  be  either  a  reattachment  node  or  a  separation  saddle.  The  local  buildup  of  oil  observed  near 
the  middle  critical  point  in  the  results  of  Keyes  k  Hains  [24]  supports  the  separation  hypothesis.  In 
this  case  a  focus  must  be  inserted  away  from  the  surface  of  the  body  to  match  the  flow  entering  the 
shock  layer  (Figure  2.4).  Streamlines  away  from  the  plane  of  symmetry  spiral  around  the  streamline 
that  emerges  out  of  the  plane  from  the  center  of  the  focus.  When  inverted,  the  streamline  pattern 
shown  in  figure  2.4  also  describes  the  impingement  of  the  type  IV  jet  above  the  stagnation  point.  We 
conclude  that  the  coexistence  of  a  stagnation  point  and  an  impingement  point  on  a  three-dimensional 
blunt  body  leads  directly  to  the  formation  of  a  saddle  point  in  the  surface  streamline  pattern  and  a 
spiral  vortex  that  wraps  around  the  body. 

Symmetry  prohibits  the  formation  of  a  saddle  in  the  surface  streamline  pattern  in  the  two- 
dimensional  shock  impingement  problem  and  so  only  a  single  critical  point  can  exist  on  the  forebody. 
The  important  conclusion  is  that  either  the  stagnation  point  and  the  type  III  shear  layer  attach¬ 
ment  /  type  IV  jet  impingement  point  must  coincide  exactly,  or  that  attachment  /  impingement 
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Figure  2.4:  Streamlines  in  the  symmetry  plane  showing  the  formation  of  a  spiral  vortex  as  a  conse¬ 
quence  of  shear  layer/jet  attachment. 


Figure  2.5:  Streamline  pattern  for  a  cylindrically  symmetric  swept  fin. 

does  not  occur.  If  the  type  III  shear  layer  or  type  IV  jet  disturbance  does  not  coincide  with  the 
stagnation  streamline  then  it  must  convect  past  the  surface  and  remain  at  a  finite  distance  from 
it  far  downstream.  This  is  in  direct  contrast  to  the  surface  streamline  pattern  for  the  undisturbed 
cylindrically  symmetric  flow  over  a  swept  fin  that  is  shown  in  figure  2.5.  The  pattern  follows  from 
the  addition  of  a  constant  transverse  velocity  component  to  the  two-dimensional  flow  field  about  a 
cylinder.  Observe  that  all  streamlines  in  the  symmetry  plane  reach  the  surface  where  they  arrive  and 
depart  tangentially  from  the  bifurcation  line  that  runs  the  length  of  the  cylinder.  The  conclusion 
is  that  disturbances  asymptote  to  the  surface  for  a  wide  range  of  impingement  locations  and  sweep 
angles.  Although  shock  impingement  disturbs  these  symmetries,  the  trends  provide  some  insight. 
Finite  length  cylinders  with  an  impinging  wave  are  strictly  topologically  equivalent  to  the  spherical 


case. 
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2.3.4  Application  of  holographic  interferometry 

Holographic  interferometry  provides  a  tool  for  directly  measuring  the  location  of  the  stagnation 
point  and  any  separation  regions  in  planar  flow.  Abel  transform  techniques  are  not  applicable 
in  the  presence  of  shock  impingement  since  the  flow  about  an  axisymmetric  body  becomes  fully 
three-dimensional.  Maxima  in  the  fringe  shift  of  the  interferogram  correspond  to  the  saddle  point 
and  center  type  features  that  are  admissible  for  two-dimensional  flow.  This  interpretation  remains 
approximately  valid  in  the  presence  of  shock  induced  vorticity  if  the  velocity  gradient  near  a  crit¬ 
ical  point  in  the  velocity  field  remains  large  with  respect  to  the  entropy  gradient  (appropriately 
normalized).  Species  concentration  gradients  also  affect  this  interpretation  in  reacting  flows.  The 
technique  remains  valid  when  the  maxima  of  density  and  dissociation  level  coincide  since  the  re- 
fractivity  of  nitrogen  atoms  is  larger  than  that  for  the  molecules.  Since  density  measurements  are 
sensitive  to  real  gas  effects,  holographic  interferometry  provides  quantitative  data  that  are  valuable 
for  discriminating  investigations  of  the  problem. 

2.3.5  Summary 

For  the  reasons  outlined  above,  the  majority  of  the  experiments  that  are  described  in  this  report 
have  concentrated  on  the  nominally  two-dimensional  mean  flow  associated  with  the  impingement 
of  an  oblique  shock  wave  on  the  leading  edge  of  a  coplanar  cylinder.  This  enabled  investigation  of 
the  simplest  possible  flow  topology  nnder  conditions  where  it  was  possible  to  directly  observe  the 
behavior. 


2.4  Apparatus 

2.4.1  Test  section  arrangement 

The  experimental  arrangement  is  shown  in  figure  2.6.  High  temperature  and  pressure  gas  was 
generated  in  the  nozzle  reservoir  by  the  T5  hypervelocity  free  piston  shock  tunnel.  The  gas  was 
expanded  into  the  test  section  using  a  contoured  axisymmetric  nozzle  with  an  area  area  ratio  of  109 
and  an  exit  diameter  of  315  mm.  The  model  was  positioned  in  the  free  jet  test  section  through  which 
the  test  gas  exhausted  into  an  evacuated  dump  tank.  A  large  shock  generator  was  partially  inserted 
into  the  uniform  nozzle  core  flow  that  extends  upstream  of  the  nozzle  exit  plane.  Variations  in  the 
shock  impingement  flow  field  were  investigated  by  vertically  translating  the  model  with  respect  to 
the  exit  of  the  shock  generator.  Care  was  taken  to  ensure  that  the  nozzle  exit  expansion  fan  and 
separated  shock  generator  boundary  layer  did  not  corrupt  the  flow  about  the  model.  A  holographic 
interferometer  was  used  to  visualize  the  shock  waves  and  figure  2.7  shows  the  vertical  arrangement  of 
the  U-shaped  interferometer  enclosure  about  the  test  section.  Optical  access  was  available  through 
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Figure  2.6:  Layout  of  apparatus  in  the  test  section  of  T5. 


two  <1)200  mm  x40  mm  BK7  glass  windows  that  were  mounted  to  flanges  on  the  sides  of  the  test 
section.  Figure  2,6  is  an  accurate  scale  drawing  of  the  arrangement.  All  apparatus  was  mounted 
symmetrically  with  respect  to  the  centerline  in  the  third  dimension  that  is  not  shown. 

2.4.2  Alignment  accuracy 

The  model,  shock  generator  and  holographic  interferometer  were  aligned  with  respect  to  the  nozzle 
exit  within  ±lmrad  about  all  three  axes.  Beam  divergence  of  approximately  2  mrad  for  the  NdYAG 
laser  light  source  obviated  the  need  for  greater  alignment  accuracy.  Optical  alignment  was  performed 
by  reflecting  the  interferometer  beam  from  a  mirror  attached  to  a  flat  reference  surface  on  the  end 
of  the  cylindrical  model.  Care  was  necessary  here  to  ensure  that  flow  features  close  to  the  wall  of 
the  cylinder  could  be  accurately  resolved.  The  positional  accuracy  of  the  test  section  apparatus  was 
estimated  to  be  ±5  mm  in  all  three  coordinate  directions.  The  vertical  position  of  the  model,  which 
controlled  the  shock  impingement  location,  was  repeatable  with  a  relative  accuracy  of  =b0.04  mm. 

2.4.3  Description  of  model 

Large  diameter  models  were  desirable  in  order  to  resolve  the  fine  flow  structure  and  to  achieve  con¬ 
ditions  of  thermochemical  nonequilibrium  at  moderate  pressures.  Low  pressure  operation  produces 
excellent  flow  quality  in  the  T5  shock  tunnel  if  the  tube  and  test  section  are  cleaned  adequately 
between  shots.  A  cylinder  with  aspect  ratio  4.5  provided  the  largest  possible  model  that  adequately 
approximated  a  two-dimensional  flow  (Sykes  [42]).  Optimization  of  the  arrangement  within  the  geo¬ 
metrical  constraints  imposed  by  the  nozzle  diameter,  nozzle  exit  expansion  fan,  shock  generator  and 
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Figure  2.7:  Arrangement  of  the  optical  apparatus  with  respect  to  the  test  section. 


optical  access  led  to  the  design  indicated  in  figure  2.6  and  a  model  diameter  of  40  mm  (Figure  F.l). 
Xhe  model  support  shown  in  figure  F.2  allowed  precise  vertical  translation  of  the  model  to  study  the 
effect  of  different  shock  impingement  locations.  A  cavity,  internal  to  the  model  support,  contained 
preamplifiers  for  the  model  instrumentation.  This  cavity  and  all  components  of  the  model  were 
sealed  from  the  external  flow  with  0-rings  and  vented  to  the  atmosphere  through  pressure  tubing 
attached  to  the  rear  of  the  model  support.  This  unusual  feature  contributed  to  signal  quality  by 
isolating  the  instrumentation  from  the  test  section  gases.  Glands  allowed  the  pressure  tubing,  which 
also  served  as  an  instrumentation  wiring  conduit,  to  pass  through  the  wall  of  the  test  section. 


2.4.4  Instrumentation 

Sensors  were  available  to  monitor  all  of  the  operating  parameters  of  the  T5  facility  and  determine 
the  free-stream  conditions  in  the  test  section.  Piezo  ceramic  pressure  transducers  installed  along  the 
length  of  the  shock  tube  were  used  to  monitor  the  speed  of  the  incident  and  reflected  shock  waves 
and  the  pressure  in  the  nozzle  reservoir. 

Thermocouple  gauges  were  used  to  record  the  surface  temperature  at  24  circumferential  locations 
around  the  model  forebody.  The  gauges  were  staggered  in  three  planes  arranged  symmetrically  with 
respect  to  the  centerline  to  achieve  maximum  resolution  (Figure  F.l).  Additional  thermocouple  ports 
were  provided  at  the  rear  of  the  body  however  these  were  not  used  for  the  experiments  described 
here.  The  sensitive  area  of  the  gauges  subtended  an  angle  of  2"  at  the  center  of  the  cylinder  and 
were  spaced  at  5°  intervals  in  the  most  densely  instrumented  portion. 
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Differential  thermocouple  preamplifiers  (Figure  E.4)  were  housed  inside  the  test  section  to  provide 
noise  levels  of  -35  -55dB  depending  on  the  signal  level  and  gain  setting.  The  lower  figure  reflects 

the  need  to  provide  headroom  for  the  anticipated  intense  heat  transfer  rates.  Mean  square  optimal 
spectral  deconvolution  techniques  were  used  to  infer  time  resolved  heat  transfer  rates  from  the 
surface  temperature  data.  Details  are  given  in  appendix  E  along  with  a  description  of  the  novel 
gauge  design  that  was  used  to  ensure  adequate  high  frequency  response.  The  signals  were  digitized 
with  ±11  bit  resolution  and  recorded  at  a  sampling  frequency  of  200  kHz,  Aliasing  was  prevented 
by  the  bandwidth  of  the  analog  electronic  circuitry  and  principally  by  the  finite  sensor  size  which 
limited  the  scale  of  the  resolvable  unsteady  flow  structures. 

2.4.5  Flow  visualization  and  interferometry 

Laser  interferometry  was  used  to  visualize  the  flow  field  and  to  provide  quantitative  measurements 
of  the  density  field.  Preliminary  experiments  utilized  a  (^200mm  field  of  view  differential  interferom¬ 
eter  with  a  conventional  folded  Z  mirror  arrangement.  Continuously  variable  differential  shear  was 
provided  by  the  photoelastic  birefringence  of  a  polycarbonate  prismatic  beam. 

The  final  experiments  utilized  a  (/>75  mm  field  of  view  high  resolution  holographic  interferometer. 
Double  pulse  holography  was  used  to  sequentially  record  the  test  section  wavefront  before  and  during 
the  flow.  When  the  two  beams  were  simultaneously  reconstructed  an  interferogram  was  produced 
that  was  equivalent  to  an  infinite  fringe  Mach-Zehnder  interferogram.  Note  that  this  procedure,  as 
it  was  applied  to  the  experiments  described  here,  produces  no  three-dimensional  information.  The 
device  that  is  described  in  appendix  D  incorporated  a  multi-element  lens  imaging  arrangement  for  the 
control  of  optical  aberrations,  Fourier  plane  conjugate  image  separation,  white  light  reconstruction 
and  Gaussian  apodizing  filters  for  non-oscillatory  imaging  of  shock  fronts. 

2.4.6  Shock  generator 

A  significant  difficulty  in  this  work  was  the  generation  of  the  impinging  shock  wave.  At  moderate 
Mach  numbers  the  oblique  shock  wave  produced  by  a  wedge  already  lies  quite  close  to  the  surface. 
The  separation  of  the  wedge  surface  and  the  shock  wave  determines  the  permissible  model  diameter. 
The  shock  generator  that  is  illustrated  in  figure  2.8  represents  a  response  to  this  difficulty  and  the 
concept  is  analogous  to  a  wave-rider  in  hypersonic  vehicle  design  (Nonweiler  [33]).  Its  shape  is 
generated  by  inclining  a  cylindrical  tube  at  some  arbitrary  flow  deflection  angle  and  sectioning  the 
tube  at  the  corresponding  oblique  shock  angle.  This  defines  an  internal  surface  that  satisfies  the 
inviscid  boundary  condition  for  the  uniform  flow  behind  a  planar  oblique  shock  wave  attached  to 
the  elliptical  planform  leading  edge.  The  argument  holds  for  non-circular  cylindrical  sections  such 
as  the  familiar  example  of  a  flat  plate  with  side  fins.  Since  the  tube  has  finite  thickness  the  rear 
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Figure  2.8:  PERSPECTIVE  view  of  shock  generator  and  cylindrical  model.  View  point  is 
from  upstream,  slightly  offset  to  one  side. 

surface  must  be  relieved  at  a  fine  angle  to  ensure  that  the  shock  on  the  exterior  surface  remains 
attached  to  the  leading  edge. 

Optimal  utilization  of  the  test  section  volume  is  achieved  when  none  of  the  nozzle  flow  spills 
underneath  the  shock  generator.  For  modest  flow  deflection  angles  a  circular  cylinder  with  diameter 
equal  to  the  nozzle  exit  diameter  is  a  good  approximation  to  the  ideal  elliptic  cross-section.  In 
practice  a  slightly  smaller  diameter  is  required  to  allow  for  the  boundary  layer  inside  the  nozzle  and 
the  finite  thickness  of  the  generator  (Figure  2.6).  This  arrangement  allows  the  size  of  the  model  to 
be  nearly  doubled  in  comparison  with  a  conventional  wedge  shock  generator.  Edge  effects  are  absent 
so  that  the  region  of  uniformly  compressed  flow  spans  the  entire  width  of  the  shock  generator. 

The  shock  generator  used  for  the  current  experiments  produced  a  flow  deflection  angle  of  6® 
and  exact  geometrical  details  are  given  in  figure  F.3.  Since  the  flow  in  the  nozzle  exhibited  some 
flow  angularity  and  Mach  number  nonuniformity  (Appendix  B)  the  actual  shock  angles  lay  slightly 
further  from  the  surface  of  the  wedge  than  indicated  by  an  inviscid  calculation  with  allowance  for  the 
boundary  layer  thickness.  The  measured  shock  angles  that  form  the  best  basis  for  any  computations 
are  listed  in  table  B.5. 
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Chapter  3  Experimental  results 
3.1  Introduction 

Holographic  interferograms  and  heat  transfer  data  were  obtained  for  three  different  test  conditions 
and  a  range  of  shock  impingement  locations.  The  data  are  contained  in  Appendix  C  along  with  an 
explanation  of  the  format  of  presentation.  This  chapter  begins  with  a  description  of  the  changes 
in  the  flow  structure  that  were  observed  as  the  location  of  the  impinging  shock  wave  was  varied. 
Following  this,  the  peak  heat  transfer  rates  for  the  type  IV  flow  field  are  examined  for  the  three 
test  conditions.  The  investigation  shows  that  increases  in  the  severity  of  the  phenomenon  at  high 
enthalpy,  owing  to  the  effects  of  dissociation,  are  less  significant  than  previously  expected.  Time 
resolved  heat  transfer  measurements  at  the  base  of  the  type  IV  jet  are  then  discussed  in  greater 
detail.  The  measurements  reveal  coherent  fluctuations  of  the  impinging  jet.  The  oscillations  are 
characterized  by  Strouhal  numbers  that  are  determined  from  the  spectra  of  the  heat  transfer  data. 
Supplementary  experiments  for  a  blunt  faced  cylindrical  model  are  used  to  reinforce  the  conclusions 
of  the  main  study.  Finally,  the  experimental  program  is  summarized  and  the  results  are  used  to 
motivate  the  development  of  improved  models  of  real  gas  effects  in  the  following  chapter. 


3.2  Description  of  the  flow  fields 

In  the  conduct  of  the  experiments,  and  in  the  description  of  them,  deliberate  emphasis  was  placed 
on  the  type  IV  jet  configuration.  This  necessarily  occurred  at  the  expense  of  other  interesting 
thermochemical  aspects  that  were  peripheral  to  the  main  objective.  As  a  minimum  requirement  it 
was  desirable  that  data  for  condition  C  should  span  the  entire  range  of  shock  impingement  locations. 
The  data  for  conditions  A  and  B  concentrated  almost  exclusively  on  the  type  IV  configuration. 

Twenty  shots  were  conducted  at  the  high  enthalpy  condition  C.  The  remaining  data  comprised 
nine  shots  at  the  low  enthalpy  control  condition  A,  and  four  shots  at  the  intermediate  condition  B 
to  aid  interpretation  of  the  results.  Vertical  translation  of  the  model  provided  a  range  of  shock 
impingement  locations  at  each  condition.  The  shock  generator  was  inclined  at  an  angle  of  6°  with 
respect  to  the  tunnel  centerline  for  all  conditions.  These  33  sequential  shots  (T5-824-856)  comprise 
the  core  holographic  and  heat  flux  data  that  are  contained  in  appendix  C. 

Several  series  of  preliminary  experiments  were  conducted  during  the  development  of  the  instru¬ 
mentation  for  this  project.  Initial  experiments  were  conducted  using  differential  interferometry  to 
investigate  the  flow  over  a  two-dimensional  cylindrical  model.  Subtle  differences  were  apparent 
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that  could  not  be  reconciled  with  the  original  data  of  Edney  [8].  Confirmatory  experiments  with  a 
spherical  model  at  low  enthalpy  successfully  reproduced  Edney’s  observations.  The  differences  were 
therefore  ascribed  to  three-dimensional  effects.  Brief  references  to  these  additional  data  are  made 
when  appropriate. 

3.2.1  Flows  without  shock  impingement 

The  flow  about  the  cylinder  without  shock  wave  impingement  is  shown  for  each  of  the  three  test 
conditions  in  figures  C.l,  C.14  and  C.20  respectively.  Appendix  B  contains  a  quantitative  assessment 
of  these  results.  In  all  three  cases  the  models  and  data  agree  tolerably  well.  Closed  fringes  visible 
in  the  interferograms  near  the  surface  of  the  body  generally  indicate  the  location  of  the  stagnation 
point.  In  these  cases  of  undisturbed  flow  they  naturally  lie  on  the  centerline. 

3.2.2  Type  I  flows 

Figure  C.21  demonstrates  a  type  I  flow.  The  bow  shock  shape  and  heat  transfer  distribution  remain 
undisturbed  since  the  interaction  occurs  downstream  of  the  sonic  line.  Preliminary  experiments 
using  differential  interferometry  indicated  no  significant  variations  with  stagnation  enthalpy  or  for 
a  spherical  model  at  low  enthalpy.  Real  gas  effects  are  known  to  influence  the  process  of  Mach 
reflection  in  a  relaxing  gas  (Hornung  et  al.  [21]).  Investigation  of  these  subtle  effects  on  the  type  I-II 
transition  process  requires  finer  variation  of  the  shock  impingement  location  than  was  appropriate 
for  the  goals  of  the  experiments  described  here. 

3.2.3  Type  II  flows 

As  the  impingement  point  was  raised  the  transition  to  type  II  flow  was  observed.  In  figure  C.22  a 
Mach  stem  connects  the  A-pattern  at  the  impingement  point  to  the  bow  shock  of  the  cylinder.  At 
the  intersection  point  a  strong  vortex  sheet  is  observed  separating  the  upper  region  of  subsonic  flow 
from  the  lower  region  of  supersonic  flow.  No  transmitted  wave  is  resolvable  at  the  upper  end  of  the 
Mach  stem  which  exhibits  a  point  of  inflection  where  it  joins  with  the  bow  shock  wave.  The  upper 
portions  of  the  fringe  pattern  and  heat  transfer  distribution  remain  symmetrical.  A  saddle  point  is 
observed  in  the  interferometric  fringe  pattern  above  the  impingement  point.  This  corresponds  to 
the  increased  density  of  the  gas  downstream  of  the  Mach  stem. 

Differential  interferometry  with  a  cylindrical  model  at  low  enthalpy  also  failed  to  resolve  a  trans¬ 
mitted  wave  at  the  upper  end  of  the  Mach  stem.  This  generic  feature  is  expected  for  type  II  flows 
(Figure  A. 4).  Figure  3.1  illustrates  the  single  instance  in  which  this  wave  was  observed  for  low 
enthalpy  flow  over  a  sphere.  This  reproduces  the  result  of  Edney  [8]  that  is  indicated  in  frame  2 
of  figure  A.l.  Curvature  of  the  bow  shock  wave  out  of  the  plane  of  symmetry  appears  to  fix  the 
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Figure  3.1:  Type  II  interaction  for  a  050mm  spherical  model  at  condition  A. 
location  of  the  Mach  stem  downstream  of  the  sonic  line. 

3.2.4  Type  III  flows 

Raising  the  impingement  point  caused  the  Mach  stem  to  merge  with  the  bow  shock.  In  figure  C.23 
the  bow  shock  above  the  impingement  point  exhibits  little  curvature  and  remains  nearly  normal 
below  the  centerline  of  the  cylinder.  A  closed  fringe  slightly  below  the  centerline  indicates  the 
location  of  the  stagnation  point.  The  saddle  point  in  the  fringe  pattern  above  the  impingement 
point  leads  to  a  secondary  maximum  in  density  away  from  the  surface  of  the  model.  A  sonic  line 
must  exist  in  the  flow  between  the  model  and  the  shear  layer.  The  density  decreases  downstream  of 
the  saddle  as  the  flow  accelerates  through  the  sonic  line.  The  shock  standoff  distance  has  increased 
only  slighty  from  its  value  without  shock  impingement.  Heat  flux  measurements  shown  in  figure 
C.23  indicate  a  slightly  asymmetric  fullness  below  the  centerline.  Although  the  Mach  stem  and  bow 
shock  have  merged,  this  flow  does  not  exhibit  the  obvious  asymmetries  of  the  shock  layer  that  are 
observed  at  higher  impingement  points  (Figure  C.24).  This  effect  is  possibly  due  to  the  behavior  of 
the  sonic  line  that  is  discussed  in  section  4.4.2. 

3.2.5  Displacement  of  the  stagnation  point 

In  the  current  experiments  Mach  numbers  of  3.0  and  0.3  were  typical  for  the  supersonic  and  subsonic 
streams  of  the  shear  layer.  On  dimensional  grounds  the  inertia  of  the  high  Mach  number  stream 
dominates  the  pressure  forces  of  the  subsonic  stream.  Consequently  the  shear  layer  does  not  deflect 
as  it  approaches  the  surface  of  the  cylinder  (Figure  C.25).  This  approach  further  reduces  the  size  of 
the  sonic  throat  that  must  exist  between  the  shear  layer  and  the  body.  The  fringe  pattern  that  is 
observed  in  figure  C.25  supports  this  conclusion.  The  fringe  shift  remains  nearly  constant  through 
the  subsonic  region  of  the  shock  layer  and  then  decreases  beyond  the  point  of  closest  approach.  The 
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mass  flux  that  exits  the  shock  layer  through  this  throat  must  decrease  as  the  shear  layer  approaches 
the  body.  The  remaining  mass  flux  that  enters  across  the  bow  shock  must  pass  above  the  cylinder 
and  so  the  stagnation  point  is  forced  downward.  This  is  evidenced  by  the  maximum  fringe  shift  in 
figure  C.25  that  occurs  30°  below  the  centerline  of  the  model.  Similarly,  the  A-shock  intersection 
point  stands  away  from  the  surface  of  the  body.  Heat  transfer  rates  downstream  of  the  sonic  throat 
remain  steady  and  increase  in  magnitude.  Temporal  fluctuations  of  the  heat  transfer  rate  would  be 
indicated  by  a  large  deviation  about  the  mean  measurement.  Equivalent  flow  features  are  observed 
in  figure  C.2  at  low  enthalpy.  In  this  case  spurious  heat  transfer  fluctuations  resulted  from  inversion 
of  the  surface  temperature  records  at  low  signal  levels. 

3.2.6  Entrainment  of  the  stagnation  streamline 

Raising  the  position  of  the  impinging  shock  wave  slightly  causes  the  stagnation  streamline  to  become 
entrained  with  the  shear  layer.  Evidence  for  this  statement  at  low  and  high  enthalpy  is  obtained 
from  the  holographic  interferograms  in  figures  C.3  and  C.30.  Figure  3.2  is  an  interpretation  of  the 
flow  field  in  the  vicinity  of  the  stagnation  point.  We  observe  closely  spaced  fringes  lying  parallel 
to  the  wall  above  the  point  where  the  shear  layer  meets  the  body.  This  region  of  high  density 
gradient  is  presumably  caused  by  a  portion  of  the  shear  layer  fluid  being  turned  upward  along  the 
wall  of  the  model.  The  stagnation  point  must  necessarily  lie  below  these  stream  lines.  There  exists 
only  one  local  maximum  in  density  along  the  surface  of  the  body  in  each  case.  Enlargements  of 
the  interferograms  allow  the  fringe  order  at  the  maxima  to  be  determined.  The  orders  are  30  =b  1 
and  11  ±  1  in  the  low  and  high  enthalpy  cases  respectively.  Both  values  exceed  the  corresponding 
undisturbed  stagnation  point  fringe  shifts  of  13  ±  1  and  9  ±  1  respectively.  This  can  only  occur  for 
streamlines  that  have  been  entrained  into  the  shear  layer  and  mixed  with  the  lower  entropy  fluid 
from  the  supersonic  stream.  Peak  heating  occurs  at  the  same  angular  position  as  the  maximum 
fringe  shift.  The  relative  heat  transfer  and  density  differences  at  low  and  high  enthalpy  cannot  be 
compared  directly  since  they  are  strong  functions  of  the  shock  impingement  location.  Noticeably 
larger  variance  in  the  heat  transfer  data  near  the  stagnation  point  is  indicative  of  mixing  that  occurs 
in  the  shear  layers.  Short  exposure  times  of  approximately  30  ns  freeze  the  motion  of  the  vortical 
structures  that  are  visible  in  the  interferograms.  These  structures  develop  over  lengths  that  are 
shorter  than  the  shock  standoff  distance  for  all  three  test  conditions. 

3.2.7  Turning  of  the  shear  layer  at  the  surface 

The  portion  of  the  shear  layer  that  is  turned  at  the  stagnation  point  and  passes  below  the  cylinder 
generates  a  compression  wave  that  coalesces  in  the  supersonic  flow  below  the  shear  layer.  At  some 
distance  from  the  wall  this  wave  overtakes  and  steepens  the  third  wave  that  originated  at  the 
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Figure  3.2:  Interpretive  sketch  of  the  entrainment  of  the  stagnation  streamline. 

impingement  point.  The  compression  wave  is  observable  in  figures  C.3  and  C.30  as  a  series  of  closely 
spaced  fringes.  The  closed  fringe  that  marks  the  local  maximum  of  density  along  the  wall  forms  a 
large  balloon  shape  in  each  interferogram.  The  balloon  shape  envelops  a  fringe  saddle  that  forms 
adjacent  to  the  wall  and  the  absolute  density  maximum  that  occurs  behind  the  compression  wave. 
This  saddle-maximum  combination  is  a  manifestation  of  the  density  gradient  across  the  shear  layer 
that  passes  below  the  cylinder.  The  density  maximum  occurs  only  in  the  supersonic  flow  beyond  the 
edge  of  the  shear  layer.  Continuing  this  progression  we  see  in  figure  C.31  that  the  saddle-maximum 
feature  moves  closer  to  the  wall  as  less  of  the  shear  layer  passes  below  the  cylinder. 

3. 2,8  Transition  to  type  IV 

Edney  [8]  hypothesized  that  transition  to  type  IV  would  occur  when  the  angle  between  the  shear 
layer  and  the  wall  exceeded  the  maximum  deflection  angle  for  the  compression  wave.  The  flows 
shown  in  figures  C.3  and  C.30  already  approach  the  maximum  deflection  condition.  Spreading  of 
the  shear  layer  near  the  impingement  point  makes  more  precise  measurements  difficult.  In  figures 
C.4  &  C.5  at  low  enthalpy  and  figures  C.26  &  C.27  at  high  enthalpy  we  see  that  the  saddle-maximum 
feature  collapses  into  the  surface  of  the  body  and  the  heat  transfer  rates  increase  rapidly.  On  the 
basis  of  continual  increases  in  the  peak  fringe  shift  (20  ±  1  in  figure  C.27)  we  conclude  that  the 
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Figure  3.3:  Interpretive  sketch  of  detachment  of  the  compression  wave. 

stagnation  streamline  moves  closer  to  the  high  speed  side  of  the  shear  layer.  This  interpretation  is 
complicated  by  the  large  entropy  gradient  across  shear  layer  whereby  the  maximum  fringe  shift  no 
longer  coincides  exactly  with  the  stagnation  point.  The  large  variance  of  the  heat  transfer  traces 
indicates  unsteadiness  of  the  flow  in  the  vicinity  of  the  stagnation  point.  The  compression  wave 
adopts  a  confused,  three-dimensional  appearance  in  the  interferograms  and  appears  to  detach  and 
curve  around  the  stagnation  point  (Figure  3.3).  The  wall  pressure  behind  this  wave  is  finally  sufficient 
to  overcome  the  inertia  of  the  supersonic  stream,  which  controls  the  location  of  the  stagnation  point, 
and  turn  the  shear  layer  substantially  above  the  cylinder.  At  several  locations  across  the  span  of  the 
image  in  figure  C.5  we  begin  to  see  the  appearance  of  an  upward  propagating  oblique  wave.  This 
shows  that  the  overtaking  intersection  of  the  compression  wave  and  the  bow  wave  has  failed  and  a 
second  A-shock  appears.  We  note  by  analogy  with  the  type  VI  transition  process  that  the  appearance 
of  this  feature  must  occur  before,  and  is  independent  of,  the  detachment  of  the  compression  wave 
from  the  wall.  The  inverted  A-shock  feature  becomes  more  prominent  in  figures  C.28  and  C.29  where 
the  maximum  fringe  shift  increases  to  24  ±  1 . 

3.2.9  Type  IV  flows 

Once  the  compression  wave  detaches  from  the  surface,  the  standolf  distances  must  be  controlled  by 
a  balance  of  the  mass  flux  entering  the  shock  layer.  These  influences  are  observed  in  figures  C.8  & 
C.9,  C.15  &  C.16  and  C.34  k  C.35  at  low,  medium,  and  high  enthalpy  respectively.  All  of  these 
flows  lie  in  the  regime  that  is  classified  as  type  IV.  As  the  impingement  point  moves  upward  the 
second  A-point  must  stand  away  from  the  wall  to  accommodate  the  increased  mass  flux  that  must 
pass  below  the  body.  Simultaneously,  the  mass  flux  passing  above  the  body  reduces  and  the  first 
A-point  moves  closer  to  the  body.  As  the  second  A-point  moves  away  from  the  body  the  type  IV 
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Test  condition 

A 

B 

c 

Undisturbed  fringe  shift 

13  ±  1 

10.5  ±  1 

9±1 

Interferogram 

Figure  C.8 

Figure  C.15 

Figure  C.34 

Type  IV  fringe  shift 

55  ±2 

29.5  ±  1 

27  ±  1 

Ratio 

4.2  i  0.5 

2.81  ±0.4 

3.0  ±0.4 

Table  3.1:  Stagnation  point  fringe  shift  for  type  IV  flows  with  cylindrical  model. 

jet  forms  with  a  third  oblique  wave  crossing  the  path  of  the  stagnation  streamline  (Figure  3.4). 
A  strong  terminating  shock,  that  originated  as  the  detached  compression  wave,  lies  close  to  the 
stagnation  point.  The  stagnation  point  fringe  shift  may  be  measured  and  compared  with  the  value 
without  shock  impingement  (Table  3.1).  These  measurements  indicate  a  trend  towards  reduced 
density  intensification  at  high  enthalpy.  Care  must  be  taken  in  comparison  of  these  highly  localized 
measurements  that  are  sensitive  to  variations  across  the  span  of  the  model.  Spanwise  nonuniformity 
is  most  pronounced  in  figures  C.32,  C.33,  C.36  k  C.37  where  multiple  images  of  shock  fronts  are 
observed  along  the  line  of  sight.  The  nonuniformity  is  ascribed  to  the  slight  out-of-plane  curvature 
of  the  incident  shock  wave. 

Comparison  of  the  type  IV  thermocouple  data  shown  in  figures  C.8,  C.15  and  C.34  indicates  a 
similar  trend  of  reduced  heat  transfer  intensification  at  the  higher  enthalpy  conditions.  In  all  cases 
the  maxima  of  density  and  heat  transfer  occur  at  the  same  angular  location. 

We  see  that  throughout  this  progression  of  impingement  locations  the  stagnation  streamline 
crosses  a  successively  increasing  number  of  weaker  waves.  In  the  undisturbed  configuration,  the 
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streamline  crosses  a  single  strong  wave.  As  the  stagnation  streamline  traverses  the  shear  layer 
it  mixes  with  lower  entropy  fluid  that  has  crossed  two  weaker  oblique  shock  waves.  Raising  the 
impingement  point  further  causes  the  stagnation  streamline  to  pass  through  two  oblique  waves  and 
the  strong  detached  compression  wave  before  it  meets  the  wall.  As  the  standoff  distances  equalize 
the  flow  passes  through  a  third  oblique  wave  which  results  in  an  additional  decrease  in  the  stagnation 
entropy.  The  holographic  shadowgraph  shown  in  figures  C.6  and  C.7  demonstrates  propagation  of 
this  third  wave  along  the  supersonic  jet  by  multiple  reflections  from  the  hounding  shear  layers.  These 
additional  reflections  alternate  isentropically  between  the  constant  pressures  of  the  subsonic  flows 
above  and  below  the  jet.  The  behavior  of  the  stagnation  point  that  is  observable  in  the  sequence  of 
holographic  interferograms  reinforces  the  discussion  of  the  streamline  patterns  in  section  2.3.3. 


3.2.10  Curvature  of  the  supersonic  jet 

The  strength  of  the  oblique  wave  that  propagates  along  the  supersonic  jet  is  determined  completely 
by  the  free-stream  conditions  and  the  incident  shock  wave  angle.  Consequently,  the  pressure  gradient 
across  the  jet  is  inversely  proportional  to  its  width.  This  is  consistent  with  the  increased  jet  curva¬ 
ture  observed  in  figures  C.IO,  C.ll,  C.17,  C.18,  C.38  and  C.39  at  low,  medium  and  high  enthalpy 
respectively.  The  corresponding  maximum  fringe  shift  ratios  are  13  i  1  >■  50  i  2,  10.5  i  1  ^  24  i  1 

and  9  ±  1  — ^  24  ±  1  respectively.  These  values,  and  the  corresponding  heat  transfer  rates,  have 
decreased  slightly  as  the  jet  narrows.  The  stagnation  point  appears  to  form  at  the  lower  edge  of 
the  jet  in  all  three  cases.  Again,  this  interpretation  is  sensitive  to  the  entropy  gradient  across  the 
lower  shear  layer  that  must  entrain  the  stagnation  streamline  as  the  jet  curls  upward.  Figures  C.IO, 
C.ll,  C.40  and  C.41  illustrate  the  manner  in  which  the  low  entropy  jet  fluid  energizes  the  boundary 
layer  on  the  body  above  and  below  the  impingement  point.  Interactions  of  the  jet  shock  waves  with 
the  shear  layers  produce  large  vortical  structures  that  are  observed  to  convect  downstream  in  the 
interferograms. 

3.2.11  Detachment  of  the  supersonic  jet 

Figures  C.13,  C.19  and  C.42  demonstrate  that  only  slight  changes  in  the  shock  impingement  lo¬ 
cation  produced  curvature  of  the  jet  sufficient  for  it  to  pass  entirely  above  the  body.  Figure  3.5 
illustrates  this  process  by  which  the  stagnation  streamline  emerges  from  the  supersonic  jet.  Note 
the  extremely  rapid  decrease  in  the  heat  transfer  rates  as  the  disturbance  detaches  from  the  surface. 
This  underscores  the  observation  that  the  behavior  of  the  stagnation  streamline  is  controlled  entirely 
by  the  inertia  of  the  supersonic  stream.  The  stagnation  point  fringe  shift  ratios  are  13  i  1  ^  35  i  2, 

10.5  ±  1  ^  24  ±  1  and  9  ±  1  20  ±  1  respectively  for  the  three  interferograms.  Since  these  values 

are  only  slightly  less  than  for  the  corresponding  type  IV  flows,  we  conclude  that  the  type  IV  fringe 
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Figure  3.5:  Emergence  of  the  stagnation  streamline  from  the  supersonic  jet. 

shifts  were  influenced  by  inevitable  variations  of  the  peak  density  location  across  the  span  of  the 
cylinder. 

3.2.12  Type  V  and  VI  flows 

Translating  the  shock  impingement  point  towards  the  upper  sonic  line  caused  the  jet  width  to  de¬ 
crease  beyond  the  resolution  of  the  interferometer  (Figures  C.43,  C.44  &  C.45).  A  point  of  inflection 
is  apparent  on  the  bow  shock  below  the  impingement  point  in  these  images.  The  characteristic  type 
V  shock  wave  was  not  observed  for  the  range  of  parameters  studied  here  and  this  is  consistent  with 
the  absence  of  the  type  II  wave  that  was  discussed  in  section  3.2.3.  Raising  the  shock  impingement 
point  beyond  the  upper  sonic  line  produced  the  conventional  type  VI  interaction  shown  in  figure 
C.46. 

3.3  Summary  of  the  heat  transfer  data 

Figures  3.6,  3.7  and  3.8  summarize  the  heat  transfer  data  at  the  low,  medium  and  high  enthalpy 
test  conditions  respectively.  Each  plot  represents  a  compilation  of  the  heat  transfer  data  for  the 
entire  sequence  of  shots  at  each  test  condition.  Two  sets  of  points  are  shown.  The  first  set,  -h,  is 
the  maximum  instantaneous  heat  flux  recorded  at  each  circumferential  location  on  the  model.  The 
maximum  value  was  determined  by  scanning  all  shots  over  an  interval  0.5  4.0  ms  after  initiation 

of  the  flow.  This  scanning  interval  excludes  flow  starting  transients  and  extends  over  the  period  of 
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Figure  3.6:  Envelope  of  maximum  heat  transfer  rates  around  the  model  at  condi¬ 
tion  A.  The  anomalous  point  at  ^  —  35^  arises  from  contributions  to  the  average  beyond  the  ends 
of  the  sampling  period. 

useful  test  time.  The  second  set  of  points,  A,  results  from  the  same  maximization  procedure  however 
the  data  were  first  averaged  over  an  interval  of  IOO/jls.  We  conclude  from  these  statistics  that  the 
peak  heat  transfer  intensification  decreased  by  a  factor  of  two  as  the  tunnel  stagnation  enthalpy  was 
increased  from  4MJ/kg  at  condition  A  to  19MJ/kg  at  condition  C. 

3.4  Comparison  with  existing  data 

Precise  comparison  of  the  data  with  existing  results  is  complicated  by  the  large  number  of  parameters 
involved.  The  data  of  Edney  [8]  that  are  shown  in  figure  3.9  may  be  compared  with  the  results  for 
the  low  enthalpy  condition  A  (Figure  3.6).  Despite  the  differences  in  the  free-stream  conditions 
and  model  geometries,  the  comparison  validates  the  current  measurements  and  peak  heating  was 
observed  at  identical  circumferential  locations.  Table  3.2  lists  the  peak  heating  rates  observed  by 
other  workers.  Differing  spatial  and  temporal  resolution  of  the  instrumentation  that  was  used  by 
the  various  authors  limits  the  utility  of  this  comparison.  The  use  of  time  averaged  values  from  the 
current  study  represented  the  best  compromise.  Interpolation  between  the  results  of  Edney  [8]  and 
Wieting  &  Holden  [48]  gives  good  agreement  with  the  low  enthalpy  data  shown  in  figure  3.6. 

The  Mach  number  of  the  flow  in  T5  decreased  as  the  stagnation  enthalpy  was  increased  (Ta¬ 
ble  B.l)  and  this  complicates  interpretation  of  the  current  results.  Coupling  of  these  two  parameters 
in  the  data  with  a  single  nozzle  area  ratio  was  unavoidable.  In  order  to  decouple  the  two  effects, 
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Figure  3.7:  Envelope  of  maximum  heat  transfer  rates  around  the  model  at  condition  B 
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Figure  3.8:  Envelope  of  maximum  heat  transfer  rates  around  the  model  at  condition  C 
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Figure  3.9:  Envelope  of  maximum  heat  transfer  rates  about  a  spherical  model  (Edney  [8]);  Moo  —  7, 
Si  =5°. 

consider  the  experiments  of  Edney  [8]  at  low  enthalpy  with  Moo  —  4.6  and  Si  =  5.®  that  produced  a 
heat  transfer  intensification  of  6.5  (Table  3.2).  The  current  high  enthalpy  experiments  at  condition 
C  with  Mco  =  5.3  and  Si  =  6.°  produced  a  suprisingly  similar  heat  transfer  intensification  of  6.0 
(Table  3.2).  Despite  the  different  model  geometries  for  the  two  studies,  we  see  that  the  prediction 
of  dramatic  real  gas  effects  is  not  supported  by  the  current  data. 

Table  3.2  also  shows  that  the  high  enthalpy  data  reported  by  Kortz  [27]  lie  below  the  heat  transfer 
rates  measured  in  the  current  experiments  at  the  same  enthalpy  and  also  below  those  of  Wieting  & 
Holden  [48]  at  the  same  Mach  number  and  wedge  angle. 

3.5  Temporal  characteristics  of  the  data 

3.5.1  Introduction 

The  instrumentation  and  signal  processing  techniques  that  are  described  in  appendix  E  provided 
unusually  high  bandwith  measurements  of  the  heat  transfer  intensification  caused  by  shock  im¬ 
pingement.  These  measurements  reveal  temporal  fluctuations  of  the  type  IV  jet  flow  that  have  not 
previously  been  reported.  Edney  [8]  noted  the  possibility  of  jet  oscillations  that  were  observed  in  the 
jet  impingement  study  of  Henderson  [17]  but  did  not  observe  them.  Recent  numerical  simulations 
of  the  type  IV  flow  field  also  indicate  the  existence  of  jet  fluctuations.  Several  authors  (Prahbu  et 
al.  [36],  Klopfer  k  Yee  [25])  have  reported  difficulty  in  obtaining  converged  steady  state  solutions 
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Reference 

Afoo 

Enthalpy 

Geometry 

q 

Edney  [8] 

4.6 

5.° 

Low 

Sphere 

6.5 

4.6 

10.° 

Low 

Sphere 

7.5 

4.6 

15.“ 

Low 

Sphere 

10. 

7.0 

5.“ 

Low 

Sphere 

8. 

Wieting  et  al.  [48] 

8.0 

10.“ 

Low 

Cylinder 

13. 

8.0 

12.5° 

Low 

Cylinder 

18. 

8.0 

15.° 

Low 

Cylinder 

25. 

Kortz  [27] 

9.9 

10.° 

20.  MJ/kg 

Cylinder 

V 

Figure  3.6 

9.9 

6.“ 

3.9  MJ/kg 

Cylinder 

13.5 

Figure  3.7 

6.3 

6.“ 

12.  MJ/kg 

Cylinder 

7.5^ 

Figure  3.8 

5.3 

6.“ 

19.  MJ/kg 

Cylinder 

6. 

Table  3.2:  Comparison  of  measured  heating  rates  from  previous  studies.  ^  These  measurements  were 
based  on  a  limited  number  of  data  points  and  may  not  resolve  the  maximum  heat  transfer  rate. 


for  type  IV  interactions.  Gaitonde  [15]  and  Zhong  [50]  addressed  the  issue  of  time  accuracy  of  the 
algorithms  that  were  used  and  independently  concluded  that  the  oscillations  are  attributable  to  the 
physics  of  the  problem.  Evidence  for  unsteadiness  is  also  observed  in  many  of  the  interferograms 
that  are  contained  in  appendix  C.  Unsteady  structures  appear  to  be  generated  at  the  jet  impinge¬ 
ment  point  and  convect  along  the  surface  of  the  body  (e.g.  Figures  C.8  &  C.IO).  Here  we  begin 
by  establishing  the  capability  of  the  current  experiment  to  resolve  the  temporal  fluctuations.  This 
is  achieved  by  studying  the  background  signal  levels  for  flows  without  shock  impingement.  Typical 
examples  of  type  IV  flows  are  then  selected  from  the  data  in  appendix  C  and  the  time  histories  of 
the  heat  transfer  rates  in  the  vicinity  of  the  jet  impingement  point  are  examined.  Frequency  spectra 
are  computed  and  the  behavior  is  characterized  by  a  dimensionless  Strouhal  number. 

3.5.2  Flows  without  shock  impingement 

Figures  3.10-3.12  demonstrate  the  constancy  of  the  test  conditions  during  the  period  of  time  over 
which  the  results  were  obtained.  Figure  3.11  reveals  some  undertailoring  at  condition  B.  This 
condition  was  used  for  only  4  shots  to  obtain  a  confirmatory  data  point  and  was  not  optimized.  The 
holographic  interferograms  were  exposed  at  ^  =  2.125  ms,  1.625  ms  and  1.275  ms  for  conditions  A, 
B  and  C  respectively.  The  heat  flux  traces  demonstrate  some  digitiser  quantisation  jitter  and  noise 
that  were  amplified  by  the  deconvolution  procedure.  Signal  levels  for  f  <  0  indicate  the  magnitude 
of  this  effect  that  was  exacerbated  by  the  requirements  of  wide  dynamic  range  and  high  bandwidth. 
Any  signals  above  this  background  level  are  real  and  indicate  flow  fluctuations  that  arose  during 
operation  of  T5. 
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Figure  3.10:  Temporal  characteristics  of  the  flow  without  shock  impingement  at 
CONDITION  A;  Shot  T5-855.  Top  row  plots  are  the  nozzle  reservoir  pressure.  Bottom  row  plots 
are  the  stagnation  point  surface  temperature  and  the  heat  flux  infered  from  it. 
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Figure  3.11:  Temporal  characteristics  of  the  flow  without  shock  impingement  at 
CONDITION  B;  Shot  T5-856.  Top  row  plots  are  the  nozzle  reservoir  pressure.  Bottom  row  plots 
are  the  stagnation  point  surface  temperature  and  the  heat  flux  infered  from  it. 
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Figure  3.12:  Temporal  characteristics  of  the  flow  without  shock  impingement  at 
CONDITION  C;  Shot  T5-824.  Top  row  plots  are  the  nozzle  reservoir  pressure.  Bottom  row  plots 
are  the  stagnation  point  surface  temperature  and  the  heat  flux  infered  from  it. 


3.5.3  Unsteadiness  of  the  type  IV  flow 

Figures  3.13“3.15  demonstrate  the  time  history  of  the  surface  heat  flux  in  the  vicinity  of  the  impinge¬ 
ment  point  for  three  representative  cases  at  low,  medium  and  high  enthalpy  respectively.  Corre¬ 
sponding  frequency  spectra  for  the  thermocouple  gauge  located  at  ^  =  —25°  and  the  nozzle  reservoir 
pressure  are  indicated  in  figures  3.16-3.18.  The  spectra  were  obtained  from  a  window  of  width  1.0  ms 
disposed  symmetrically  about  the  laser  firing  time.  Sharp  peaks  in  the  spectra  at  high  frequencies 
are  a  manifestation  of  digitizer  jitter  and  the  sensitivity  of  the  results  to  the  deconvolution  filter 
parameters  was  negligble. 

The  low  enthalpy  results  shown  in  figures  3.13  and  3.16  demonstrate  an  oscillotory  behaviour  that 
has  its  peak  power  density  at  a  frequency  of  6.5  kHz.  Noise  levels  in  the  deconvoluted  signals  appear 
at  least  an  order  of  magnitude  below  these  peaks.  Each  fluctuation  cycle  is  observable  as  a  sawtooth 
shaped  perturbation  of  the  surface  temperature  record.  This  shape  reflects  equation  E.7  which 
describes  the  response  of  the  gauge  substrate  to  an  impulsively  applied  thermal  load.  Additional 
realizations  of  this  flow  that  are  shown  in  figures  C.6  and  C.19  produced  somewhat  broader  peaks  that 
extended  from  5-15  kHz.  Computer  graphical  0-t  color  plots  of  the  heat  flux,  that  are  not  reproduced 
here,  demonstrated  that  the  fluctuations  of  the  individual  gauges  are  spatially  correlated.  Further 
evidence  of  the  coherent  nature  of  these  fluctuations  is  provided  by  the  repeatability  of  the  peak 
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heat  flux  for  each  cycle.  Although  the  power  spectral  density  of  the  heat  flux  signal  exceeds  that  of 
the  nozzle  reservoir  gauge  by  at  least  two  orders  of  magnitude  over  the  range  of  frequencies  at  which 
oscillations  were  observed  (Figure  3.16),  data  from  a  transducer  placed  immediately  upstream  of 
the  model  would  be  required  to  conclusively  demonstrate  the  independence  of  the  fluctuations  from 
free-stream  perturbations. 

Figures  3.14  and  3.15  demonstrate  that  the  peak  heating  was  fortunately  achieved  for  only  a  short 
period  of  time  at  the  higher  enthalpy  test  conditions.  This  occured  because  the  shock  generator  was 
long  with  respect  to  the  diameter  of  the  model.  Small  variations  in  the  shock  angle  during  the  test 
time  produced  a  large  effect  on  the  position  of  the  impinging  wave.  The  flow  structure  is  extremely 
sensitive  to  this  location.  A  sequence  of  interactions  was  produced  during  each  run  as  the  shock 
wave  moved  downward  with  respect  to  the  model.  For  this  reason  the  type  III  interaction  in  figure 
C.25  also  indicates  the  peak  heat  transfer  profile  for  a  type  IV  interaction  that  occured  earlier  during 
the  run.  Similarly,  figure  C.42  represents  a  type  IV  flow  at  later  times.  The  observed  invariance  of 
the  peak  heating  rate  at  different  test  times  provides  further  confidence  in  the  flow  quality. 

Frequency  spectra  of  the  surface  heat  flux  at  conditions  B  and  C  are  shown  in  figures  3.17  and 
3.18.  The  quality  of  the  spectra  was  severly  limited  by  the  small  number  of  cycles  that  were  observed. 

3.5.4  Comparison  with  existing  data. 

Using  the  free-stream  data  from  table  B.l,  the  frequency  of  oscillation  for  condition  A  may  be 
expressed  as  a  Strouhal  number,  Sh  =  fD/uoo  ^  0.16±0.08.  Best  estimates  of  the  Strouhal  number 
are  0.11  ±  100%  at  condition  B  and  0.06  ±  100%  at  condition  C.  Caution  should  be  exercised  in  the 
interpretation  of  this  trend  since  stagnation  enthalpy  was  not  the  sole  parameter.  Oscillations  were 
observed  in  the  recent  time  resolved  computations  of  Gaitonde  [15]  and  Zhong  [50]  who  simulated 
the  experimental  results  of  Wieting  &  Holden  [48].  Both  computations  independently  predicted  the 
result,  5/i  «  1.3.  It  should  be  noted  that  convergence  of  the  oscillation  cycle  was  not  achieved. 
Grid  refinement  produced  continually  increasing  modulation  of  the  fluctuation  cycle  and  this  was 
much  less  than  that  shown  in  figure  3.13.  Contour  plots  of  the  results  demonstrated  features  that 
are  remarkably  similar  to  those  observed  in  the  holographic  interferograms.  Gaitonde  [15]  cites 
additional  experiments  by  Holden  that  indicated  a  Strouhal  number  in  the  range,  Sh  =  0.13  ^  0.45. 
This  agrees  tolerably  well  with  the  current  data  at  condition  A. 


3.6  Supplementary  experiments 


Holographic  interferometry  yields  an  image  that  is  integrated  across  the  width  of  the  test  section. 
Whilst  the  use  of  a  cylindrical  model  simplified  quantitative  analysis  of  the  data,  the  inevitable  three- 
dimensional  variations  blurred  the  fine  details  of  the  images.  To  overcome  this  difficulty  additional 
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Figure  3.13:  Temporal  variation  of  the  surface  heat  flux  for  type  IV  flow  at  condi¬ 
tion  A;  Shot  T5-846.  Top-bottom:  0  =  —15^,  —20°, —25°, -30°. 


Figure  3.14:  Temporal  variation  of  the  surface  heat  flux  for  type  IV  flow  at  condi 
TION  B;  Shot  T5-853.  Top-bottom:  0  =  -20°,  —25°, -30°, -35°. 


Figure  3.15:  Temporal  variation  of  the  surface  heat  flux  for  type  IV  flow  at  condi¬ 
tion  C;  Shot  T5-839.  Top-bottom:  9  =  -10°,  -15°, -20°, -25°. 


Figure  3.16:  Frequency  spectra  of  nozzle  reservoir  pressure  and  surface  heat  flux  (gauge  #18)  for 
type  IV  flow  at  condition  A;  Shot  T5-846. 


Figure  3.17:  Frequency  spectra  of  nozzle  reservoir  pressure  and  surface  heat  flux  (gauge  #18)  for 
type  IV  flow  at  condition  B;  Shot  T5-853. 
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Figure  3.18:  Frequency  spectra  of  nozzle  reservoir  pressure  and  surface  heat  flux  (gauge  #18)  foi 
type  IV  flow  at  condition  C;  Shot  T5-839. 
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tests  were  conducted  with  a  <^50  mm  blunt  cylinder  model.  The  flat  face  of  the  cylinder  was  oriented 
directly  towards  the  free-stream  flow.  This  configuration  reduced  the  bow  shock  curvature  so  that 
the  interaction  region  was  enlarged  and  the  short  path  length  across  the  blunt  face  of  the  model 
provided  high  spatial  resolution  and  acceptable  sensitivity.  Figures  C.47-C.50  successfully  resolve 
the  oblique  wave  structure  of  the  jet,  the  terminating  normal  shock  and  the  region  of  elevated  density 
behind  it.  They  provide  further  confirmation  of  the  decrease  in  stagnation  density  intensification 
with  increasing  tunnel  stagnation  enthalpy. 

3.7  Summary  of  the  experimental  results 

The  experimental  data  indicated  a  clear  trend  towards  reduced  stagnation  density  and  heat  transfer 
intensification  as  the  tunnel  stagnation  enthalpy  was  increased.  Critical  assessment  of  the  flow 
conditions  confirmed  the  validity  of  these  measurements.  Comparison  with  existing  data  allowed  a 
crude  decoupling  of  the  effects  of  Mach  number  and  stagnation  enthalpy.  It  was  concluded  that  the 
observations  do  not  support  the  existing  predictions  (Figure  1.7)  of  greatly  increased  heat  transfer 
intensification  at  high  enthalpy.  Care  must  be  taken  in  interpreting  this  conclusion  since  the  data 
have  been  reduced  to  Stanton  number  form,  St  =  q/pooUooho^.  The  dimensional  heat  transfer  rate 
therefore  increases  considerably  with  stagnation  enthalpy. 

Holographic  interferometry  demonstrated  the  controlling  influence  of  the  inertia  of  the  supersonic 
stream  on  the  flow  topology.  Quantitative  interpretation  of  the  fringe  shift  at  the  base  of  the  type  IV 
jet  was  complicated  by  three-dimensional  averaging  along  the  line  of  sight.  Additional  interferograms 
successfully  resolved  the  internal  structure  of  the  type  IV  jet  for  the  flow  over  a  blunt  faced  cylinder. 

Time  resolved  heat  flux  measurements  demonstrated  the  existence  of  coherent  fluctuations  of  the 
type  IV  jet  and  allowed  the  estimation  of  Strouhal  numbers  for  the  phenomenon. 
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Chapter  4  Models  of  the  influence  of  real  gas  effects 
4.1  Overview 

The  experimental  investigation  indicates  that  further  increases  in  shock  impingement  heat  transfer 
rates  due  to  real  gas  effects  are  less  significant  than  indicated  by  the  existing  variable  7  model 
(Section  1.2).  The  model  that  is  described  here  follows  the  pattern  of  the  existing  analysis  but 
incorporates  refined  models  of  the  influence  of  real  gas  effects.  Development  of  the  model  proceeds 
in  four  stages; 

1.  Normal  shock  solutions  are  obtained  using  the  ideal  dissociating  gas  (IDG)  model  due  to 
Lighthill  [29]. 

2.  Extension  of  solutions  to  oblique  shock  waves. 

3.  Mapping  of  oblique  shock  waves  into  the  p-6  plane  to  obtain  solutions  for  the  A-shock  interac¬ 
tion  points  that  control  the  formation  of  the  type  IV  jet. 

4.  Modeling  of  the  impingement  heating  using  the  equilibrium  reacting  flow  similarity  solution 
of  Fay  &  Riddell  [11]. 

Sections  4.2  and  4.3  are  devoted  to  the  preliminary  steps  1  and  2.  Careful  nondimensionalization  of 
the  problem  leads  to  a  tractable  formulation  and  facilitates  identification  of  the  important  parameters 
in  the  subsequent  steps. 

Three  different  solutions  of  the  above  model  are  obtained  in  order  to  study  the  influence  of 
nonequilibrium  thermochemistry; 

1.  A  frozen  solution  that  applies  when  A  <C  1. 

2.  An  equilibrium  solution  that  applies  when  A  1. 

3.  A  nonequilibrium  approximation  that  applies  when  A  1. 

The  dissociation  rate  parameter,  A,  is  defined  as  the  initial  rate  of  the  dissociation  reaction  down¬ 
stream  of  a  normal  shock  normalized  with  respect  to  the  rate  of  flow  over  the  body  (Equation  4.25). 
The  frozen  and  equilibrium  solutions  are  obtained  by  using  appropriate  thermodynamic  models  for 
the  shock  jump  calculations.  The  nonequilibrium  solution  arises  from  the  observation  that  the  shock 
waves  produced  at  the  A-points  differ  greatly  in  strength  (Figure  1.8).  It  is  assumed  that  the  weak 
type  IV  jet  waves  remain  frozen  (Shocks  1,  3  &  5  in  figure  4.10)  whilst  the  bow  wave  and  the 
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terminating  jet  shock  relax  to  equilibrium  (Shocks  2,  4  &  6  in  figure  4.10).  The  validity  of  the  ap¬ 
proximation  is  assessed  by  computing  the  reaction  rates  (Freeman  [14])  immediately  downstream  of 
the  translational  discontinuities.  All  three  solutions  are  computed  for  a  range  of  parameters  and  the 
results  are  used  to  demonstrate  generic  influences  of  real  gas  effects  on  shock  interaction  problems. 

Specific  predictions  of  the  model  are  presented  for  the  test  conditions  used  in  the  experimental 
investigation.  The  most  appropriate  choice  between  the  three  solutions  discussed  above  is  determined 
by  the  dissociation  rate  parameter,  A,  for  each  free  stream  condition.  The  predictions  are  then 
compared  with  the  experimental  data. 


4.2  Normal  shock  for  an  ideal  dissociating  gas 

We  begin  by  deriving  the  jump  conditions  for  a  normal  shock  wave  in  a  relaxing  gas.  The  problem  is 
known  in  the  literature  (Vincenti  &  Kruger  [45])  however  the  compact  dimensionless  form  presented 
here  encompasses  a  nonequilibrium  upstream  state  and  is  particularly  convenient  for  the  solution  of 
multiple  wave  interaction  problems.  Consider  the  relaxation  to  equilibrium  of  the  system, 

TVo-j-M^iV  +  A-fM  (4.1) 

where  here  Nitrogen  is  the  prototypical  gas  and  M  represents  a  non-reacting  third  body.  The 
translational  and  rotational  shock  thickness  is  neglected  and  we  consider  only  the  relaxation  of  the 
internal  modes  over  a  length  scale  comparable  to  that  of  the  fluid  motion.  The  mass,  momentum 
and  energy  conservation  equations  that  apply  across  the  one-dimensional  translational  discontinuity 
and  throughout  the  downstream  relaxation  zone  are, 


P2U2  =  piUi  (4.2) 

P2  +  /»«2  =  Pi  +  P«1  (4-3) 

h2  +  2^2  —  ^1  4"  2^1  (^•'^) 

where  p  is  the  density  of  the  gas,  u  is  the  velocity  normal  to  the  shock,  p  is  the  pressure  and  h  is 
the  enthalpy  per  unit  mass.  The  subscripts,  1  and  2,  refer  here  to  the  upstream  and  downstream 
states  respectively.  In  dimensionless  form  equations  4.2  and  4.3  become, 


p  =  1  + 


pit  =  1 

1  p-1 


P  p 


(4.5) 


(4.6) 


where  the  variables  are  normalized  with  respect  to  their  upstream  values  and  generally  the  notation 
^  zz  (\)2/<^\  applies.  The  dimensionless  parameter,  P,  arises  and  it  behaves  as  P  ^  for  a 


perfect  gas; 


We  use  Lighthilhs  model  [29]  for  the  thermodynamics  of  the  system  4.1, 

+  (4.8) 

/,.=  T[(4  +  a)r  +  a^d]  (4.9) 

Im 

where  k  is  Boltzmann’s  constant,  m  is  the  mass  of  one  atom  of  the  gas,  T  is  the  temperature, 
is  a  temperature  characterizing  the  dissociation  energy  and  a  is  the  dissociated  mass  fraction 
determined  from  the  number  densities,  a  =  •  Using  the  caloric  equation  of  state  4.9  and 

the  dimensionless  continuity  equation  4.5  the  dimensionless  energy  conservation  equation  becomes, 

{4  +  a2)~  +  a2  +  ^^Ho  (4.10) 

where  the  dimensionless  parameters  are  defined  as; 

0=^,  (4-11) 


The  conserved  stagnation  enthalpy,  Hq,  is  normalized  with  respect  to  the  dissociation  energy  of  the 
gas  and  K  is  the  normalized  specific  kinetic  energy  of  the  upstream  flow. 

Importantly,  for  a  nonequilibrium  binary  mixture  of  thermally  perfect  components,  three  param¬ 
eters  are  sufficient  to  define  the  state  of  the  upstream  gas.  Here  we  specify  P,  Hq  and  ay.  The 
remaining  parameters,  0  and  K,  are  then  given  by  the  following  identities,  obtained  from  ho  written 
in  terms  of  the  conditions  upstream  of  the  discontinuity; 


Ho-K  -ay' 


Ho  —  ay 

1  +  2^^' 
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From  equation  4.8  the  dimensionless  form  of  the  thermal  equation  of  state  becomes; 


p(l  +  Ol) 

p{l  +  02) 


(4.16) 


Eliminating  f  and  ©  from  equation  4.10  with  the  pressure  given  by  4.6  leads  to  the  following 
quadratic  equation  for  the  density, 


4  +  02 


+  02 


(4.17) 


One  solution  describes  the  variation  of  density  throughout  the  relaxation  region  downstream  of  the 
translational  shock.  The  remaining  solution  describes  relaxation  of  the  possibly  nonequilibrium 
upstream  state  without  any  discontinuity.  Both  solutions  are  parameterized  in  terms  of  the  dissoci¬ 
ation  mass  fraction,  02  ,  whose  rate  of  reaction  remains  indeterminate.  Note  that  this  formulation 
encompasses  the  case  of  a  recombination  shock  wave  where  oi  >  0  and  02  <  cvi  although  this  does 
not  arise  in  the  current  problem.  The  solutions  are  illustrated  for  typical  values  of  the  parameters 
in  figures  4.1  and  4.2. 


4.2.1  Frozen  shock  solutions 

Consider  the  frozen  solutions  that  apply  directly  downstream  of  the  translational  discontinuity.  Here 
a2  =  oci  and  the  solutions  of  the  quadratic  expression  4.17  are, 


PJ  ^ 


_ Z±2j - 

(l+ai)+2F(4+ai)* 


(4.18) 


In  the  strong  shock  limit,  P  — >  0,  this  simplifies  to. 


These  solutions  are  consistent  with  the  limiting  cases  of  a  partially  vibrationally  excited  diatomic 
gas  and  a  monatomic  gas  respectively.  Note  also  the  value  Pmax  —  which  produces  pf  =  1  and 
corresponds  to  sonic  flow  upstream  of  the  shock  wave. 


4.2.2  Equilibrium  shock  solution 

Equation  4.17  may  be  combined  with  an  equilibrium  expression  to  obtain  the  asymptotic  state  far 
downstream  from  the  translational  shock.  This  determines  the  endpoints  for  the  trajectories  shown 
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a2 

Figure  4.1:  Strong  shock  solutions,  P  — +  0,  of  the  conservation  equations  downstream 
OF  A  NORMAL  SHOCK  WAVE  IN  AN  IDEAL  DISSOCIATING  GAS.  The  upstream  conditions  are  taken 
as  P  ==  0,  ai  =  0  and  solutions  are  plotted  for  differing  values  of  Pq-  The  two  frozen  solutions  (A) 
are  marked  on  the  vertical  axis.  The  solid  curves  represent  the  relaxation  of  the  density  downstream 
of  the  translational  discontinuity.  For  moderate  values  of  Hq  the  density  increases  monotonically  as 
the  flow  relaxes.  For  Hq  >  the  density  initially  decreases.  As  Po  oo  the  solution  asymptotes 
to  p  =  4  which  is  the  strong  shock  limit  for  a  monatomic  gas.  At  these  enthalpies  the  model  is 
inaccurate  since  it  does  not  model  the  electronic  excitation  that  would  occur.  Each  of  the  solutions 
exhibits  a  singularity  at  Po-a2  =  0.  Beyond  this  singularity  the  specific  kinetic  energy  downstream 
of  the  relaxation  zone  becomes  negative  as  may  be  observed  from  equation  4.15.  The  second  set 
of  solutions  shown  by  the  —  •  •  •  —  curves  represent  the  relaxation  of  the  upstream  state  without 
a  translational  discontinuity.  This  effect  becomes  larger  at  high  enthalpies  as  the  upstream  state 
Ofi  =  0  deviates  further  from  equilibrium. 
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Figure  4.2:  Relaxation  of  a  nonequilibrium  upstream  state  behind  a  normal  shock 
WAVE  IN  AN  IDEAL  DISSOCIATING  GAS.  The  upstream  state  is  taken  as  P  =  0.01,  =  0.25  and 

Ho  ~  0.2  1.0.  We  see  that  the  frozen  solutions  (A)  are  displaced  from  the  vertical  axis  as  a  result 

of  the  upstream  dissociation  level.  The  interpretation  of  the  two  sets  of  solutions  that  are  indicated 
by  the  solid  and  dot-dash  curves  is  qualitatively  similar  to  figure  4.1  for  the  cases  Ho  =  0.6,  0.8, 1.0. 
The  case  Ho  —  0.4  corresponds  to  a  recombination  shock  where  the  exothermic  relaxation  process 
causes  the  density  to  decrease  downstream  of  the  shock.  The  case  i/o  =  0.2  is  physically  unrealistic 
since  for  values  of  Hq  <  c^i  the  specific  kinetic  energy  of  the  flow  is  negative. 
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Figure  4.3:  EQUILIBRIUM  SHOCK  SOLUTION  IN  THE  p-a2  PLANE.  The  free-stream  conditions  are 

P  =  0.025,  Ho  =  0.55,  =  1  x  10^  ai  =  0.05. 

in  figures  4.1  and  4.2.  For  the  equilibrium  ideal  dissociating  gas  (Lighthill  [29]), 

_  Pd 

1  ^  <^2eq  Peq 

Writing  this  in  dimensionless  form,  using  the  identity  4.14  and  expressing  the  temperature  ratio  in 
terms  of  the  density  ratio  leads  to  a  second  transcendental  equation  relating  ao  ^'Hd  p; 

Ofo  PH  ^^+^^eq)Plq 

(4.19) 

1  “  peq 

An  additional  dimensionless  equilibrium  constant  is  introduced, 

Pd  =  PdlP\-  (4-20) 

Numerically  solving  equation  4.17  with  02  =  c^2eq  given  by  4.19  gives  the  equilibrium  shock  solution 
for  specified  upstream  conditions  P,  Hq  and  oi.  This  solution  is  represented  graphically  in  figure  4.3 
by  the  intersection  of  the  two  curves  in  the  p-a  plane  defined  by  equations  4.17  and  4.19.  When 
solving  these  equations  numerically  the  solution  is  bracketed  in  the  intervals  G  [p/ ,  co)  and 
^  [^i3  !]•  This  enables  the  development  of  a  robust  solver  with  known  convergence  properties 
in  the  troublesome  weak  shock  limit.  For  the  routine  construction  of  p-6  loci  this  is  a  significant 


concern. 
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4.2.3  Nonequilibrium  solution 

The  rate  at  which  the  dissociation  reaction  proceeds  for  the  IDG  model  is  given  by  Freeman  [14], 

da2  doi2  dt  CpT’^ 


dx  dt  dx 


X  p  O 

(1  -  a2)e  ^  -  —ai 
Pd 


(4.21) 


where  C  is  the  reaction  rate  constant  and  rj  represents  its  temperature  dependence.  Defining  the 
characteristic  reaction  length, 

and  applying  the  previous  dimensionless  notation  this  becomes; 


dx  dx 


(1  -  a2)e  r  - 


^  P  2 


Pd 


(4.23) 


Eliminating  the  temperature  dependence  using  equations  4.16,  4.6  and  4.14  we  obtain, 

0  _  p^(l  +  a2)  u  24) 

f  2/f(p(P+l)-l)' 

Equations  4.23  and  4.24  represent  a  third  expression  in  terms  of  and  p  that,  when  combined  with 
equation  4.17,  determines  the  relaxation  of  the  density  between  the  frozen  state  given  by  equation 
4.18  and  the  equilibrium  solution  given  by  equations  4.17  and  4.19  simultaneously.  Integrating  this 
ordinary  differential  equation  in  the  downstream  direction  produces  the  relaxation  profiles  that  are 
shown  in  figure  4.4. 


4.2,4  Reaction  rate  parameter 

The  reaction  rate  immediately  downstream  of  the  translational  discontinuity,  when  normalized  with 
respect  to  the  rate  of  flow  over  a  blunt  body,  represents  a  useful  indication  of  the  departure  of  the 
flow  field  from  equilibrium.  This  leads  to  definition  of  the  dissociation  rate  parameter, 


A 


da2 

D  dao 

D  dan 

dt 

j  ui  dx 

U2 —  =  -j:- 

j  ui  dx 

D 

Pf^d' 


where  D  is  the  diameter  of  the  body,  id  is  the  characteristic  reaction  length  and  the  subscript,  /, 
refers  to  the  frozen  state  immediately  downstream  of  the  translational  discontinuity.  The  initial 
dissociation  rate,  is  determined  by  the  exponential  term  of  equation  4.23  for  the  frozen  state 

downstream  of  the  shock,  where  a2  —  ai*  Equation  4.18  gives  the  frozen  solution  for  the  density 
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Figure  4.4:  Nonequilibrium  solution  for  the  relaxation  downstream  of  a  trans¬ 
lational  DISCONTINUITY  IN  AN  IDEAL  DISSOCIATING  GAS.  The  free-stream  conditions  are 
P  =  0.0254,  Ho  =  0.56,  pd  =  8.3  x  10®,  ai  =  0.051,  rj  =  -2.5.  These  conditions  are  equivalent  to 
the  high  enthalpy  test  condition  C  that  was  used  in  the  experimental  portion  of  this  investigation. 
For  these  conditions  the  characteristic  length,  la  =  0.544mm,  so  that  for  a  model  diameter  of  40 
mm  the  flow  is  close  to  equilibrium. 
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Od  =  113200  K 
pd  =  130000  kg 

m  =  14.0  X  10-76.023  x  10^^  kg 
C  —  2.1  X  10^^  kg-^ 

7]  =  -2.5 _ 


Table  4.1:  Constants  for  IDG  model. 


and  the  temperature  dependence  follows  from  4.24; 


This  parameter  is  linearly  dependent  on  the  ratio  of  the  body  size,  D,  to  the  characteristic  reaction 
length,  £d-  The  remaining  terms  are  most  strongly  dependent  on  the  stagnation  enthalpy,  Hq.  This 
may  be  demonstrated  in  the  hypersonic  limit,  P  0,  with  ai  =  0.  The  frozen  density  ratio  reduces 
to  pf  -r  7/(1  +  8P)  7  and  hence  Q/f  7-/12K  7^/12i7o-  In  this  limit  we  obtain. 


A  2^ 


-r] 

^~49/12Ho 


(4.26) 


This  is  a  good  approximation  for  a  wide  range  of  the  parameters  since  the  general  form,  given  by 
equations  4.25  and  4.24,  depends  only  weakly  on  P  and  a^.  The  rationale  for  binary  scaling  that 
is  used  in  hypervelocity  shock  tunnels  is  apparent  from  the  form  of  equation  4.26  since  this  gives 
similarity  in  terms  of  and  A. 


4.2*5  Physical  constants  for  the  IDG  model 

The  physical  constants  for  the  IDG  model  of  nitrogen  that  were  used  for  all  computations  in  this 
report  are  given  in  table  4.1. 
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4.3  Extension  to  oblique  shock  waves 

4.3.1  Normal  components  of  the  free  stream  parameters 

The  extension  to  oblique  waves  is  obtained  by  the  conventional  technique  of  applying  the  previous 
solutions  to  the  component  of  velocity  normal  to  the  shock  front.  If  jS  is  the  angle  of  the  oblique 
shock  wave,  then  from  the  definitions  of  P  and  K  we  obtain  the  normal  components, 


piy\  sin^  P 


Pi 

sin^  P  ’ 


mufsinV  r.  ^  2o 

te. 


The  identity  4.15  then  gives  the  normal  component  of  the  enthalpy, 


(4.27) 

(4.28) 


Given  Pi^  and  ai  the  variation  of  density  downstream  of  the  oblique  shock  is  then  deter¬ 

mined  by  equation  4.17.  Solving  simultaneously  with  equation  4.19  gives  the  equilibrium  solution 
and  integrating  the  ODE  4.23  determines  the  intermediate  states.  Once  the  density  is  known  the 
momentum  equation,  4.6,  gives  the  pressure  variation  along  the  streamline.  The  variation  of  the 
flow  deflection  angle,  6,  is  determined  from  continuity  and  the  conservation  of  velocity  tangential  to 
the  shock  to  be, 

tan(/?  —  S)  =  — 7^.  (4.30) 

P 

This  result  carries  over  directly  from  the  perfect  gas  case  and  here  it  applies  throughout  the  relaxation 
zone. 


4.3.2  Variation  of  parameters  across  an  oblique  shock  wave 

To  facilitate  the  computation  of  multiple  shock  jumps  for  the  interaction  problem  we  must  find 
a  convenient  means  of  determining  P  and  Ho  downstream  of  the  oblique  shock.  Since  stagnation 
enthalpy  is  conserved, 

Ho.^  =  Ho,,  (4.31) 

In  order  to  determine  P2  observe  that  U2^  =  s^nd  so  K2  —  ^^"27  +  Further,  since 

uij,  =  U2t^  and  pu^  =  I,  it  follows  that 


H2  =  Hi 


(4.32) 
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Identity  4.15  then  yields  P2  and  the  computation  for  the  downstream  shock  wave  proceeds  as  before 
using  P2,  Ho^  and  02  as  the  new  upstream  condition.  Downstream  of  the  shock  we  must  adjust  the 


equilibrium  constant  according  to, 


(4.33) 


Again,  these  relations  apply  throughout  the  relaxation  zone. 


4.3.3  Features  of  IDG  oblique  shock  waves 

The  parametric  dependence  of  some  arbitrary  quantity,  (j),  downstream  of  an  oblique  shock  wave 
may  then  be  summarized  as, 

(f)  =  pd,x/id,r],l3]] 

where  P  behaves  as  P  ~  Pq  is  the  dimensionless  stagnation  enthalpy,  ai  is  the  upstream 

dissociation  level,  pd  is  the  equilibrium  constant,  x/id  is  the  dimensionless  distance  normal  to  the 
shock  front,  rj  is  the  temperature  exponent  of  the  reaction  rate  and  (3  is  the  shock  wave  angle.  To  gain 
some  insight  into  the  behavior  of  the  solutions,  consider  the  limit  of  infinite  Mach  number  with  the 
upstream  kinetic  energy  equal  to  the  dissociation  energy;  P  —  0,  Po  —  =  0.  We  will  consider 

only  the  frozen  and  equilibrium  limits  and  so  the  dependence  on  x/ld  and  r]  disappears.  The  density 
ratio  across  the  shock  wave  for  such  a  condition  is  illustrated  in  figure  4.5.  The  dual  mechanisms 
of  compressibility  and  dissociation  cause  non-monotonic  behavior  of  the  shock  density  ratio  with 
increasing  shock  angle.  Even  for  very  weak  shock  waves  the  Mach  number  normal  to  the  shock  wave 
is  sufficient  to  cause  the  density  to  approach  the  strong  shock  limit  for  frozen  chemistry.  Only  at 
larger  shock  angles  do  the  frozen  and  equilibrium  solutions  deviate  as  the  exponential  dependence  of 
the  reaction  rate  and  the  equilibrium  dissociation  fraction  on  the  normal  component  of  the  enthalpy 
becomes  significant.  For  the  typical  shock  tunnel  test  condition  that  is  illustrated  in  figure  4.6  the 
distinction  between  the  two  mechanisms  is  less  obvious  however  the  conclusion  remains  valid.  For 
weak  waves  the  density  rises  faster  than  the  dissociation  level  and  the  dissociation  rate  rises  at  a 
still  slower  rate.  The  important  observation  is  that  for  waves  of  moderate  strength  that  produce  low 
dissociation  levels,  the  reaction  rate  remains  small. 


4.4  IDG  oblique  shock  loci  in  the  p-6  plane 

4.4.1  Features  of  oblique  shock  loci 

Shock  interaction  problems  are  conventionally  represented  graphically  in  the  p-6  plane.  The  influence 
of  real  gas  effects  on  this  construction  is  illustrated  for  a  single  oblique  shock  wave  in  figure  4.7. 
Real  gas  effects  cause  the  frozen  and  equilibrium  loci  to  be  separated  in  proportion  to  the  enthalpy. 
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Figure  4.5;  Frozen  and  equilibrium  density  ratio,  equilibrium  dissociation  level,  and  initial  dissocia¬ 
tion  rate  as  functions  of  shock  angle  for  an  oblique  wave  in  an  ideal  dissociating  gas.  The  free-stream 
conditions  are  P  =  1  x  10~®,  Ho  =  1-0,  ai  =  0,  pd  =  1  x  10®. 


Figure  4.6;  Frozen  and  equilibrium  density  ratio,  equilibrium  dissociation  level,  and  initial  dissocia¬ 
tion  rate  as  functions  of  shock  angle  for  an  oblique  wave  in  an  ideal  dissociating  gas.  The  free-stream 
conditions  are  P  =  0.025,  Ho  =  0.55,  oi  =  0.05,  /3rf  =  1  x  10~.  The  discontinuity  in  slope  at  /?  =  33“ 
arises  intentionally  as  a  consequence  of  the  frozen  free-stream  dissociation  level,  ai. 


58 


For  endothermic  reactions  both  the  pressure  and  flow  deflection  angle  increase  across  the  relaxation 
zone.  The  relative  insensitivity  of  the  pressure  to  real  gas  effects  is  evident.  Nonequilibrium  effects 
are  most  pronounced  in  the  flow  deflection  behind  strong  waves  where  the  normal  component  of  the 
enthalpy  is  large.  For  nearly  normal  waves  the  tangential  component  of  velocity  is  of  the  same  order 
as  the  normal  velocity  and  so  relaxation  produces  a  large  change  in  the  flow  deflection  angle.  This 
may  be  observed  in  the  streamlines  that  are  mapped  into  the  p-6  plane  in  figure  4.7.  Streamlines 
are  defined  by  equation  4.17  and  are  parameterized  in  terms  of  0^2;  where  0:2  6  [ai,<^2eg]  for  each 
j3.  The  pressure  and  flow  deflection  angle  are  given  by  equations  4.6  and  4.30.  In  the  weak  shock 
limit  the  origin  of  the  loci  in  figure  4.7  do  not  coincide  when  ai  is  nonzero  and  sufficiently  large. 
This  is  a  consequence  of  relaxation  of  the  nonequilibrium  upstream  state  that  is  demonstrated  by 
the  behavior  of  the  second  solution  branches  in  figures  4.1  and  4.2. 

4.4*2  Mapping  of  the  frozen  sonic  line 

The  definition  of  the  speed  of  sound  for  a  reacting  gas  is  ambiguous  since  the  manner  in  which 
the  composition  changes  across  the  disturbance  must  be  specified  (Vincenti  &  Kreuger  [45]).  The 
governing  partial  differential  equations  for  the  steady  flow  of  an  inviscid  reacting  gas  exhibit  a  change 
from  elliptic  to  hyperbolic  type  when  the  local  flow  velocity  exceeds  the  frozen  speed  of  sound.  Thus 
it  is  appropriate  to  consider  the  Mach  number  with  respect  to  the  frozen  speed  sound  in  both  the 
frozen  and  equilibrium  limits  behind  the  shock  front.  Using  the  current  notation  the  frozen  Mach 
number  is, 

where  the  variation  of  P  through  the  shock  layer  is  given  by  equations  4.31,  4.32  and  4.15.  The 
mapping  of  the  sonic  line  into  the  p-6  plane  is  indicated  in  figure  4.7.  A  streamline  that  intersects 
this  sonic  line  describes  an  initially  subsonic  flow  that  relaxes  to  a  supersonic  equilibrium  state  as 
the  temperature  falls  behind  the  shock  front.  This  fact  was  exploited  by  Hornung  &  Smith  [22]  to 
explain  observations  of  the  detachment  of  an  oblique  shock  wave  from  a  wedge.  It  is  also  indicative 
of  a  downstream  stretching  of  the  sonic  line  behind  the  bow  shock  on  a  blunt  body. 

4.4.3  Shock  wave  interaction  problems  in  the  p-S  plane 

The  use  of  the  p-6  plane  representation  to  solve  shock  interaction  problems  has  been  discussed 
in  section  A. 3.  Here  solutions  may  be  constructed  using  loci  for  either  the  frozen  or  equilibrium 
equations  of  state.  The  resulting  solutions  apply  directly  at  the  intersection  point  or  far  from 
the  shock  fronts  respectively.  Near  and  far  are  interpreted  relative  to  the  characteristic  reaction 
length,  id-  The  nonequilibrium  solution  4.23  that  applies  for  length  scales  intermediate  to  the  frozen 
and  equilibrium  limits  remains  valid  in  regions  where  the  flow  is  one-dimensional  normal  to  the  shock 
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Figure  4.7:  p-6  PLANE  representation  of  an  oblique  shock  wave  for  an  ideal  dissociating 

GAS.  Solid  curve  is  the  frozen  locus.  Dashed  curve  is  the  equilibrium  locus.  — - are  streamlines 

for  shock  angles  /?  =  40° ,50°, 60° ,68. 5°, 75°,  80° ,85°, 87. 5° ,89°  and  89.75°  respectively.  —  •  —  is 
the  mapping  of  the  frozen  sonic  line  connecting  the  frozen  and  equilibrium  limits  (+).  The  free 
stream  conditions  approximate  test  condition  C  that  was  used  in  the  experimental  investigation; 
P  =  0.025, Tfo  =  0.55,  pd  =  1  x  lO’’,  oi  =  0.05. 


fronts  (Figure  4.8).  A  typical  computation  of  the  loci  illustrating  the  asymmetrical  interaction  of 
two  oblique  shock  waves  is  illustrated  in  figure  4.9.  Although  chemical  nonequilibrium  causes  subtle 
differences  in  the  Mach  reflection  transition  mechanisms,  these  effects  are  small  since  the  waves 
remain  relatively  weak  up  to  the  transition  point. 


Solution  unknown  under  influence  of  vorticity 
generated  by  curvature  of  the  transmitted 
shock  waves. 


Figure  4.8:  Extent  of  validity  of  p-S  plane  representation  for  the  interaction  of  oblique  shock  waves 
in  an  ideal  dissociating  gas. 
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Figure  4.9:  p-S  plane  representation  of  the  interaction  of  oblique  shock  waves  in  an 
IDEAL  DISSOCIATING  GAS.  The  shock  angles  are  /?i  —  16°,  /?2  =  "51.5°  and  the  numbering  of  the 
states  is  given  in  figure  A. 3.  The  free-stream  conditions  are  P  =  0.025jHo  =  0.55,  pd  =  1  x  10^, 
ai  =  0.05.  Solid  curves  are  the  frozen  shock  loci.  —  ...  —  are  the  equilibrium  shock  loci. 
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Figure  4.10:  Schematic  of  flow  field  and  shock  wave  notation  for  the  type  IV  interaction  in  the 
vicinity  of  the  impingement  point. 

4.5  A-shock  configuration  for  an  ideal  dissociating  gas 

4.5.1  Introduction 

The  A-pattern  of  shock  waves  that  is  observed  to  form  where  the  shock  wave  impinges  on  the  leading 
edge  of  the  blunt  body  is  one  of  the  most  clearly  recognizable  features  of  the  type  IV  flow  field.  For 
a  given  free-stream  condition  and  incident  shock  angle,  the  remaining  quantities  are  determined 
entirely  by  local  satisfaction  of  the  conservation  equations  at  the  impingement  point.  This  is  the 
fundamental  mechanism  that  causes  the  formation  of  the  type  IV  jet  and  so  analysis  of  the  influence 
of  real  gas  effects  on  the  shock  angles  and  density  field  is  appropriate.  The  idealized  flow  field  and 
the  notation  for  the  shock  waves  are  shown  in  figure  4.10.  We  proceed  by  mapping  the  interaction 
problem  into  the  p-6  plane  in  order  to  obtain  the  frozen,  equilibrium  and  approximate  nonequilibrium 
solutions. 


4.5.2  Formulation  of  the  problem 

A  methodology  for  determination  of  the  generic  states  A  and  B  across  a  shock  wave  was  described 
in  section  4.3.  Here  it  is  described  functionally  by  the  relationship, 

The  shock  jump  conditions  ,  J,  are  obtained  from  equation  4.17  combined  with  equation  4,18  in 
the  frozen  limit  or  with  equation  4.19  in  the  equilibrium  limit.  The  equilibrium  solution  is  only 
valid  sufficiently  far  from  the  intersection  point  that  curvature  of  the  shock  waves  due  to  relaxation 
may  be  neglected  (Section  4.4.3).  Normal  components  of  the  parameters  are  used  as  determined  in 
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section  4.3.  The  solution  of  the  A-shock  problem  is  then  given  by  the  system  of  equations, 


ip2,p2,62,P2,Ho^,a2)  =  J{Poo,Ho^,Olco,Pd,l32),  (4-34) 

{pi,Pi,h,Pi,Hoi,o:i)  =  J{Poo,Ho^,ocoo,Pd,l3i),  (4.35) 

{P3,P3,S3>P3,Ho3,o‘3)  =  JiPi,Ho^,ai,Pd/pi,P3)-  (4.36) 

Matching  the  pressure  and  flow  deflection  angle  across  the  shear  layer  that  originates  at  the  inter¬ 
section  point  we  have, 

P2  ~  (4*37) 

^2  —  +  ^3-  (4.38) 


Collectively,  equations  4.34,  4.35,  4.36,  4.37  &  4.38  represent  20  equations  in  the  20  variables  pi,2,3, 
Pi, 2, 3,  <^1,2,3,  7^1, 2, 3,  7foi,2.35  «1,2,3,  y^2,3  and  for  convenience  we  choose  /?i,  Poo,  «co  and  pd 

as  parameters.  For  a  given  free-stream  condition  and  shock  angle  this  system  describes  a  path  of 
solutions  for  each  from  which  we  may  study  the  influence  of  real  gas  effects.  The  equations  are 
solved  using  Newton’s  method  for  the  cases  of  frozen  and  equilibrium  chemistry  and  the  solution 
is  continued  in  Solution  curves  for  the  density,  dissociation  levels,  shock  wave  angles,  flow 

deflection  angles,  and  reaction  rates  are  plotted  in  figure  4.11.  The  uniqueness  of  the  solution  that 
is  obtained  is  not  guaranteed  however  the  solution  branch  that  we  study  is  consistent  with  the 
experimental  observations. 

Consider  next  the  computation  of  the  initial  reaction  rates  downstream  of  the  translational 
discontinuities.  Observe  that  the  free  stream  conditions  for  the  computation  of  shock,  3,  given 
by  equation  4.36,  are  based  on  state  1.  The  normalization  of  the  reaction  rates  must  therefore  be 
adjusted  to  a  common  basis  for  comparison  and  this  is  chosen  to  be  the  characteristic  reaction  length 
for  the  free  stream,  For  the  i-th  shock  we  then  have, 

dcxi  doti  dxij^  dXoo  39) 

dx OQ  dxij^  dxjj^  dxoo 

where  the  subscript,  N,  refers  to  the  distance  normal  to  the  shock  front.  The  first  term  on  the 

right-hand  side,  is  given  by  equation  4.23  using  the  normal  components  of  the  parameters  for 

the  i-th  shock.  The  remaining  terms  that  adjust  the  normalization  are; 

dxj^  _  1  _ 

^i  —  \N 

dXoo  _  .  _  Mqo 

diZ~  PooO^dC' 
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^  =  sin(ft+  Si_,). 

;  =  1 

The  notation,  i-j,  refers  to  the  sequence  of  upstream  shocks  involved  in  the  normalization  (e.g. 
for  shock  i  =  6;  G  [5,3,1]  with  i_i  =  5,  z_2  =  3  etc.).  Substituting  these  relations  into 
equation  4.39  and  using  the  definition  of  the  dimensionless  specific  kinetic  energy,  /i ,  we  obtain  the 
consistently  normalized  reaction  rate; 


docj 

dxoo 


daj 

dXij^r 


i-j~l  i^j  —  1 

JJ  /5i_ .  sin(/?i  + 


j=o 


j=l 


(4.40) 


4,5.3  Discussion  of  the  solution 

Interpretation  of  the  solutions  that  are  illustrated  in  figure  4.11  follows  from  the  results  presented  in 
figure  4.6.  We  see  that  the  influence  of  real  gas  effects  increases  strongly  with  increasing  stagnation 
enthalpy,  Hq.  Dissociation  is  negligible  behind  the  weak  impinging  wave  and  dissociation  levels 
increase  proportionally  for  the  moderate  strength  wave,  3,  and  the  nearly  normal  wave,  2.  The 
corresponding  increases  in  the  density  are  the  most  important  real  gas  effects. 

A  further  similarity  with  the  conclusions  drawn  from  figure  4.6  is  the  behavior  of  the  normalized 
reaction  rates.  The  reaction  rate  behind  the  impinging  wave,  dai/dxco,  is  negligibly  small  in  both 
the  frozen  and  equilibrium  limits.  The  angles  of  the  stronger  shock  waves  ,  /?2  and  /?3,  vary  in  the 
frozen  and  equilibrium  limits  and  so  produce  different  initial  reaction  rates.  Behind  shock  3  the 
initial  reaction  rate  in  the  equilibrium  limit  is  an  order  of  magnitude  smaller  than  that  behind  the 
nearly  normal  wave,  2.  This  is  a  significant  effect  because  it  controls  the  relaxation  of  the  density 
downstream  of  the  interaction. 

The  strongest  effect  on  the  measurable  wave  angles  is  observed  in  the  deflection  angle,  62,3  = 
S2  =  Si  63,  of  the  shear  layer  generated  at  the  impingement  point.  We  see  from  the  mapping  of 
the  streamlines  into  the  p-S  plane  in  figure  4.7  that  the  strongest  effects  on  the  flow  deflection  angle 
are  apparent  for  nearly  normal  waves.  This  causes  the  strong  initial  effect  on  the  flow  deflection 
angle  as  Hq  increases.  The  effect  reduces  as  the  dissociation  level  behind  shock  3  rises  at  still  larger 
values  of  Hq. 


4.5.4  Representation  of  the  nonequilibrium  approximation  in  the  p^S  plane 

Since  the  reaction  rates  downstream  of  the  weaker  oblique  waves,  1  and  3,  at  a  A-point  are  signifi¬ 
cantly  less  than  that  downstream  of  the  strong  wave,  2,  we  introduce  the  nonequilibrium  approxima¬ 
tion,  da^/dxoo  ^  das/dxoo-  It  is  assumed  that  waves  1  and  3  remain  frozen  whilst  wave  2  is  allowed 
to  relax  to  equilibrium.  This  provides  a  simple  means  of  assessing  the  impact  of  thermochemical 
nonequilibrium  on  the  flow  field,  within  the  constraints  of  the  current  local  model  that  is  predicate 


Figure  4.11:  DENSITY  variations,  dissociation  levels,  shock  wave  angles,  streamline 

DEFLECTION,  AND  REACTION  RATES  AS  FUNCTIONS  OF  STAGNATION  ENTHALPY,  AT  A  A- 

SHOCK  POINT  IN  AN  IDEAL  DISSOCIATING  GAS.  The  free-stream  conditions  are  Poo  =  0.025,  pd  — 
1  X  10'^,  ttoo  =  0,  A  =  10°.  Solid  curves  are  the  frozen  solutions.  —  •  •  •  —  are  the  equilibrium 
solutions. 
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on  the  assumption  of  one- dimensional  flow  normal  to  the  shock  fronts. 

Although  the  strong  wave,  2,  is  curved  we  can  still  map  certain  regions  of  the  flow  into  the  p-6 
plane  with  validity.  Three  valid  solutions  may  be  obtained; 

1.  The  frozen  solution  that  applies  at  the  intersection  point. 

2.  The  intermediate  solution  that  applies  under  the  above  approximation  where  only  the  strong 
wave,  2,  is  allowed  to  relax  to  equilibrium. 

3.  The  equilibrium  limit  that  applies  far  from  the  intersection  point. 

These  limits  are  mapped  into  the  p-6  plane  in  figure  4.12  along  with  their  representation  in  the 
physical  plane  that  is  shown  in  figure  4.13.  The  approximation  dao/dx^o  ^  das/dxoo  amounts  to 
the  requirement  that  the  shock  fronts  remain  straight  over  length  scales  that  are  intermediate  to  the 
two  disparate  relaxation  lengths.  Therefore  the  one-dimensional  solution  for  the  relaxation  behind 
a  straight  shock  front  may  be  used  to  map  the  streamlines  in  this  region.  The  extent  of  validity 
of  this  solution  is  limited  by  the  subsequent  relaxation  of  the  supersonic  region,  3.  Sufficiently  far 
downstream  the  equilibrium  solution  is  obtained. 

Since  P2  ^  P3  it  is  also  possible  to  infer  the  shear  layer  departure  angle  from  the  frozen  limit  in 
the  p~S  plane.  Note  the  interpretation  of  P  —  pf  pu^  as  the  ratio  of  pressure  force  to  inertia  force  (c.f. 
P  ^  1/AP).  Since  p2  =  Ps  and  Pn  >  P3  it  follows  that  the  inertia  of  the  frozen  supersonic  flow,  3, 
dominates  the  pressure  of  the  relaxing  subsonic  flow  behind  the  strong  shock.  This  determines  the 
vertical  departure  of  the  stream  line  from  the  frozen  intersection  point  in  figure  4.12.  Extrapolating 
upward  towards  the  intermediate  limit  we  see  that  this  effect  produces  a  further  curvature  of  the 
strong  shock,  2,  towards  normal.  This  is  addition  to  that  caused  by  the  difference  in  the  relaxation 
lengths. 

Numerical  calculations  of  the  three  limits  are  shown  in  figure  4.14.  Points  of  inflection  are 
produced  on  both  transmitted  wave  fronts  whilst  the  shear  layer  curves  monotonically  at  this  level 
of  approximation.  These  arguments  lead  to  the  flow  field  that  is  sketched  in  figure  4.13.  These  slight 
changes  in  the  shock  angles  are  of  secondary  importance  in  comparison  to  the  density  and  reaction 
rate  variations. 

The  aim  of  the  nonequilibrium  approximation  is  to  provide  an  indication  of  the  influence  of 
finite  rate  effects  without  recourse  to  full  numerical  simulation.  An  additional  possibility,  that  is  not 
pursued  here,  follows  from  the  observation  that  the  partial  differential  equations  that  describe  the 
steady  flow  in  the  supersonic  jet  are  hyperbolic.  Numerical  simulations  of  the  nonequilibrium  jet 
flow  using  space  marching  techniques  are  possible  if  the  shear  layers  are  treated  as  constant  pressure 
boundaries.  This  approach  shares  the  limitations  of  the  current  local  model  whereby  the  scale  of 
the  interaction  region  and  the  location  of  the  terminating  jet  shock  remain  indeterminate. 
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Deflection  angle  (deg.) 


Figure  4.12:  Representation  of  the  IDG  A-shock  configuration  in  the  p-6  plane.  Solid 
curves  are  the  frozen  shock  loci.  —  •  •  •  —  curves  are  the  equilibrium  shock  loci.  Solution  points 
are;  A-frozen  limit,  o-nonequilibrium  approximation,  square-equilibrium  limit.  -  and  —  • 

—  are  the  streamlines  behind  shock  2  in  the  nonequilibrium  approximation  and  equilibrium  limit 
respectively.  Free-stream  conditions  are  P  —  0.025,  =  0.56,  Oqo  =  0,  =  1  x  10^,  /?i  —  16  . 


Figure  4.13:  Schematic  representation  of  the  effects  of  finite  rate  dissociation  on  the  A-shock  config¬ 
uration. 
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Figure  4.14:  Effect  of  finite  rate  dissociation  on  the  shock  wave  angles  and  flow 

DEFLECTION  DOWNSTREAM  OF  A  A-SHOCK  CONFIGURATION.  Solid  curves  are  the  frozen  solutions, 

- curves  are  the  nonequilibrium  approximations,  —  ...  —  curves  are  the  equilibrium  solutions. 

Free-stream  conditions  are  P  =  0.025,  a^o  =  0,  =  1  x  10^,  /?i  =  16°. 

4.6  Type  IV  inverted  A-shock  for  an  ideal  dissociating  gas 

The  preceding  analysis  for  the  A-shock  configuration  may  be  extended  to  account  for  the  second 
inverted  A-shock  that  is  observed  to  appear  in  the  type  IV  interaction.  The  flow  geometry  and 
additional  notation  are  described  in  figure  4.10,  The  system  of  equations  4.34-4.38  is  extended  by 


the  addition  of, 

hy  P4y  <^4.)  =  Pd/pl,M-  (4-41) 

{hyhy  hy  P^y  Hosy  =  J (^3,  ^^63^  Pd/ (PlP3),  /^s)-  (4.42) 

P4  =  P3P5,  (4-43) 

64  =  ^3  -h  ^5-  (4.44) 


Equations  4.41-4.44  represent  an  additional  14  equations  in  the  14  variables  p4  5,  ^4^5,  64^5,  ^4^5, 
P04  51  <^4,5  Rfid  /?4^5  that  are  decoupled  from  the  solution  of  the  preceding  system  4.34-4.38.  Again, 
this  system  describes  a  path  of  solutions  for  each  Hq^  where  the  remaining  parameters,  Poo, 
pd  ^  1^1  y  are  held  constant.  Since  the  two  systems  of  equations  decouple,  the  A-shock  solutions 
illustrated  in  figure  4.11  apply  at  the  impingement  point  for  the  type  IV  configuration.  Solutions 
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of  the  additional  system  of  equations  4.41-4.44  that  apply  at  the  downstream  inverted  A-shock  are 
illustrated  using  the  same  format  in  figure  4.15. 

The  behavior  is  somewhat  different  from  that  shown  in  figure  4.11  for  the  particular  values  of  the 
parameters  that  have  been  chosen.  Here  the  incident  wave,  3,  is  stronger  than  the  reflected  wave,  5. 
In  the  equilibrium  case  the  reaction  rate  behind  wave  5  is  very  small  since  the  wave  is  weak  and  the 
flow  behind  the  stronger  wave  upstream  of  it  has  already  relaxed.  In  the  frozen  case  the  reaction 
rate  is  much  higher  because  of  the  departure  from  equilibrium  caused  by  the  now  frozen  wave,  3, 
upstream  of  it.  In  either  case  the  reaction  rates  for  the  weaker  waves,  3  and  5,  remain  an  order 
of  magnitude  smaller  than  that  for  the  nearly  normal  wave,  4.  Figures  4.11  and  4.15  combine  to 
describe  a  particular  instance  of  a  type  IV  jet  in  which  the  intermediate  wave,  3,  is  stronger  than 
the  surrounding  waves,  1  and  5.  We  will  return  to  this  important  point  later. 


«  0.20 


Figure  4.15:  Density  variations,  dissociation  levels,  shock  wave  angles,  streamline 

DEFLECTION,  AND  REACTION  RATES  AS  FUNCTIONS  OF  STAGNATION  ENTHALPY,  AT  A  TYPE 

IV  INVERTED  A-shock  POINT  IN  AN  IDEAL  DISSOCIATING  GAS.  The  free-stream  conditions  are 
Poo  —  0.025,  pd  =  I  10^,  Ooo  =  0,  /?i  =  16*^.  Solid  curves  are  the  frozen  solutions.  —  •  •  •  —  are 
the  equilibrium  solutions. 
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4.7  Heat  transfer  model 

4.7.1  Introduction 

The  preceding  sections  describe  the  use  of  the  IDG  model  to  study  the  thermochemical  processes 
that  occur  in  the  inviscid  type  IV  jet  flow.  The  solutions  that  are  obtained  determine  the  state  of 
the  gas  that  impinges  on  the  surface  of  the  body.  The  following  section  describes  the  reacting  flow 
stagnation  point  model  that  is  used  to  estimate  the  heat  transfer  rates.  In  subsequent  sections  the 
two  models  are  combined  to  obtain  the  final  result. 


4.7.2  Fay  Riddell  stagnation  point  heat  transfer  solution 

The  stagnation  point  heat  transfer  solution  of  Fay  &  Riddell  [11]  represents  the  canonical  example 
of  the  influence  of  finite  rate  chemistry  on  blunt  body  heat  transfer.  By  applying  the  Illingworth 
transformation  to  the  boundary  layer  equations  for  a  compressible,  non-equilibrium  reacting  flow, 
Fay  k  Riddell  demonstrated  that  self  similar  solutions  are  admitted  at  the  stagnation  point  for 
arbitrary  chemical  reaction  rates.  Fay  k  Riddell  implicitly  formulated  the  problem  including  the 
source  terms  and  outer  solution  for  axisymmetric  bodies.  This  affects  the  coefficients  of  the  similarity 
equations  and  therefore  the  dimensionless  velocity  and  enthalpy  profiles  across  the  boundary  layer. 
To  rework  the  Fay  k  Riddell  solution  and  re-correlate  the  computations  for  planar  flows  of  nitrogen 
is  non  trivial  and  has  not  been  undertaken  here.  The  transport  properties  for  nitrogen  are  assumed 
to  be  sufficiently  similar  to  Fay  k  Riddell’s  binary  mixture  approximation  for  air  that  the  results 
are  adequate  to  describe  any  flow  phenomena  observed  here. 

The  stagnation  point  heat  transfer  rate  when  the  boundary  layer  is  in  chemical  equilibrium  is 
correlated  by  an  expression  of  the  form  (Fay  k  Riddell  [11]); 


^  —  Clipw  fJ'Vj)  {pepe) 


\0.4 


du 

dx  I® 


10(^6  “  hyi)) 


1  +  (Le 


0.52 


1) 


hd 


(4.45) 


Here  q  is  the  heat  transfer  per  unit  area  and  time,  fi  is  the  viscosity,  is  the  transverse  velocity 

gradient  at  the  outer  edge  of  the  boundary  layer,  Le  is  the  Lewis  number,  hd  is  the  chemical  enthalpy 
of  the  dissociated  gas  at  the  edge  of  the  boundary  layer  and  the  remaining  symbols  are  defined  as 
before.  The  subscripts,  iv  and  e,  denote  conditions  at  the  wall  and  at  the  outer  edge  of  the  boundary 
layer  respectively.  The  constant  of  proportionality  is  Ci  ==  0.94  for  axisymmetric  bodies.  For  planar 
bodies  the  value  Ci  =  0.70  is  obtained  by  extrapolation  of  results  for  non-reacting  flow  (White  [47]). 
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4.7.3  Outer  flow  solutions 

In  order  to  apply  this  result  we  must  assess  the  velocity  gradient  term,  du^jdx^  at  the  stagnation 
point.  Conventionally  this  is  done  by  application  of  the  Newtonian  approximation  for  thin  shock 
layers.  Since  the  density  of  the  outer  flow  at  the  stagnation  point  is  constant  we  have, 

Pe{x)  +  ^PeUe{xf  =  Pe(0). 

The  X  coordinate  is  measured  along  the  surface  of  the  body  with  its  origin  at  the  stagnation  point. 
Applying  L’Hopital’s  rule  in  the  limit  as  x  — >  0  gives  the  conventional  result, 

2  _  (Py  I 

_  dx~  lo 

Pe 

The  Newtonian  approximation  gives  the  pressure  distribution  as,  pe(^)  =  Poo^^  siil^  f3{x).  For  thin 
shock  layers  where  the  shock  inclination,  /?,  is  closely  approximated  by  the  local  surface  slope,  x/R^ 
the  stagnation  point  velocity  gradient  becomes; 


due  I 
_  dx 


(4.46) 


where  R  is  the  radius  of  curvature  of  the  planar  or  axisymmetric  body.  Writing  the  result  in  this 
form,  the  influences  of  real  gas  effects  are  evident  in  terms  of  the  free-stream  velocity  and  the  shock 
density  ratio. 

The  thin  shock  layer  approximation  used  to  obtain  equation  4.46  fails  in  the  limit  R  oo.  This 
case  is  important  since  it  is  required  for  analysis  of  the  blunt  cylinder  heat  transfer  data  and  also  as 
an  approximation  for  the  impingement  of  the  type  IV  jet  where  the  jet  width  is  much  smaller  than 
the  body  diameter.  On  dimensional  grounds  the  dependence  of  the  velocity  gradient  must  be; 

=  func  ( Moo,  — )  • 

^oo  \  Poo  ) 


In  the  perfect  gas  case,  that  is  discussed  in  some  detail  by  White  [47],  the  density  ratio  is  determined 
uniquely  by  the  Mach  number  and  in  the  hypersonic  limit  this  asymptotes  to  White  (Figure 

7-7)  cites  computational  and  experimental  results  that  give, 

dUe  I  D 

^^^  =  0.15, 

'W-OO 


in  the  limit  A/  — >  oc  for  flow  of  a  perfect  gas  (7  =  1.4)  over  a  blunt  faced  body.  On  the  basis  of  the 
density  dependence  exhibited  by  equation  4.46  we  extrapolate  to  the  case  of  a  dissociating  gas  at 


high  Mach  numbers; 
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(4.47) 


with, 


C2  — 


1,  l//i!>0, 

0.15  X  0.26,  l//?.  =  0. 


(4.48) 


The  discontinuous  limit  that  is  indicated  in  equation  4.48  is  unrealistic  and  the  question  of  appro¬ 
priate  scalings  for  the  velocity  gradient  in  this  regime  remains  open.  This  level  of  approximation  is 
considered  adequate  for  the  current  conceptual  model.  Only  ratios  of  the  parameters  (1  and  (^2  are 
contained  in  the  final  result  and  they  arise  in  a  form  that  is  not  critical  to  the  conclusions  of  the 


model. 


A.7A  Stagnation  point  recombination  rate  parameter 

The  key  conclusion  drawn  from  the  nonequilibrium  solutions  of  Fay  k  Riddell  [11]  is  that  if  the 
recombination  rate  is  sufficiently  high,  so  that  the  gas  remains  close  to  equilibrium  throughout  the 
boundary  layer,  or  if  the  wall  is  catalytic  and  the  Lewis  number  remains  close  to  unity,  so  that 
the  diffusive  and  conductive  mechanisms  of  energy  transport  in  the  layer  are  similar  in  magnitude, 
then  the  influence  of  recombination  on  the  boundary  layer  is  adequately  represented  by  the  chemical 
component  of  the  total  enthalpy  that  is  used  as  the  forcing  gradient  in  correlations  of  the  heat  flux. 
Nonequilibrium  effects  were  observed  as  a  decrease  in  the  heat  transfer  rate  to  the  wall  only  under 
conditions  when  the  recombination  rate  was  slow  and  the  wall  was  non-cat aly tic.  Wall  catalycity 
was  unimportant  at  high  recombination  rates. 

Following  Fay  k  Riddell  [11],  the  departure  of  the  stagnation  point  boundary  layer  from  equilib¬ 
rium  may  be  assessed  by  comparing  the  lifetime  of  an  atom  at  the  outer  edge  of  the  boundary  layer 
with  the  diffusion  time  scale  across  the  boundary  layer.  For  the  ideal  dissociating  gas  model  that  we 
use  here  the  time  rate  of  change  of  the  mass  of  atoms  per  unit  mixture  mass  due  to  recombination 
is  given  by  the  second  term  of  the  rate  equation  4.21, 


da 


recomh 


-Cp'^T^—  =  -q^fP 

Pd 


(4.49) 


The  dimensionless  notation  and  parameters  that  are  used  here  are  developed  in  section  4.2.  The 
lifetime  of  an  atom  in  the  boundary  layer  scales  with  If  this  is  evaluated  at  the  outer 

edge  of  the  boundary  layer  and  normalized  by  the  velocity  gradient  (equation  4.47)  at  the  stagnation 
point  we  define  the  recombination  rate  parameter, 
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Figure  4.16:  Contour  plot  of  the  recombination  rate  parameter,  as  a  function 

OF  dimensionless  stagnation  enthalpy,  Ho,  and  equilibrium  density  ratio The  con¬ 
tours  values  are  indicated  on  the  plot  and  the  remaining  parameters  are  taken  as  P  ==  0,  (Too  =  0, 
T]  -  -2.5.  The  dependence  on  P  for  large  Mach  numbers  is  weak.  At  low  enthalpies  where  the 
dissociation  levels  are  low  the  contours  asymptote  to  slope  rj.  At  higher  enthalpies  E  is  insensitive  to 
Hq  and  exhibits  an  inverse  power  dependence  on  pd>  The  three  test  conditions  used  for  the  current 
experiments  are  indicated. 


Only  the  dimensional  terms  of  equation  4.49  multiplying  the  dissociation  level,  have  been  in¬ 
cluded  since  this  leads  to  the  form  that  arises  in  the  similarity  formulation.  The  temperature  ratio  at 
the  edge  of  the  boundary  layer  is  approximated  by  the  solution  for  the  equilibrium  shock  layer  4.24, 

0  /32(l  +  a) 

f  ~  2K[p{P^l)-iy 

The  boundary  layer  tends  to  equilibrium  for  E  >>  1  and  is  frozen  for  E  1-  The  nonequilibrium 
solutions  of  Fay  &  Riddell  [11]  indicate  that  nonequilibrium  effects  become  significant  only  for 
values  of  E  <  0.1  when  the  wall  is  non  catalytic  to  the  recombination  reaction.  Contours  of  the 
recombination  rate  parameter  in  the  universal  form,  E^^,  are  indicated  in  figure  4.16  along  with 
the  values  for  the  three  test  conditions  that  were  used  in  the  experimental  investigation.  Table  4.2 
indicates  the  values  of  the  recombination  rate  parameter  for  the  current  experiments  and  in  all  cases 
the  boundary  layer  is  justifiably  treated  as  equilibrium  since  E  ^  0.1. 
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Test  condition 

A 

B 

C 

C2^d 

0.080 

0.018 

0.025 

C2 

1 

1 

1 

id 

0.19  mm 

0.46  mm 

0.54  mm 

D 

40.6  mm 

40.6  mm 

40.6  mm 

s 

17. 

1.6 

1.9 

Table  4.2:  Predicted  recombination  rate  parameters  for  test  conditions  A,  B  &  C. 
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Further,  the  pressure  across  the  boundary  layer  is  constant,  —  l/F^oo?  Q-nd  the  thermal 

equation  of  state  gives  the  wall  density  ratio, 


pw  — 


1  ^oo 
PoqTw 


Since  the  wall  temperature  is  low  and  the  boundary  layer  is  in  equilibrium  we  have  taken  a^j  =  0. 
Substituting  these  dependencies  into  equation  4.51  we  obtain  the  Stanton  number  as  a  function  of 
the  relevant  dimensionless  parameters, 


Sihody  ~ 


^  C1C2  (  1  +  CTc 

RC(y 


0.38 


(1  4-  Cl'e) 


-0.2Sm-0.03 -~0.13 


*'  w 


Pe 


1  +  0.2 


(4.52) 


The  complete  dimensionless  dependence  of  the  stagnation  point  heating  problem  as  it  is  modeled 
here  becomes  SUody{Pco,  pd^  D / R^oo.Tw)^  The  remaining  variables,  ae  and  pe,  are 

determined  from  the  equilibrium  normal  shock  solution  in  terms  of  the  same  set  of  parameters. 
Variation  of  the  flow  variables  as  the  subsonic  flow  decelerates  to  the  stagnation  point  is  neglected 
since  the  Mach  number  is  low. 
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Figure  4.17:  Adaption  of  the  Fay  &  Riddell  stagnation  point  solution  to  the  jet  impingement  problem. 


4.8  Type  IV  jet  heat  transfer  solution 


4.8.1  Introduction 

In  this  section  we  combine  the  stagnation  point  heat  transfer  model  with  the  free-stream  conditions 
provided  by  the  inviscid  type  IV  jet  solution.  The  results  of  the  model  will  be  expressed  in  terms 
of  the  ratio  of  the  heat  transfer  rate  at  the  jet  impingement  point  with  respect  to  the  undisturbed 
stagnation  point  value.  Whilst  equation  4.52  applies  at  the  undisturbed  stagnation  point,  its  validity 
at  the  jet  impingement  point  requires  further  consideration. 

If  the  local  vorticity  and  out-of-plane  strain  rate  at  the  type  IV  jet  stagnation  point  remain 
small,  then  the  deformed  stagnation  streamline  meets  the  wall  symmetrically  at  90°.  Under  these 
restrictions  we  may  adapt  the  Fay  &  Riddell  solution  4.52  to  study  the  jet  impingement  problem 
as  shown  in  figure  4.17.  The  conditions  of  zero  vorticity  and  out-of-plane  strain  rate  aie  sensibly 
satisfied  for  streamlines  that  pass  through  the  jet  in  a  two-dimensional  shock  impingement  problem 
(Figure  4.10).  Vorticity  produced  by  curvature  of  the  oblique  jet  shocks  or  the  terminating  normal 
wave  is  neglected. 

Differences  in  the  flow  topology  for  three-dimensional  shock  impingement  configurations  were  dis¬ 
cussed  in  section  2.3.3.  Here  we  note  an  additional  difference  whereby  out-of-plane  vortex  stretching 
and  diffusion  alter  the  form  of  the  stagnation  point  at  the  base  of  the  impinging  jet.  The  current 
results  cannot  therefore  be  applied  to  non-planar  configurations  such  as  the  shock-on-fin  problem. 
Despite  this,  it  is  reasonable  to  assume  that  the  generic  influences  of  real  gas  effects  are  similar  for 
both  two  and  three-dimensional  flows. 

The  following  model  also  assumes  quasi-steady,  slowly  varying  conditions  at  the  stagnation  point. 
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4.8.2  Heat  transfer  intensification  due  to  jet  impingement 

We  now  proceed  to  apply  the  Fay  k  Riddell  [11]  solution  at  the  jet  impingement  point.  The 
parameters  describing  the  free-stream  conditions  are  given  by  the  inviscid  shock  interaction  solution 
(Sections  4.5  k  4.6).  The  conditions  upstream  of  the  normal  shock  at  the  base  of  the  jet  are  taken  to 
be  those  existing  at  state,  5,  behind  the  third  wave  produced  at  the  inverted  A-shock.  Any  additional 
waves  in  the  supersonic  jet  are  neglected  as  discussed  in  section  1.2.  We  then  have, 


=  23/4, 


Re^ 


(1  +  as) 


H-0.03 --0.13 


1+0.2- 


(4.53) 


The  Reynolds  number,  Res,  is  based  on  the  width  of  the  supersonic  jet,  Aa;.  As  before,  subsonic 
variations  of  the  equilibrium  state,  6,  behind  the  terminating  jet  shock  are  neglected  as  the  flow 
decelerates  to  the  stagnation  point.  The  intensification  of  the  heat  transfer  as  a  result  of  jet  im¬ 
pingement  is,  _ 

^  E,  (4.54) 

Qbody  body  Poo  ^oo  ^^body  V  -^^oo 

All  of  the  terms  in  equations  4.52  and  4.53  that  are  necessary  to  evaluate  this  expression  are  available 
from  the  inviscid  jet  solution  with  the  exception  of  Re^.  The  ratio,  arises  and  this  is 

determined  by  the  thermal  equation  of  state,  4.16,  to  be. 


_  hhh  f  1  +  0^5  ^ 

Twoo  Ts  P^PsPl  \  1  +  ttcxD  / 

The  ratio  of  the  Reynolds  numbers  is  similarly  determined  since  we  have  assumed  a  power  law 
dependence  for  the  viscosity; 


Re^  _  p^u^PoqAx 
R^oo  poo^ooP^R^ 


hhh 


I  /is  Ax 

y  Koo  D 


0.7 


(4.56) 


Substituting  equations  4.53,  4.52,  4.55  and  4.56  into  equation  4.54  we  obtain  the  following  expression 
for  the  heat  transfer  intensification; 


ij  et 
Qbody 


/  j.^  \  0.63 

(tS) 

/l  +  aeV”'"^  rRo  +  0.2a6' 

\l  +  ae/  \pe )  _i7o  +  0.2ae. 


(4.57) 


Observe  that  although  the  heat  transfer  rate  varies  inversely  with  the  square  root  of  the  Reynolds 
number  in  equation  4.53,  the  intensification  of  the  heat  transfer  due  to  jet  impingement  in  equation 
4.57  is  independent  of  the  Reynolds  number.  This  is  consistent  with  the  experimental  observations 
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of  Wieting  &  Holden  [48]  at  low  enthalpies.  No  variation  in  the  heat  transfer  amplification  was 
observed  as  the  Reynolds  number  was  changed  with  the  remaining  parameters  held  nearly  constant. 

The  last  three  factors  of  equation  4.57  are  close  to  unity.  This  is  a  consequence  of  the  conservation 
of  stagnation  enthalpy  at  all  points  in  the  flow  and  the  weak  additional  dependence  of  the  stagnation 
point  dissociation  level  on  density.  To  a  very  good  approximation  ae.  The  ratio  of  the  normal 


shock  density  ratios,  4^,  is  somewhat  less  that  unity  however  the  small  exponent,  —0.13,  reduces  its 

C2 

influence.  The  square  root  dependence  on  the  geometrical  factors,  and  >2  ,  is  explicit. 

The  remaining  two  factors  in  the  center  demonstrate  intensification  caused  by  increases  in  the  jet 
specific  kinetic  energy  and  density  respectively. 


4,8,3  Numerical  solutions 

The  factors  of  equation  4.57  representing  jet  density  and  kinetic  energy  intensification  control  most 
of  the  interesting  effects.  Their  behavior  is  determined  by  the  two  A-shock  solutions  of  sections  4.5 
and  4.6.  Equations  4.34-4.38,  4.41-4.44  and  4.57  determine  the  solution  which  has  the  form, 


4j  et  ^Ibody^^^ody 

Qbody]l  D 


func{Poo ,  ,Pd,otoo,0i)‘ 


(4.58) 


The  equations  are  solved  using  Newton’s  method  and  the  solution  is  continued  in  one  parameter 
(T/Ooo  A)  typical  values  of  the  remaining  parameters.  The  term,  ^  >2  " ^ — 
the  result  for  a  particular  geometry  and  remains  indeterminate. 


^^^body^^body  ^  gp^^^j^ligeS 


The  influence  of  dissociation  chemistry  appears  principally  through  the  parameter  i^Ooo-  Finite 
rate  effects  are  considered  in  the  approximate  intermediate  limit,  <C 

This  approximation  leads  to  three  solutions; 


1.  the  frozen  solution, 

2.  the  nonequilibrium  solution  under  the  above  approximation  where  only  shocks  2,  4  and  6  are 
allowed  to  come  to  equilibrium, 


3.  the  full  equilibrium  solution. 

The  validity  of  the  nonequilibrium  approximation  is  assessed  from  the  requirement  that  the  reaction 
rate  downstream  of  the  jet  shocks  (1,  3  fc  5)  be  small  relative  to  that  behind  the  undisturbed  bow 
shock  (i.e.  <  !)•  The  influences  of  aco  and  pd  are  numerically  significant  but  secondary 

in  terms  of  demonstrating  the  essential  behavior.  The  claim  therefore  is  that  the  solution  curves 
plotted  using  the  universal  form,  4.58,  in  figures  4.18-4.22,  illustrate  the  main  mechanisms  involved. 
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Figure  4.18:  Influence  of  dimensionless  stagnation  enthalpy,  Mach  number,  ^ 

l/y/P^  and  reaction  RATE  ON  TYPE  IV  HEAT  TRANSFER  INTENSIFICATION.  From  top:  Heat 

I  C2 

transfer  intensification,  ;  stagnation  density  intensification,  ^  and  jet  reaction 

Qbody  y 

rate  ratio,  4^/4^-  . curves  are  the  frozen  solutions.  —  •  ♦  •  —  curves  are  the  equilibrium 

solutions. - curves  are  the  nonequilibrium  solutions  in  the  limit 

The  parameters  are  Pqo  =  0.002,  0.010,  0.025,  =  1  x  10^,  Ooo  =  0  and  /?,  =  16°. 
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Figure  4.19:  Influence  of  impinging  shock  wave  angle,  /?i,  and  reaction  rate  on  type 
IV  heat  transfer  intensification  at  very  high  Mach  number;  Poo  =  0.002.  From  top:  Heat 


transfer  intensification, 


(}hody 


the  frozen  solution. 


Aa; 


-To — : -  and  iet  reaction  rate  ratio,  . curve  is 

axoo 

—  curve  is  the  equilibrium  solution. - curve  is  the  nonequilibrium 

solution  in  the  limit  ^he  parameters  are  i?o„  =  0.6,  pd  =  1  x  10^ 

and  Ooo  =  0. 


Figure  4.20:  Influence  of  impinging  shock  wave  angle,  /?i,  and  reaction  rate  on  type 
IV  heat  transfer  intensification  at  high  Mach  number;  Poo  =  0.010.  The  remaining 
parameters  and  format  are  the  same  as  figure  4.19. 
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Figure  4.21:  Influence  of  impinging  shock  wave  angle,  /?i,  and  reaction  rate  on  type 
IV  heat  transfer  intensification  at  moderate  Mach  number;  Poo  =  0.025.  The  remaining 
parameters  and  format  are  the  same  as  figure  4.19. 
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Figure  4.22:  Influence  of  impinging  shock  wave  angle,  /?i,  and  reaction  rate  on  type 
IV  heat  transfer  intensification  at  low  supersonic  Mach  number;  Poo  =  0,050.  The 
remaining  parameters  and  format  are  the  same  as  figure  4.19, 
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4*8.4  Discussion 

The  influence  of  real  gas  effects  on  the  shock  impingement  problem  may  be  observed  in  flgure  4.18. 
The  frozen,  nonequilibrium  and  equilibrium  solutions  are  plotted  as  a  function  of  dimensionless 
stagnation  enthalpy,  Naturally,  the  frozen  solutions  do  not  vary  with  stagnation  enthalpy. 

The  difference  between  the  equilibrium  and  frozen  solutions  increases  with  stagnation  enthalpy  as 
dissociation  produces  increasing  density  ratios  across  the  jet  shocks,  3  and  5.  This  effect  is  absent 
in  the  nonequilibrium  solution  that  in  fact  decreases  slightly  with  respect  to  the  frozen  solution. 
Since  cve  and  as  =  0  in  the  nonequilibrium  approximation,  the  effects  of  dissociation  on  the 

normal  wave,  6,  and  the  undisturbed  bow  shock  are  similar.  This  limits  the  deviation  of  the  frozen 
and  nonequilibrium  heat  transfer  ratios.  The  signiflcance  of  real  gas  effects  is  seen  to  increase  with 
decreasing  Mach  number  (or  increasing  Poo)*  The  lower  plot  in  figure  4.18  demonstrates  that  the 
validity  of  the  nonequilibrium  approximation  improves  with  increasing  Mach  number  and  decreasing 
stagnation  enthalpy. 

The  existing  perfect  gas  model  (Edney  [8],  Keyes  &  Hains  [24])  considers  the  effects  of  three 
parameters; 

1.  the  free-stream  Mach  number, 

2.  the  impinging  shock  angle  ,/?i, 

3.  the  ratio  of  specific  heats,  j. 

This  model  predicts  that  the  heat  transfer  intensification  increases  with  the  free-stream  Mach  number 
and  peaks  at  a  modest  value  of  the  impinging  shock  strength.  These  trends  are  reproduced  by  the 
frozen  solutions  in  figures  4.19-4.22  that  correspond  to  free-stream  Mach  numbers  of  Moo  =  19., 
8.7,  5.5  and  3.9  respectively  (7  =  4/3).  As  the  Mach  number  increases,  and  as  /?i  approaches  the 
maxima  in  figures  4.19-4.22,  the  influence  of  real  gas  effects  declines.  The  relative  influence  of  real 
gas  effects  is  quite  large  for  shock  angles  either  side  of  the  maxima. 

The  reaction  rate  curves  shown  in  figures  4.19-4.22  indicate  that  nonequilibrium  effects  are  most 
significant  for  weak  (but  finite)  strength  waves.  For  vanishingly  weak  waves  the  jet  reaction  rate 
increases  and  this  effect  becomes  more  abrupt  as  the  Mach  number  increases.  At  higher  Mach 
numbers  (Figures  4.19  &  4.20)  the  reaction  rate  in  the  jet  is  depressed  over  a  wide  range  of  shock 
angles. 

End  points  of  the  various  solution  curves  correspond  to  conditions  at  which  the  flow  downstream 
of  shock  5  becomes  subsonic  with  respect  to  the  frozen  speed  of  sound.  The  differing  end  points 
of  the  frozen,  nonequilibrium  and  equilibrium  solutions  are  a  manifestation  of  the  behavior  of  the 
sonic  line  that  was  discussed  in  section  4.4.2.  The  limited  extent  of  the  nonequilibrium  solution 
in  figure  4.22  indicates  that  finite  width  subsonic  regions  become  important  features  of  the  flow  at 
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Figure  4.23:  Variation  of  type  IV  jet  shock  strengths  with  impinging  shock  angle,  /?i. 
Free  stream  conditions  for  the  frozen  solutions  are  Poo  =  0.010,  Ooo  =0;  Pi  =  8°,  15°,  20°,  40°. 
The  shock  angles  and  flow  deflection  angles  shown  are  to  scale. 


lower  Mach  numbers. 


4.8.5  Interpretation  of  the  mechanisms 

The  thermochemical  effects  that  are  observed  in  figures  4.18-4,22  may  be  understood  from  the 
variations  of  the  jet  shock  strengths  that  are  demonstrated  in  figure  4.23.  At  the  intermediate  shock 
angles,  Pi  =  15°  and  Pi  =  20°,  the  strengths  of  the  jet  shocks  (1,3  &  5)  are  balanced  to  produce  the 
maximum  compression  and  heat  transfer  intensification.  For  very  weak  impinging  waves,  pi  —  8°, 
we  see  that  the  shock  connecting  the  two  A-points  becomes  disproportionately  strong.  In  the  case  of 
a  strong  impinging  wave,  pi  =  40°,  the  compression  across  the  incident  wave  dominates.  Maximum 
heating  results  therefore  from  a  balancing  of  the  jet  shock  strengths.  The  key  conclusion  is  that  this 
balance  simultaneously  reduces  the  relative  influence  of  equilibrium  real  gas  effects  and  depresses 
the  jet  reaction  rate.  This  follows  from  consideration  of  figures  4.5  and  4.6  which  demonstrate  the 
limited  extent  and  rate  of  the  dissociation  reaction  downstream  of  moderate  strength  shock  waves. 
Real  gas  effects  become  significant  at  lower  Mach  numbers  and  away  from  the  peak  heating  shock 
angle  when  at  least  one  of  the  jet  shocks  causes  significant  dissociation. 
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4.9  Comparison  with  existing  models 

Edney  [8]  presented  the  type  IV  model  data  in  terms  of  the  jet  pressure  ratio  because  this  quantity 
is  sensitive  to  variations  in  the  Mach  number  and  ratio  of  specific  heats.  For  the  current  model, 
density  and  heat  transfer  ratios  were  better  measures  of  the  influence  of  real  gas  effects.  Provided 
that  the  jet  Mach  number  remains  high  we  can  approximate  the  intensification  of  the  stagnation 
pressure  by, 

P6o  Pzul  _  .  ^  ~  A'S 

—  = - ^  -  P^PapijT-- 

Pe  Poo  W 00 

Using  this  form,  figure  1.7  compares  the  current  frozen  solution  with  the  results  of  Edney  [8].  The 
real  gas  solution  has  been  verified  against  independent  calculations  for  a  single  oblique  shock  wave. 

The  curve  for  7  =  1.2  in  figure  1.7  demonstrates  Edney ’s  prediction  of  greatly  increased  heat 
transfer  at  high  enthalpies.  This  value  is  typical  for  stagnation  temperatures  in  the  range  8000- 
12000K  that  are  studied  here.  Comparing  figure  1.7  with  figure  4.20  we  see  that  the  variable  7 
model  considerably  over-estimates  the  difference  between  the  frozen  and  equilibrium  solutions  of  the 
IDG  model  at  the  same  Mach  number.  This  occurs  because  the  density  ratio  across  all  shock  waves 
is  increased  uniformly  as  7  is  reduced.  This  includes  the  jet  shock  waves  that  determine  the  heat 
transfer  intensification  and  across  these  waves  the  influence  of  real  gas  effects  is  limited. 


4.10  Comparison  with  experimental  data 

Tables  4.3-4. 5  contain  the  detailed  predictions  of  the  model  for  the  free-stream  conditions  that  were 
used  in  the  experimental  portion  of  this  work.  The  parameters  for  the  computation,  Poo,  Pd, 

Uoo  and  /?!,  are  given  in  tables  B.l  and  B.5.  No  other  adjustable  parameters  appear  in  the  model, 

other  than  the  indeterminate  geometrical  terms  that  are  estimated  below .  _ _ 

The  heat  transfer  rate  is  predicted  by  the  model  in  the  universal  form,  ^  ' 

Both  the  undisturbed  body  flow  and  the  impinging  jet  are  planar  so  that,  =  Cij.f  =  0-70 

(Equation  4.45).  Equation  4.48  gives  =  1  and  <2,.*  =  0-26  •  During  discussion  of  the  experi¬ 

mental  results  it  was  noted  that  the  maximum  heat  transfer  occurs  when  the  inverted  A-shock  lies 
close  to  the  surface  of  the  body  (Figures  3.3  &  3.4).  Noting  this  we  propose  the  following  as  an 
unbiased  scaling  of  the  jet  width  to  body  diameter  ratio; 

^  =  ^sin(|^3 -I- -  ^21).  (4.59) 

This  relationship  scales  the  jet  width  with  the  undisturbed  shock  standoff  distance,  x?  ^nd  the 
divergence  angle  of  the  oblique  wave  and  shear  layer  at  the  A-shock  (Figure  4.10).  All  of  the 
quantities  appearing  in  equation  4,59  may  be  predicted  apriori.  The  undisturbed  shock  standoff 
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Frozen 

Nonequilibrium 

Equilibrium 

A 

14.25° 

14.25° 

14.25° 

-87.1° 

-87.1° 

-87.1° 

h 

-33.1° 

-33.1° 

-33.1° 

61 

10.2° 

10.2° 

10.2° 

62 

-15.5° 

-15.5° 

-15.5° 

£6 

21.5 

21.5 

21.5 

Pe 

das  f  dae 
dofy  d^ 

- 

- 

- 

/  .2  ^ 

ij&i  1  Ax  ^''^body^^i>ody  \ 

gbody  D  / 

4.63 

4.63 

4.63 

x/D 

0.218 

0.218 

0.218 

1^3  +  ^1  -  ^2] 

7.4° 

7.4° 

7.4° 

Qiei 

tjbodv 

14.1 

14.1 

14.1 

Table  4.3:  Model  predictions  for  condition  A. 


Frozen 

Nonequilibrium 

Equilibrium 

Pi 

15° 

15° 

15° 

P2 

-84.7° 

-86.3° 

-86.2° 

Pz 

-43.5° 

-45.1° 

-44.9° 

Si 

8.12° 

8.12° 

8.12° 

S2 

-24.3° 

-25.5° 

-25.7° 

h 

Pe 

7.05 

6.53 

7.13 

das  j  dae 

d^  rio  dx 

- 

0.0089 

- 

f  c'^  C  \ 

jj&i  I  Ax  ^^body^'^body  \ 

2.66 

2.54 

2.75 

ibody  1  D  C?  ) 

xlD 

0.186 

0.186 

0.186 

\Pz  +  ^1  —  ^2  1 

11.1° 

11.5° 

11.1° 

Qjei 

7.16 

6.73 

7.41 

^bodv 

Table  4.4:  Model  predictions  for  condition  B. 


distances  were  obtained  from  table  B.6 

Examination  of  the  jet  divergence  angle,  |/?3  +  (5i  —  62 1,  in  table  4,5  reveals  an  additional  influence 
of  finite  rate  effects  whereby  the  decreased  density  of  the  jet  causes  a  divergence  of  the  shock  and 
shear  layer.  This  translates  into  a  further  reduction  in  the  heat  flux  that  is  not  encompassed  by  the 


universal  form, 


4body 


Aa? 


Although  the  calculations  were  based  on  the  measured  incident  shock  angles,  /?i,  the  flow  deflec¬ 
tion  angles  computed  from  them  are  somewhat  larger  than  the  actual  shock  generator  angle  of  6°. 
The  discrepancy  is  consistent  with  the  discussion  of  nozzle  flow  nonuniformity  and  Mach  number 
over-prediction  in  appendix  B.  Observe  the  decline  of  the  discrepancy  with  increasing  stagnation 
enthalpy  that  is  consistent  with  figures  B.1-B.3. 
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Frozen 

Nonequilibrium 

Equilibrium 

Pi 

le" 

16° 

16° 

P2 

-ss.r 

-86.1° 

-85.7° 

/?3 

-48.4'’ 

-51.6° 

-49.3° 

^1 

7.31° 

7.31° 

^2 

-27.2° 

-29.2° 

-31.5° 

4.12 

3.59 

5.61 

dor 5  /  dae 

dx  r'^  dx  f-wi 

- 

0.119 

- 

/  c  \ 

/  Ao;  \ 

ibody  yo  J 

XID 

2.03 

0.173 

1.82 

0.173 

2.52 

0.173 

IPs  +  ^1  -  ^2 

IS.O"* 

15.1° 

10.5° 

iiei 

^bodv 

5.07 

4.37 

7.23 

Table  4.5:  Model  predictions  for  condition  C. 

4.10.1  Heat  flux  data 

At  condition  A  the  dissociation  levels  are  negligible  and  so  the  frozen,  nonequilibrium  and  equilibrium 
solutions  contained  in  table  4.3  are  indistinguishable.  The  peak  heat  transfer  rate  agrees  remarkably 
with  the  heat  transfer  measurements  that  are  summarized  in  table  4.6. 

Nonequilibrium  effects  are  most  pronounced  for  condition  B  where  the  reaction  rate  in  the  su¬ 
personic  jet  is  two  orders  of  magnitude  less  than  that  behind  the  undisturbed  bow  shock  (Table  4.4). 
The  influence  of  this  on  the  measurable  quantities  is  only  slight  since  the  dissociation  levels  remain 
small.  Despite  the  limited  number  of  data  points  obtained  at  this  condition,  the  trends  predicted 
by  the  model  are  confirmed  by  the  data  in  table  4.6. 

The  effect  of  dissociation  chemistry  is  most  obvious  at  condition  C  (Table  4.5).  Significant  dis¬ 
sociation  occurs  and  this  is  reflected  in  the  differing  peak  density  and  heat  transfer  ratios.  Chemical 
nonequilibrium  produces  a  decrease  in  these  quantities  however  the  depression  of  the  jet  reaction  rate 
is  only  moderate.  Since  the  reaction  rate  parameter  is  large  for  the  undisturbed  flow  (A  =  21.),  the 
jet  should  remain  close  to  equilibrium.  Table  4.6  indicates  experimental  heat  transfer  rates  that  are 
consistent  with  this  interpretation.  Although  the  model  appears  to  predict  the  heating  rates  for  the 
current  conditions  quite  reliably,  further  experiments  spanning  a  range  of  reaction  rate  parameters 
would  be  required  to  discern  the  nonequilibrium  behavior. 

4.10.2  Density  data 

The  peak  interferometric  fringe  shift  for  the  type  IV  flow  was  measured  from  the  holograms  shown  in 
figures  C.8,  C.15  and  C.34.  Table  4.7  lists  the  ratios  of  the  peak  fringe  shifts  relative  to  the  stagnation 
point  values  without  shock  impingement.  Estimates  of  the  error  include  only  the  precision  of  the 
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Test  condition 

A 

B 

C 

Peak  measured  4^ 

Qbody 

Instantaneous 

15.2 

10.8 

7.9 

100  jjLS  averaged 

13.5 

7.6 

6.0 

Predicted  4^ 

^body 

Frozen 

14.1 

7.2 

5.1 

Nonequilibrium 

14.1 

6.7 

4.4 

Equilibrium 

14.1 

7.4 

7.2 

Table  4.6:  Comparison  of  measured  and  predicted  heat  transfer  intensification  for 
TYPE  IV  flows  with  CYLINDRICAL  MODEL.  The  experimental  heat  transfer  rates  are  obtained 
from  figures  3. 6-3. 8 


Test  condition 

A 

B 

C 

Undisturbed  fringe  shift 

13  ±  1 

10.5  ±  1 

9±1 

Type  IV  fringe  shift 

55  ±2 

29.5  zb  1 

27zbl 

Ratio 

4.2±0.5 

2.81  ±0.4 

3.0  ±0.4 

Predicted  ^ 

Pe 

Frozen 

21.5 

7.05 

4.12 

Nonequilibrium 

- 

6.53 

3.59 

Equilibrium 

7.13 

5.61 

Table  4.7:  Comparison  of  measured  fringe  shift  and  predicted  stagnation  density  intensification  for 
type  IV  flows  with  cylindrical  model. 


measurement.  Additional  systematic  error  was  introduced  through  line  of  sight  integration  of  span- 
wise  variations  of  the  jet  structure.  The  peak  fringe  shift  ratios  are  therefore  smaller  than  the  density 
ratios  predicted  by  the  model. 

Figures  C.48  and  C.50  illustrate  the  flow  over  a  blunt  ended  cylinder  at  conditions  A  and  C. 
Since  the  chord  length  across  the  face  of  the  model  was  similar  in  both  cases,  the  peak  fringe  shifts 
of  17  zb  1  and  8  db  1  can  be  meaningfully  compared.  Weighted  by  the  free-stream  densities,  the  ratio 
of  these  values  is  3.0.  Although  three-dimensional  effects  were  still  present,  this  ratio  approaches 
that  of  the  peak  densities  predicted  by  the  model  in  the  equilibrium  limit;  3.8 . 


4.11  Dissociation  and  recombination  rate  parameters  for 
the  type  IV  jet. 

As  a  final  check  on  the  validity  of  the  assumptions  of  the  model  we  consider  the  dissociation  and 
recombination  rate  parameters  for  the  type  IV  jet.  The  definitions  of  these  parameters  follow  from 
equations  4.25  and  4.50  where  the  free-stream  conditions  are  now  taken  to  be  those  downstream  of 
the  jet  shock  system  (Figure  4.17).  It  is  convenient  to  represent  these  parameters  as  ratios  of  the 
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form, 

As  D  _  dae/dx5 
Aco  Ax  dcxe/dxoo 

where  the  characteristic  reaction  lengths  have  been  renormalized  according  to  equation  4.40.  The 
reaction  rates,  dae/dx^  and  da^dioo,  along  with  the  frozen  shock  density  ratios,  and  pf^,  are 
obtained  from  the  type  IV  jet  solution  given  in  sections  4.5  and  4.6. 

The  ratio  of  the  recombination  rate  parameters  follows  from  equation  4.50; 


A' 


-P5P3P1 


fUl 

Ph 


(4.60) 


Ss  A>C2,e. 

Sqo 


^OO 


Ac 


-P5P3Pl , 


(4.61) 


where  as  =  pT  /pd^  and  a^  =  I pd^o  are  determined  from  the  type  IV 

jet  solution.  The  temperature  ratios  in  these  expressions  are  given  by  equation  4.24. 

Figure  4.24  shows  the  jet  dissociation  rate  parameter  for  the  free-stream  conditions  that  cor¬ 
respond  to  figure  4.20.  Increases  in  the  dissociation  rate  due  to  the  higher  density  at  the  jet  im¬ 
pingement  point  are  offset  by  the  ratio  A(r/D.  The  net  effect  is  that  the  jet  shock  dissociation  rate 
parameter  is  similar  to  the  undisturbed  bow  shock  rate  parameter.  Due  to  relaxation  that  occurs 
in  the  jet,  the  equilibrium  curve  lies  below  the  nonequilibrium  curve.  This  also  explains  the  strong 
dip  in  the  equilibrium  solution  for  weak  impinging  waves  where  the  unbalanced  oblique  jet  shock 
system  causes  significant  dissociation. 

The  recombination  rate  parameter  ratios  are  observed  to  vary  with  the  square  of  the  dissociation 
rate  parameter  ratios.  This  is  consistent  with  the  scaling  of  the  recombination  rate  with  the  square 
of  the  density.  We  conclude  that  the  jet  boundary  layer  lies  significantly  closer  to  equilibrium  than 
that  on  the  undisturbed  blunt  body.  For  conditions  where  Eqo  ^  1  but  ^  this  may  produce  an 
increase  in  the  heat  transfer  rate  if  the  body  surface  is  non-cat alytic  to  the  recombination  reaction. 


4.12  Summary 

The  current  experiments  exposed  deficiencies  in  existing  variable  j  models  of  the  shock  impinge¬ 
ment  phenomenon.  These  models  over-predict  the  experimentally  observed  heat  transfer  rates  and 
densities  at  high  enthalpy.  Improved  models  of  the  inviscid  jet  shock  system  and  the  impingement 
point  boundary  layer  were  introduced  in  order  to  explain  the  discrepancies. 

The  main  real  gas  effects  on  the  inviscid  jet  shock  system  arise  because  of  the  widely  differing 
shock  strengths  produced  at  the  A-points.  Rapid  dissociation  occurs  downstream  of  the  strong 
waves  whereas  the  extents  and  rates  of  dissociation  for  the  weak  jet  shocks  are  depressed.  It  was 
further  demonstrated  that  peak  heating  occurs  when  the  strengths  of  the  oblique  shock  waves  in  the 
supersonic  jet  are  balanced.  This  minimizes  the  entropy  rise  along  the  stagnation  streamline.  These 
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Figure  4.24:  Reaction  rate  parameters  for  the  type  IV  jet  stagnation  point.  Upper  plot 
is  the  ratio  of  the  dissociation  rate  parameters  in  the  form,  ^  Lower  plot  is  the  ratio  of  the 

recombination  rate  parameters,  ‘  “  —  curves  are  the  equilibrium  solutions, 

curves  are  the  nonequilibrium  solutions.  The  free-stream  conditions  are  Poo  ==  0.010,  Hq^  ~  0.6, 
PdoQ  “lx  10  ,0^00  —  0. 
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effects  combine  to  greatly  reduce  the  influence  of  thermochemistry  at  the  peak  heating  condition. 
Real  gas  effects  were  shown  to  be  important  when  the  balance  of  the  oblique  shock  waves  is  disturbed. 
This  occurs  for  impinging  shock  strengths  greater  or  less  than  that  which  produces  maximum  heating 
and  for  lower  free  stream  Mach  numbers. 

A  nonequilibrium  approximation  was  introduced  that  applies  when  the  oblique  jet  shocks  are 
weak  with  respect  to  the  undisturbed  bow  shock.  Within  the  scope  of  the  approximation,  non¬ 
monotonic  behavior  with  the  reaction  rate  is  predicted. 

The  model  successfully  reproduced  the  experimentally  observed  heat  transfer  rates.  Accurate 
predictions  of  the  heat  transfer  rates  should  be  possible  for  conditions  where  turbulent  diffusion  of 
the  supersonic  jet  is  unimportant.  Further  experimentation  is  necessary  to  resolve  the  prediction  of 
non-monotonic  behavior  with  the  reaction  rate  parameter. 
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Chapter  5  Conclusions 


A  series  of  experiments  was  conducted  in  the  T5  free  piston  driver  shock  tunnel  to  determine  the 
quantitative  effects  of  high  enthalpy  thermochemistry  on  the  shock  impingement  heating  problem. 
By  studying  the  nominally  two-dimensional  mean  flow  about  a  cylinder  with  a  coplanar  impinging 
shock  wave  the  topology  of  the  problem  was  simplified  so  that  it  could  be  observed  directly  by 
holographic  interferometry.  Interferograms  obtained  for  three  different  test  conditions  and  a  range 
of  shock  impingement  locations  demonstrated  the  controlling  influence  of  the  inertia  of  the  supersonic 
jet  and  indicated  that  existing  predictions  of  greatly  increased  stagnation  density  at  high  enthalpy 
were  not  realized.  Time  resolved  heat  transfer  measurements  confirmed  this  observation  and  agreed 
accurately  with  existing  low  enthalpy  experimental  investigations. 

Temporal  fluctuations  of  the  impinging  jet  were  observed  at  all  test  conditions  and  Strouhal 
numbers  were  estimated  from  the  power  spectra  of  the  heat  transfer  data.  Variations  of  the  tunnel 
test  conditions  during  the  test  time  limited  the  quality  of  the  spectra  at  high  enthalpies.  Whilst  it 
was  not  possible  to  conclusively  discount  the  effect  of  free-stream  perturbations,  the  data  exhibited 
features  that  have  been  observed  in  recent  numerical  studies.  The  Strouhal  numbers  were  consistent 
with  the  limited  existing  data. 

New  instrumentation  developed  for  this  project  has  allowed  critical  assessment  of  the  flow  quality 
produced  by  the  T5  facility.  Measurements  of  the  flow  over  blunt  bodies  without  shock  impingement 
agreed  with  established  theories  within  experimental  uncertainty.  Flow  contamination  that  was 
detectable  with  this  instrumentation  did  not  become  apparent  until  times  beyond  the  period  of  test 
flow  that  was  used  to  obtain  the  data.  This  confirmed  the  validity  of  the  current  measurements. 

In  order  to  reconcile  the  experimental  observations  with  existing  analyses  the  influence  of  ther¬ 
mochemistry  on  the  jet  shock  system  was  modeled  using  the  ideal  dissociating  gas  approximation. 
Careful  nondimensionalization  of  the  problem  led  to  a  tractable  formulation  and  identification  of  the 
important  parameters.  The  results  were  used  to  demonstrate  generic  influences  of  real  gas  effects 
on  the  shock  impingement  problem.  The  central  conclusion  was  that  peak  heating  arises  from  a 
balancing  of  the  jet  shock  strengths.  This  balance  was  demonstrated  to  simultaneously  reduce  the 
relative  influence  of  equilibrium  real  gas  effects  and  depress  the  jet  reaction  rate.  Real  gas  effects 
were  shown  to  be  important  at  lower  Mach  numbers  (<  7.5)  and  for  impinging  shock  angles  away 
from  that  which  produces  peak  heating.  The  model  accurately  reproduced  the  experimentally  ob¬ 
served  heat  transfer  rates  using  apriori  estimated  parameters.  Interferometric  comparisons  were 
hindered  by  variations  in  the  impingement  location  across  the  span  of  the  model.  Relative  trends  of 
the  interferometric  data  at  the  different  conditions  were  consistent  with  the  model  predictions. 


92 


A  non-equilibrium  approximation  was  introduced  that  applies  under  conditions  where  the  jet 
shocks  remain  weak  with  respect  to  the  undisturbed  bow  shock.  Within  the  scope  of  the  approxi¬ 
mation  non-monotonic  behavior  was  predicted  with  the  reaction  rate  parameter.  Discrimination  of 
finite  rate  effects  w^as  not  possible  in  the  current  investigation  since  the  reaction  rate  was  not  varied 
as  an  independent  parameter. 
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Appendix  A  Review  of  shock  impingement  flow  field 
models 

A.l  Introduction 

It  is  appropriate  to  include  a  brief  review  of  the  seminal  observations  of  Edney  [8]  since  these  form 
the  basis  of  what  follows.  Schlieren  cinematography  of  shock  wave  impingement  flow  fields  allowed 
Edney  to  categorize  six  distinct  flow  regimes  known  as  types  I- VI.  These  six  interaction  patterns 
arise  for  different  locations  of  the  impinging  shock  wave  with  respect  to  the  blunt  body,  as  indicated 
in  figure  1.3.  The  sequence  of  interactions,  obtained  by  Edney  for  a  hemispherical  body,  is  illustrated 
in  figures  A.l  and  A. 2.  These  flow  fields  were  observed  for  a  variety  of  model  geometries  and  so  the 
classification  appears  to  be  universal. 

A. 2  Shock  loci  methods  for  a  perfect  gas 

The  flow  in  smooth  regions  of  the  shock  impingement  flow  field  and  across  shock  discontinuities  is, 
in  principle,  understood.  The  predominant  concern  is  then  points  where  surfaces  of  discontinuity 
intersect  and  here  local  analysis  about  intersection  points  in  the  pressure  -  flow  deflection  angle 
[p-6)  plane  provides  a  useful  tool.  These  techniques  were  developed  in  the  study  of  Mach  reflection 
problems  and  reviews  are  given  by  Courant  h  Friedrichs  [6]  and  Hornung  [20]  who  reviews  recent 
developments.  In  this  method  the  solution  is  represented  graphically  by  the  intersection  in  the  p-b 
plane  of  the  loci  of  possible  downstream  states  for  the  interacting  waves. 

In  the  absence  of  viscous  effects,  the  usual  Mach  reflection  problems  may  be  mirrored  about  the 
surface  of  reflection.  They  are  therefore  degenerate  symmetrical  cases  in  which  two  shock  surfaces 
of  opposite  families  interact.  Although  the  asymmetric  case  has  received  relatively  little  attention, 
many  of  the  ideas  in  the  literature  carry  over  in  modified  form  to  the  shock  on  leading  edge  problem. 
Edney  [8]  used  these  methods  to  interpret  the  shock  impingement  flow  fields  and  justify  the  type  I-VI 
classification.  The  key  conclusion  drawn  from  this  analysis  was  the  role  of  the  three  shock  solutions 
(or  A-shocks).  Our  approach  here  is  to  briefly  review  the  existing  methods  before  extending  them 
to  study  the  influence  of  real  gas  effects. 

For  a  perfect  gas  the  mapping  of  an  oblique  shock  wave  into  the  p-b  plane  is  determined  by 
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Figure  A. 2:  Interference  types  IV-VI  for  supersonic  flow  about  a  hemisphere;  Moc  =  4.6  (Edney  [8]). 
Sequence  of  interactions  is  left-right,  top-bottom. 
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conservation  of  mass,  momentum  and  energy  normal  to  the  shock  front  to  be, 

.  _  (7  +  l)Mg,sin^/3 

^  (7  -  l)M‘^  sin^  /?  —  2  ’ 

"-14.  P~  ^ 

^  sin^  /?  ’ 

tan  ,8  _  , 

tan(,5  —  6) 

where  p  and  p  are  the  density  and  pressure  ratios  across  the  shock.  The  free-stream  Mach  number 
is  Moc  and  7  is  the  ratio  of  specific  heats.  For  a  given  free-stream  condition  the  above  equations 
determine  the  locus  of  possible  downstream  states.  This  forms  a  curve  in  the  p-6  plane  that  is 
parameterized  in  terms  of  the  shock  angle,  /?. 

A. 3  The  type  I- VI  classification 

A. 3.1  Type  I 

Type  I  interference  is  observed  in  the  first  frame  of  figure  A.l  and  this  is  the  well  known  regular 
solution  for  the  interaction  of  weak  oblique  shock  waves.  The  type  I  flow  field  is  illustrated  schemat¬ 
ically  in  figure  A.3  along  with  its  mapping  into  the  p-6  plane.  Enhanced  heating  occurs  as  a  result 
of  the  shock  wave  boundary  layer  interaction  that  occurs  at  the  body  surface.  Results  from  the 
literature  (Delery  &  Marvin  [7])  may  be  expected  to  yield  reasonable  estimates  of  the  heating  rates. 

A.3.2  Type  II 

As  the  impingement  point  moves  upward  towards  the  lower  sonic  line  a  type  II  interaction  is  observed 
in  the  second  frame  of  figure  A.l.  In  this  configuration  the  oblique  shock  loci  (Figure  A.4)  do  not 
intersect  and  a  subsonic  region  is  observed  behind  a  Mach  stem  connecting  the  two  mutually  inverted 
A-pattern  shock  intersection  points.  The  transition  mechanism  from  type  I  is  similar  to  the  Mach 
reflection  problem  however  the  strengths  of  the  interacting  waves  are  unequal.  Edney’s  criterion  for 
transition  was  equivalent  to  the  maximum  deflection  condition  in  Mach  reflection  (Hornung  [20]). 
Mechanisms  for  enhanced  heating  in  type  II  flows  are  similar  to  those  for  type  I  flows  however  the 
shock  that  interacts  with  the  boundary  layer  is  generated  at  the  upper  A-point  (Figure  A.4). 

A. 3. 3  Type  III 

When  the  Mach  stem  merges  with  the  subsonic  portion  of  the  bow  shock  ahead  of  the  body  a  type  III 
flow  results  (Figure  A.l,  frames  3  &  4).  During  the  type  II-III  transition  the  locus  representing  the 


'ressure  ratio 


-40  -20  0  20  40 

Deflection  angle  (deg) 


Figure  A. 3:  Schematic  of  type  I  interference  flow  and  its  representation  in  the  p-6  plane. 
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Figure  A. 5:  Schematic  of  type  III  interference  flow  and  its  representation  in  the  p-S  plane. 

upper  A-point  in  Figure  A. 4  becomes  vanishingly  small.  The  influence  of  the  impinging  wave  is  now 
global  since  it  impinges  on  the  subsonic  portion  of  the  bow  wave.  Both  the  radius  of  curvature 
and  the  standoff  distance  of  the  asymmetrical  upper  portion  of  the  bow  shock  increase  to  match 
the  A-point  solution  (Figure  A. 5).  In  frame  5  of  figure  A.l  the  shear  layer  that  is  generated  at  the 
shock  impingement  point  is  seen  to  deflect  where  it  meets  the  body.  The  shear  layer  is  assumed  to 
attach  to  the  surface  of  the  body  and  a  compression  wave  is  generated  that  overtakes  and  steepens 
the  third  wave  produced  at  the  impingement  point.  As  this  secondary  interaction  becomes  stronger 
(Figure  A.l,  frame  6)  the  transition  to  type  IV  begins.  Transition  to  type  IV  is  presumed  to  occur  at 
the  detachment  angle  for  the  compression  wave.  Observe  that  a  parallel  exists  between  this  process 
and  the  mechanism  of  complex  Mach  reflection  of  a  propagating  shock  from  the  surface  of  a  wedge. 

A. 3. 4  Type  IV 

The  supersonic  jet  that  characterizes  the  type  IV  flow  field  is  observed  to  form  in  the  first  frame 
of  figure  A. 2.  An  extensive  description  of  the  type  IV  flow  field  is  contained  in  section  1.2.  Type 
IV  interactions  produce  severe  pressure  loads  and  heat  transfer  rates  where  the  supersonic  jet  im- 
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Figure  A. 6:  Schematic  of  type  V  interference  flow. 


pinges  on  the  body.  The  width  of  the  type  IV  jet  reduces  as  the  impingement  point  moves  upward 
(Figure  A. 2,  frames  2“-4)  and  a  point  of  inflection  is  induced  on  the  bow  shock  below  the  interaction 
region.  At  the  highest  impingement  point  the  supersonic  jet  is  barely  distinguishable. 

A. 3. 5  Type  V 

The  width  of  the  supersonic  jet  decreases  rapidly  as  the  impingement  point  moves  toward  the  upper 
sonic  line.  Although  the  jet  structure  was  not  resolvable  Edney  [8]  assumed  that  this  remains 
unchanged  from  the  type  IV  solution.  No  other  local  solution  appears  to  exist.  Since  the  bow  shock 
angle  immediately  below  the  jet  is  fixed  by  the  local  solution,  a  point  of  inflection  is  induced  on  the 
wave  front  for  sufficiently  high  impingement  points  (Figure  A. 6).  The  type  V  flow  is  characterized 
by  the  appearance  of  a  third  A-point  that  causes  a  weak  wave  to  impinge  on  the  surface  of  the  body 
(Figure  A. 2,  frame  5). 

Type  V  flows  are  observed  more  clearly  in  shock-on-fln  impingement  problems  where  the  curva¬ 
ture  of  the  bow  shock  in  the  plane  of  symmetry  is  greatly  reduced. 

A. 3. 6  Type  VI 

For  impingement  points  lying  well  above  the  upper  sonic  line  a  type  VI  flow  is  observed  (Figure  A. 2, 
frame  6).  This  is  analogous  to  the  conventional  overtaking  intersection  of  oblique  shock  waves  with 
a  weak  reflected  expansion  fan.  The  representation  of  the  type  VI  flow  in  the  p-b  plane  is  shown  in 
figure  A. 7.  Relatively  little  is  known  about  the  transition  mechanism  for  this  overtaking  intersection 
in  comparison  with  the  body  of  literature  that  exists  for  the  Mach  reflection  problem.  Edney  used 
the  limit  of  tangential  intersection  of  the  p-b  loci  as  the  transition  mechanism  between  regimes  V 
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Figure  A. 7:  Schematic  of  type  VI  interference  flow  and  its  representation  in  the  p-6  plane. 


and  VL  Minor  heating  occurs  where  the  reflected  expansion  interacts  with  the  boundary  layer.  The 
interaction  occurs  downstream  of  the  sonic  line  so  that  flow  over  the  upstream  part  of  the  body  is 
equivalent  to  that  for  a  blunt  body  immersed  in  the  rotated  flow  behind  the  impinging  shock. 
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Appendix  B  Test  section  conditions 
B.l  Summary 

Three  different  test  conditions  were  used  in  the  experiments  described  here  and  these  were  chosen 
to  highlight  the  influence  of  real  gas  effects  on  the  shock  impingement  phenomenon.  Speciflc  values 
for  the  flow  parameters  are  given  in  table  B.l.  The  first  case,  A,  represents  a  low  enthalpy  control 
condition  where  the  only  real  gas  effects  were  a  partial  excitation  of  the  vibrational  modes  in  the 
stagnant  regions  of  the  flow.  Control  of  the  T5  piston  dynamics  was  optimized  for  high  enthalpy  flows 
and  prevented  the  operation  of  the  facility  at  lower  enthalpies  with  the  current  120  kg  piston.  Case 
C  represents  a  high  enthalpy  condition  that  was  chosen  to  highlight  the  influence  of  thermochemistry 
on  the  problem.  Case  B  is  an  intermediate  enthalpy  condition  that  was  designed  to  aid  interpretation 
of  the  differences  in  the  flow  mechanisms  for  cases  A  and  C. 


B.2  Prediction  of  the  test  section  conditions 

T5  was  operated  as  a  reflected  mode  shock  tunnel  whereby  the  reservoir  conditions  upstream  of  the 
nozzle  were  generated  by  the  reflection  of  the  incident  shock  wave  from  the  end  wall  of  the  shock 
tube.  The  volume  of  double  shock  heated  gas  was  then  expanded  through  a  contoured  nozzle  into 
the  test  section.  Prediction  of  the  test  section  conditions  was  based  on  a  two-step  process  whereby 
the  reservoir  conditions  and  nozzle  flow  were  computed  sequentially. 

Given  the  initial  filling  conditions  of  the  shock  tube  and  the  measured  velocity  of  the  incident 
shock  wave,  the  reservoir  condition  was  calculated  using  an  inviscid  equilibrium  chemistry  model  of 
the  incident  and  reflected  waves.  The  stagnant  reservoir  state  was  then  expanded  isentropically  to 
the  measured  reservoir  pressure  in  order  to  account  for  subsequent  weak  waves  reflected  from  the 
driver  gas  interface,  propagated  influences  that  arose  from  the  unsteady  motion  of  the  piston  and 
the  expansion  of  the  reservoir  volume  supplying  the  nozzle.  This  computation  was  performed  using 
the  ESTC  code  described  by  Macintosh  [31]  where  the  reader  may  obtain  further  details  of  the  tech¬ 
nique.  The  thermodynamics  of  the  system  were  modeled  using  harmonic  oscillator  approximations 
at  temperatures  below  5000K  and  polynomial  fits  above  this  threshold. 

Flow  through  the  nozzle  was  calculated  using  the  2D  supersonic  inviscid  reacting  nozzle  flow 
code  (SURF)  developed  by  Rein  [38].  The  thermodynamic  models  were  similar  to  those  described 
above  with  the  addition  that  the  nonequilibrium  gas  was  described  as  a  mixture  of  thermally  perfect 
components.  Reaction  rates  were  described  using  Arrhenius  forms  and  vibrational  excitation  was 
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Test  condition 

A 

B 

C 

Shock  tube 
conditions 

Fill  pressure 
(kPa  ±0.5kPa) 

75 

32.5 

20 

Fill  temperature 
(K  ±4K) 

297 

297 

297 

Shock  speed 

1890 

3500 

4360 

(m/s) 

±100 

±245 

±255 

Nozzle 

Pressure 

14.6 

22.3 

28.3 

reservoir 

(MPa) 

±0.7 

±1.3 

±1.9 

conditions 

Temperature 

(K) 

3210 

7550 

8960 

Enthalpy 

(MJ/kg) 

3.88 

12.0 

19.1 

Test  section 

Velocity  (m/s) 

2540 

4450 

5350 

conditions 

Density  (kg/m^) 

0.0218 

0.0155 

0.0157 

Pressure  (kPa) 

1.03 

5.48 

11.4 

N  concentration  (mole/kg) 

9.9  X  10-® 

6.9  X  10-1 

3.65  X  10° 

Test  section 

IDG  Poo 

7.32  X  10“^ 

1.79  X  10-"* 

2.54  X  10-3 

dimensionless 

IDG  Po„ 

0.102 

0.346 

0.560 

parameters 

IDG  ttoo 

1.4  X  10-® 

9.7  X  10-3 

0.051 

IDG  pd 

6.0  X  10® 

8.4  X  10® 

8.3  X  10® 

IDG  A  per  mm 

4.  X  10-1^ 

3.4  X  10-3 

0.525 

Re  per  mm 

5540 

1630 

1350 

Table  B.l:  Summary  of  the  free-stream  conditions.  Test  conditions  were  computed  on  the 
basis  of  measurements  of  the  initial  shock  tube  fill  pressure  and  temperature,  incident  shock  speed 
prior  to  reflection,  and  nozzle  reservoir  pressure  after  shock  reflection.  The  dimensionless  parameters 
for  the  ideal  dissociating  gas  (IDG)  model  are  defined  in  section  4.2.  Error  estimates  given  for  the  fill 
conditions  represent  the  accuracy  of  the  pressure  gauge  and  the  variation  of  the  ambient  temperature. 
Error  estimates  for  the  measured  shock  speed  and  reservoir  pressure  are  the  standard  deviation  of 
the  quantities  sampled  over  the  entire  sequence  of  shots.  The  error  in  the  computed  quantities  may 
be  inferred  from  these  estimates. 
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ESTC/SURF 

IDG/NENZF 

Shock  tube 
conditions 

Fill  pressure 
(kPa) 

25 

25 

Fill  temperature 

(K) 

300 

300 

Shock  speed 
(m/s) 

5600 

5600 

Nozzle 

reservoir 

Pressure 

(MPa) 

90. 

90.1 

conditions 

Temperature 

(K) 

11350 

11880 

Enthalpy 

(MJ/kg) 

32.9 

33.4 

Test  section  : 

Velocity  (m/s) 

6650 

6930 

conditions 

Density  (kg/m^) 

0.034 

0.033 

Pressure  (kPa) 

55.5 

50.0 

Atomic  N  (mole/kg) 

9.3 

17.5 

Test  section 

IDG  Poc 

3.7  X  lO-'-" 

3.1  X  lO--" 

dimensionless 

IDG  Pooo 

0.98 

0.99 

parameters 

IDG  Qoo 

0.130 

0.245 

IDG  pd 

3.82  X  10® 

3.94  X  10® 

Table  B.2:  Verification  of  ESTC/SURF  computations  by  comparison  with  simplified  IDG/NENZF 
model. 


modeled  using  a  simple  switch  from  frozen  to  equilibrium  at  3000K.  Verification  of  the  results  was 
obtained  by  comparison  with  a  simplified  calculation  using  the  IDG  model  for  the  reflected  shock 
and  the  one- dimensional  nonequilibrium  nozzle  flow  code  of  Lordi  et  al.  [30]  (Table  B.2). 


B.3  Variation  of  the  free-stream  conditions  across  the  test 
section 

Since  the  flow  in  the  nozzle  was  chemically  reacting,  a  contour  optimized  to  produce  axial  flow  at 
high  enthalpy  produced  nonuniform  flow  at  other  off  design  conditions.  The  variations  across  the 
test  section  were  obtained  at  the  nozzle  exit  plane  from  the  SURF  code  and  these  computations  are 
shown  in  figures  B.l,  B,2  and  B.3.  Spatial  nonuniformity  of  the  flow  increases  at  low  enthalpies  as 
expected.  Further  evolution  of  the  profile  must  be  expected  downstream  of  the  nozzle  exit  plane, 
however  it  was  not  possible  to  compute  this  with  the  current  implementation  of  the  SURF  code. 
The  conditions  given  in  table  B,1  were  based  on  conditions  35mm  above  the  nozzle  centerline.  This 
represented  the  approximate  location  of  the  model  for  the  type  IV  impingement  flow. 


Figure  B.l:  Flow  nonuniformity  in  the  nozzle  exit  plane  for  test  condition  A 
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Figure  B.2:  Flow  nonuniformity  in  the  nozzle  exit  plane  for  test  condition  B 
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Figure  B.3:  Flow  nonuniformity  in  the  nozzle  exit  plane  for  test  condition  C. 

B.4  Experimental  verification  of  the  flow  density 

B. 4.1  Interferometric  measurements  for  a  finite  aspect  ratio  cylinder 

The  accuracy  of  the  prediction  of  the  test  conditions  and  the  quality  of  the  flow  in  the  facility  may 
be  assessed  from  the  holograms  and  stagnation  point  heat  transfer  measurements  obtained  without 
shock  impingement.  Consider  first  the  stagnation  density  data  contained  in  figures  C.l,  C.14  and 

C. 20  that  correspond  to  test  conditions  A,  B  and  C  respectively.  The  variation  of  density  behind 
the  normal  portion  of  the  bow  shock  was  computed  for  condition  C  using  the  IDG  model  presented 
in  section  4.2.  The  curvature  of  the  shock  was  neglected  and  the  results  are  shown  in  figure  4.4. 
Observe  that  the  density  asymptotes  to  its  equilibrium  value  in  a  distance  that  is  small  compared 
to  the  shock  standoff  distance.  In  the  strong  shock  limit  for  a  perfect  gas  with  7  —  7/5  the  density 
rises  by  a  further  factor  of  1.08  between  the  bow  shock  and  the  stagnation  point.  We  extrapolate 
this  result  to  the  high  enthalpy  cases  with  sufficient  accuracy  for  the  present  comparison.  We  may 
therefore  compare  the  fringe  shift  calculated  from, 

=  {l.OSpoopeq  ”  Pco)  [^>3(1  “  C^eq)  +  KNO'eq]  y,  (BT) 

with  that  observed  in  the  interferograms  (Table  B.3).  Here  Ocaid^^'^  is  the  predicted  fringe  shift,  poo 
is  the  free-stream  density  predicted  by  SURF,  peq  is  the  equilibrium  IDG  post  shock  density  ratio, 
aeq  is  the  equilibrium  dissociated  mass  fraction,  and  Kn  are  the  Gladstone-Dale  constants  for 
molecular  and  atomic  nitrogen,  L  is  the  span  of  the  cylinder  and  A  is  the  wavelength  of  light.  The 
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Test  condition 

A 

B 

C 

Poc  (SURF  kg/m^) 

0.0218 

0.0155 

0.0157 

pjr  (IDG  model) 

6.61 

6.08 

5.61 

peq  (IDG  model) 

6.61 

8.59 

9.67 

afr  (IDG  model) 

1.4  X  10-® 

0.010 

0.051 

aeq  (IDG  model) 

2.0  X  10-® 

0.114 

0.285 

Calculated  fringe  shift  (fr) 

11.1 

7.17 

6.67 

Calculated  fringe  shift  (eq) 

11.1 

11.0 

13.2 

Measured  fringe  shift 

13.  ±  1 

10.5  ±1 

9.±1 

Table  B.3:  INTERFEROMETRIC  VERIFICATION  OF  THE  STAGNATION  DENSITY.  Physical  constants 

nsed  were;  L  =  0.180m,  A  =  532nm,  Kj^  —  3.1  x  10“^  Kn2  —  2.45  x  10  m^kg 

slight  effect  of  variation  of  the  refractivity  of  the  gas  mixture  at  the  free-stream  conditions  has  been 
neglected. 

Table  B.3  shows  that  at  low  enthalpy  the  stagnation  density  is  under-predicted.  This  was  prob¬ 
ably  due  to  the  strong  variation  in  density  and  pressure,  and  hence  Mach  number,  across  the  nozzle 
exit.  At  an  enthalpy  of  12MJ/kg  the  measured  stagnation  density  agrees  with  the  computation 
within  the  established  uncertainty.  Condition  C,  which  corresponds  to  an  enthalpy  of  19.1MJ/kg, 
demonstrated  a  stagnation  density  somewhat  lower  than  that  predicted  by  the  ESTC/SURF  com¬ 
putation.  The  uncertainty  of  ±1  fringe  indicated  in  table  B.3  is  based  on  the  difficulty  of  sub-fringe 
resolution  and  the  slight  deviations  from  infinite  fringe  conditions  that  are  observed  in  the  holo¬ 
grams.  Additional  computations  with  SURF  indicated  that  the  error  of  ±250  m/s  in  the  measured 
shock  speed  and  ±2  MPa  in  the  measured  reservoir  pressure  constitute  an  error  of  ±0.5  fringes  in  the 
calculated  value  for  high  enthalpy  condition  C.  Three-dimensional  flow  over  the  ends  of  the  cylinder 
introduces  additional  systematic  error. 

B.4.2  Abel  transform  interferometric  measurements  for  a  blunted  cylin¬ 
der 

For  axisymmetric  flow  fields  the  true  density  field  may  be  deconvolved  from  the  holographic  data  by 
use  of  Abel  transform  methods.  A  recursive,  direct  solution  algorithm  was  used  to  reconstruct  the 
density  field  in  a  stepwise  manner  from  the  edge  of  the  flow-field  (Emrich  et  al.  [10]).  Holographic 
data  was  available  for  the  flow  over  a  blunt  faced  cylinder  at  test  conditions  A  and  C  (Figures  C.4l 
&  C.49).  The  fringe  shift  was  measured  along  the  face  of  the  body  at  the  quarter  fringe  points  and 
cubic  splines  were  used  to  interpolate  between  them.  The  density  fields  deconvolved  from  the  spline 
fitted  data  are  shown  in  figures  B.4  and  B.5.  A  comparison  with  the  conditions  computed  using 
the  ESTC/SURF  codes  is  indicated  in  table  B.4.  The  measurements  and  predictions  agree  within 
experimental  uncertainty  if  potential  error  in  the  prediction  of  the  test  conditions  is  included.  These 
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Figure  B.4:  Abel  transformation  of  density  field  for  flow  over  the  blunt  faced  cylinder  at  test  condi¬ 
tion  A;  Shot  T5-819. 


Test  condition 

A 

C 

Poo  (SURF  kg/m^) 

0.0218 

0.0157 

pfr  (IDG  model) 

6.61 

5.61 

peq  (IDG  model) 

6.61 

9.67 

Calculated  ^pjr  (kg/m^) 

0.122 

0.072 

Calculated  Ape^  (kg/m^) 

0.122 

0.136 

Measured  Ap  (kg/m^) 

0.117  ±0.02 

0.110±0.02 

Table  B.4:  Interferometric  verification  of  the  stagnation  density  using  Abel  trans¬ 
form  METHODS.  The  calculated  values  are  obtained  from  the  ESTC/SURF  computations.  Error 
estimates  for  the  measured  values  are  based  on  a  sub-fringe  resolution  of  ±0.5  with  respect  to  the 
fringe  shift  at  the  stagnation  point. 


results  confirm  those  for  the  nominally  two-dimensional  cylinder  flow  discussed  above. 


B.4. 3  Density  variation  across  the  incident  shock  wave 

Further  information  on  the  free-stream  conditions  can  be  obtained  from  the  jump  in  density  that 
was  observed  across  the  incident  shock  wave.  The  results  of  section  4.5  indicate  that  the  extent  of 
chemical  reaction  across  the  weak  impinging  shock  wave  was  negligible.  If  we  also  assume  vibrational 
freezing  then  the  density  ratio  across  the  wave  is  given  by  the  perfect  gas  result, 

„  _  (7  ±  l)Af^sin^P 

(7-  l)M^sm2/?  +  2' 
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Figure  B.5:  Abel  transformation  of  density  field  for  flow  over  the  blunt  faced  cylinder  at  test  condi¬ 
tion  C;  Shot  T5-820. 


Test  condition 

A 

B 

c 

Peo  (SURF  kg/m") 

0.0218 

0.0155 

0.0157 

II 

9.9 

6.33 

5.31 

Shock  angle;  (3  (±1°) 

14.25 

15.0 

16.0 

P 

3.3 

2.1 

1.8 

Calculated  fringe  shift 

5.6 

1.9 

1.5 

Measured  fringe  shift  ±0.5 

4.5 

2. 

2. 

Table  B.5:  Density  rise  across  the  impinging  shock  wave. 


This  level  of  approximation  is  sufficient  given  the  large  uncertainty  in  the  interferometric  data 
when  the  fringe  shift  is  small.  Measurements  of  the  fringe  shift  that  are  indicated  in  table  B.5 
may  be  repeatably  obtained  from  any  of  the  holograms  contained  in  appendix  C.  Cases  A  and  B 
are  consistent  with  the  earlier  measurements  of  the  stagnation  point  density.  At  the  low  enthalpy 
condition  the  density  is  slightly  over-predicted  and  the  medium  enthalpy  condition  is  predicted  very 
accurately.  In  the  high  enthalpy  case  the  density  is  under-predicted  whereas  the  stagnation  point 
density  was  over-predicted.  This  suggests  that  the  nature  of  the  discrepancy  is  more  complicated 
than  a  simple  over-prediction  of  the  free-stream  density.  It  is  not  possible  to  draw  any  further 
conclusions  on  the  basis  of  the  current  measurements. 
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Test  condition 

ABC 

Moo  (t-1-4) 

Pfr  (IDG  model) 
peq  (IDG  model) 

9.9  6.33  5.31 

6.61  6.08  5.61 

6.61  8.59  9.67 

D  ref 

-X_ 

D  CQ 

Measured  ^  (±0.0025) 

0.202  0.217  0.228 

0.202  0.154  0.132 

0.218  0.186  0.173 

Table  B.6:  Experimental  verification  of  the  shock  standoff  distance. 

B.5  Experimental  verification  of  the  shock  standoff  distance 

An  additional  check  on  the  density  ratio  across  the  bow  shock  wave  is  provided  by  measurements  of 
the  shock  standoff  distance.  For  hypersonic  flows,  where  the  shock  layer  ahead  of  the  body  is  thin 


and  approximately  constant  density; 


D  2p 


Here  ^  is  the  shock  standoff  distance  normalized  with  respect  to  the  cylinder  diameter  and  p  is 
the  density  ratio  across  the  normal  portion  of  the  bow  shock.  This  correlation,  which  arises  from 
conservation  of  the  mass  flux  entering  the  shock  layer  across  the  bow  shock,  has  been  successfully 
extended  to  the  case  of  multicomponent  reacting  mixtures  by  Wen  [46].  Wen  used  the  classical 
perfect  gas  correlation  of  Ambrosia  &  Wortman  [1]  to  determine  the  similarity  constant  in  equation 
B.2.  Using  an  approximate  solution  for  the  nonequilibrium  flow  along  the  stagnation  streamline, 
Wen  [46]  determined  an  average  value  for  the  density  in  the  shock  layer  as  a  function  of  dimensionless 
kinetic  energy  and  reaction  rate  parameter. 

Here  we  consider  only  the  frozen  and  equilibrium  limits  of  such  a  correlation  for  the  flow  over 
cylinders.  Ambrosia  &  Wortman  [1]  accurately  correlated  the  shock  standoff  distance  from  a  cylinder 


by  the  expression, 


0.193  X 

Dre  f 


This  applies  in  the  frozen  limit  for  the  reacting  flow.  Accounting  for  the  increase  in  density  caused 
by  relaxation  using  the  IDG  solution  we  combine  equations  B.2  and  B.3  to  obtain, 


X  _  X  Pfr 

Deq  pQq 


(B.4) 


The  predicted  and  measured  values  for  the  shock  standoff  distance  are  indicated  in  table  B.6. 

At  low  enthalpies,  where  the  original  correlation  should  be  quite  accurate,  the  discrepancy  in  the 
shock  stand-off  distance  is  attributed  to  nonuniformity  of  the  nozzle  flow.  Reference  to  the  numerical 
results  of  Wen  [46]  suggests  that  the  reaction  rate  parameter  for  test  condition  B  (A  =  0.034  mm“^ 
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Test  condition 

ABC 

Stanton  number  (Equation  4.52) 

0.0123  0.0172  0.0167 

Table  B.7;  Predicted  stagnation  point  Stanton  numbers. 

X  40.6  mm  =  1.4)  is  sufficiently  large  at  moderate  enthalpy  to  cause  the  measured  standoff  distance 
to  approach  the  equilibrium  limit.  This  is  confirmed  here  by  the  measured  shock  standoff  distance 
and  is  consistent  with  the  holographic  data.  The  reaction  rate  parameter  for  the  high  enthalpy 
case  C  (A=21.),  combined  with  the  prediction  of  a  further  increase  in  the  shock  layer  density  ratio, 
should  cause  a  further  reduction  in  the  stand-off  distance.  Whilst  the  the  measured  shock  position 
moves  proportionally  closer  to  the  equilibrium  limit,  the  standoff  distance  remains  somewhat  larger 
than  that  predicted  by  the  correlation.  These  facts  are  consistent  with  the  over-prediction  of  the 
stagnation  point  fringe  shift  that  has  been  discussed  earlier. 

B.6  Stagnation  point  heat  transfer  rate 

The  stagnation  point  heat  transfer  rate  was  modeled  using  the  Fay  k  Riddell  [11]  equilibrium 
similarity  solution  in  the  form  discussed  in  section  4.7.  Predictions  of  the  Stanton  number  are  shown 
in  table  B.7.  Figures  C.l,  C.14  and  C.20  demonstrate  tolerable  agreement  with  this  correlation  at 


all  three  test  conditions. 
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Appendix  C  Holographic  interferograms  and  heat  transfer 
data 

The  data  that  are  presented  in  this  appendix  comprise  the  main  results  obtained  from  the  experi¬ 
mental  investigation  conducted  in  the  T5  free  piston  shock  tunnel  Experiments  were  conducted  at 
the  three  test  conditions  that  are  described  in  appendix  B,  For  each  shot  a  holographic  interferogram 
is  presented  along  with  the  heat  transfer  distribution  around  the  circumference  of  the  model.  The 
data  have  been  arranged  by  test  condition  and  in  order  according  to  the  position  of  the  model  Due 
to  slight  variations  in  the  incident  shock  angle,  some  of  the  interferograms  appear  to  be  out  of  order. 
The  position  of  the  model  was  measured  by  gauging  the  distance  between  flat  reference  surfaces  on 
the  model  traverse.  Values  of  this  distance,  g/D,  that  are  indicated  for  each  of  the  figures,  have  been 
normalized  by  the  model  diameter.  The  zero  point  of  this  quantity  was  arbitrary  but  repeatable  to 
within  ±0.001. 

The  effects  of  variation  of  the  shock  impingement  location  were  visualized  using  the  double 
pulse  holographic  interferometer  that  is  described  in  appendix  D.  Quantitative  interpretation  of 
the  density  field  is  possible  and  the  necessary  physical  constants  are  contained  in  table  C.l.  The 
reader  will  observe  some  limitations  in  the  quality  of  images  that  were  obtained  during  the  early 
development  of  the  interferometer.  They  are  presented  here  for  completeness.  Repeated  magnified 
views  are  also  presented  where  appropriate. 

Time  resolved  heat  transfer  measurements  were  made  at  24  locations  around  the  circumference 
of  the  model.  Geometrical  details  of  the  test  section  apparatus  and  sensor  locations  are  contained 
in  appendix  F.  Sensor  locations  have  been  overlayed  on  the  holographic  interferograms  and  the 
dimensions  of  the  symbols  correspond  to  the  sensitive  area  of  the  thermocouples.  The  instantaneous 
heat  transfer  distribution  was  plotted  for  each  interferogram  along  with  an  average  value  taken  over 
an  interval  of  width  100  fis.  The  instantaneous  profile  and  the  center  of  the  averaging  interval 
correspond  to  the  time  at  which  the  hologram  was  recorded.  Superimposed  on  these  plots  are  the 
standard  deviations  of  the  time  histories  over  the  averaging  interval  This  gives  an  indication  of  the 


Wavelength 
Path  length 

Gladstone-Dale  constants 


A  =  532nm 
L  =  0.180m 

AA.v  =  3.1  X  10“^  m-n^g-^ 

=  2.45  X  IQ-^^  m’^kg-^ 


Table  C.l:  Physical  constants  for  holographic  interferograms. 
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temporal  fluctuations  of  the  heat  transfer  rates  at  various  locations  on  the  model  surface. 

The  final  set  of  points  on  each  plot  is  the  envelope  of  the  maximum  heat  transfer  rate  recorded 
at  each  circumferential  location  during  the  shot.  The  first  500  /xs  of  the  time  records  were  excluded 
from  this  maximization  so  that  tunnel  starting  transients  did  not  influence  the  profiles.  The  records 
were  scanned  until  4  ms  after  primary  shock  reflection.  This  interval  spans  the  period  of  useful  test 
flow.  Whilst  the  utility  of  a  maximum  level  measurement  inevitably  depends  on  the  temporal  and 
spatial  resolution  of  the  device,  these  factors  remained  constant  for  all  data  presented  here.  This 
statistic  has  proven  useful  for  the  interpretation  of  the  effect  of  variation  of  the  stagnation  enthalpy. 
Further  details  of  the  specialized  techniques  used  to  obtain  this  data,  along  with  an  assessment  of 
the  accuracy  of  the  sensors,  are  contained  in  appendix  E. 
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Figure  C.2:  Shot  T5-844;  condition  A,  g/D  =  2.128. 
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Figure  C.3:  Shot  T5-845;  condition  g/D  ~  2.254. 
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Shot  T5-849 
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Figure  C.4:  Shot  T5-849;  condition  A,  g/D  =  2.339. 
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Figure  C.5:  Magnified  view  of  holographic  interferogram;  Shot  T5-849 
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Figure  C.7;  Magnified  view  of  holographic  shadowgraph;  Shot  T5-854 
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Figure  C.IO:  Shot  T5-850;  condition  k,  g/D  —  2.396. 
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Figure  C.ll:  Magnified  view  of  holographic  interferogram;  Shot  T5-850 
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Figure  C.14:  Shot  T5-856;  condition  B,  g/D  —  2.185. 
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Figure  C.17:  Shot  T5-852;  condition  B,  gjD  =  2.197. 
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Figure  C.18;  Magnified  view  of  holographic  interferogram;  Shot  T5-852 


Shot  T5-851 


Figure  C.19:  Shot  T5-851;  condition  B,  g/D  =  2.281. 
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Figure  C.21:  Shot  T5-840;  condition  C,  gjD  =  1.361. 


C-24 


Shot  T5-834 
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Figure  C.22:  Shot  T5-834;  condition  C,  g/D  =  1.485. 
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Figure  C.23:  Shot  T5-833;  condition  C,  g/D  =  1.622. 
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Figure  C.24:  Shot  T5-825;  condition  C,  g/D  =  1.926. 
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Figure  C.25:  Shot  T5-826;  condition  C,  gjO  =  2.036. 


Figure  C.26;  Shot  T5-838;  condition  C,  g/D  =  2.113. 
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Figure  C.27:  Magnified  view  of  holographic  interferogram;  Shot  T5-838 
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Figure  C.29:  Magnified  view  of  holographic  interferogram;  Shot  T5-831 
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Figure  C.30:  Shot  T5-836;  condition  C,  g/D  =  2.162. 
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Figure  C.31:  Shot  T5-827;  condition  C,  g/D  =  2.179. 
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Figure  C.32;  Shot  T5-837;  condition  C,  g/D  =  2.189. 
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Figure  C.34:  Shot  T5-830; 
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C,  g/D  =  2.198. 
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Figure  C.36:  Shot  T5-839;  condition  C,  gjD  =  2.211. 
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Figure  C.38:  Shot  T5-829;  condition  C,  g/D  =  2.246. 
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Figure  C.39:  Magnified  view  of  holographic  interferogram;  Shot  T5-829 
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Figure  C.40:  Shot  T5-828;  condition  C,  g/D  =  2.274. 
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Figure  C.42:  Shot  T5-832;  condition  C,  g/D  =  2.307. 
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Figure  C.44:  Shot  T5-841;  condition  C,  g/D  =  2.620. 


Figure  C.47:  Flow  over  a  ^50mm  blunt  cylinder  without  shock  impingement.  Shot  T5-819;  condi¬ 
tion  A. 
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Figure  C.49:  Flow  over  a  (/>50mm  blunt  cylinder  without  shock  impingement.  Shot  T5-820;  condi¬ 
tion  C. 
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Figure  C.50;  Flow  over  a  (/)50mm  blunt  cylinder  with  type  IV  shock  impingement.  Condition  C. 
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Appendix  D  Holographic  interferometer 

The  main  features  of  the  design  of  this  device  were  a  consequence  of  the  desire  to  provide  maximum 
resolution  of  features  at  the  test  section.  Consequently,  short  focal  length  multi-element  lenses  were 
used  for  the  imaging  optics  and  these  were  placed  as  close  as  possible  to  the  test  section.  The 
device  also  incorporates  an  unusually  small  crossing  angle  between  the  test  section  and  reference 
beams.  This  allowed  the  use  of  relatively  low  resolution  photographic  media  and  short  coherence 
length  lasers.  Diffraction  minima  about  the  images  of  shock  waves  were  suppressed  by  the  use  of 
Fourier  plane  conjugate  image  separation  and  Gaussian  apodizing  filters.  White  light  reconstruction 
suppressed  the  appearance  of  laser  speckle. 

D.l  Description  of  optical  apparatus 

The  optical  arrangement  is  shown  in  Figure  D.l  and  the  technique  used  was  off  axis  reference  beam 
image  plane  holography.  A  frequency  doubled  NdYAG  laser  (A=532  nm)  with  a  coherence  length  of 
10  mm  and  30  ns  pulse  duration  was  used  as  the  light  source.  The  source  beam  was  passed  vertically 
downward  into  the  top  of  the  interferometer  by  a  periscope  arrangement.  A  dielectric  cube  beam 
splitter  then  separated  the  reference  and  test  section  beams. 

The  test  section  beam  was  expanded  through  a  diverging-converging  singlet  lens  pair  and  folded 
to  pass  through  the  test  section.  A  doublet  lens  focused  the  modulated  beam  so  that  it  passed 
through  an  aperture  stop  of  diameter  d?ii5  mm.  Since  this  stop  was  an  image  at  infinity  with 
respect  to  the  test  section,  the  /-number  was  large  (f/d?5il20)  with  respect  to  the  diffuse  test  section 
luminosity.  An  interference  filter  (AA=10  nm)  further  reduced  the  background  level.  A  second 
doublet  objective  re-collimated  and  minified  the  beam  before  it  was  recombined  with  the  reference 
beam  at  the  plate  beam  splitter. 

The  unexpanded  reference  beam  was  reflected  through  a  path  length  adjuster  before  being  passed 
around  the  test  section  through  the  interior  of  the  interferometer  support  beam.  Several  beam  steer¬ 
ing  prisms  were  used  to  align  the  beam  with  the  singlet  pair  reference  beam  expander.  The  reference 
beam  formed  a  5®  angle  as  it  interfered  with  the  test  section  beam.  Image  plane  holograms  were  ob¬ 
tained  by  focusing  the  test  section  onto  the  film  plane.  The  rejected  beams  emerging  from  the  plate 
beam  splitter  were  arranged  to  fall  on  a  ground  glass  screen  placed  in  front  of  an  exposure  control 
photo  diode.  A  rotating  dichroic  polarizer  was  mounted  on  the  camera  body  to  regulate  the  ratio 
of  the  beam  intensities.  Custom  built  programmable  digital  timing  circuitry  was  used  to  thermally 
stabilize  the  laser  prior  to  the  shot  and  to  synchronize  the  laser  pulses  and  electromechanical  shutter 
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when  the  tunnel  was  fired.  The  two  interferometric  pulses  were  generated  from  separate  flash-lamp 
cycles. 

Holograms  w’ere  recorded  on  Kodak  SO-253  plastic  based  film  that  was  mounted  in  a  modified  35 
mm  camera  body.  Film  development  was  optimized  to  best  accommodate  the  double  pulses  within 
the  linear  portion  of  the  emulsion  characteristic  curve.  The  developed  images  were  bleached  in  a 
solution  of  potassium  permanganate  and  sulfuric  acid  to  produce  phase  modulated  holograms  with 
improved  diffractive  efficiency. 

By  replacing  the  developed  hologram  at  its  position  during  exposure,  the  reference  beam  was 
used  to  reconstruct  the  hologram.  A  folding  mirror  passed  the  reconstructed  beams  through  a  pair 
of  lenses  that  imaged  the  hologram  at  unity  magnification  onto  a  secondary  camera  mount  where  the 
interferogram  was  recorded.  Multiple  diffracted  orders  produced  by  the  hologram  were  separated 
by  an  aperture  placed  at  the  focus  of  the  two  lenses.  Low  power  densities  during  reconstruction 
allowed  the  use  of  an  aperture  with  a  Gaussian  transmittance  profile.  These  filters  were  produced 
photographically.  It  may  be  shown  that  this  profile  produces  non-oscillatory  images  of  the  shock 
fronts  that  are  free  from  diffraction  effects.  The  interferograms  were  reconstructed  using  white  light 
to  reduce  the  appearance  of  laser  speckle.  Conventional  photographic  methods  were  used  to  record 
and  print  the  reconstructed  interferograms. 

The  interferometer  was  assembled  in  a  U-shaped  plywood  cabinet  to  provide  appropriate  damp¬ 
ing  characteristics.  The  long  term  static  stability  of  the  device  is  unimportant  since  double  pulse 
holography  is  sensitive  only  to  path  length  differences  during  recording.  An  aluminum  truss  across 
the  top  of  the  box  allowed  access  to  the  internal  components  whilst  maintaining  the  stiffness  of  the 
housing.  Adjustable  rails  were  used  to  position  the  interferometer  over  the  test  section  of  the  shock 
tunnel.  Optical  components  were  rigidly  attached  to  the  baseboard  by  custom  made,  nonadjustable 
mounts  that  were  permanently  aligned  during  the  initial  set  up  of  the  device.  The  interferometer  is 
easily  transportable  and  self  contained  except  for  the  external  laser. 

D.2  Modes  of  operation 

This  system  was  used  in  a  single  shot  mode  for  schlieren  and  shadowgraph  reconstruction  and  princi¬ 
pally  in  a  double  pulse  holographic  interferometry  mode.  In  this  second  technique  an  initial  hologram 
is  recorded  with  no  flow  through  the  test  section.  An  additional  hologram  is  recorded  on  the  same 
film  at  an  appropriate  time  during  the  flow.  The  two  waveforms  are  reconstructed  simultaneously 
to  produce  an  interferogram  that  displays  the  path  length  difference  due  to  the  refractivity  of  the 
gas  in  the  test  section.  Note  that  this  procedure  does  not  provide  a  three-dimensional  capability. 
Here  holography  is  used  to  facilitate  the  sequential  recording  and  simultaneous  reconstruction  of 
two  wavefronts.  The  results  that  are  produced  are  functionally  equivalent  to  a  Mach-Zehnder  inter- 
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Figure  D.l:  Optical  and  mechanical  layout  of  the  holographic  interferometer. 

ferogram  whereby  the  information  is  integrated  along  the  line  of  sight.  The  only  three-dimensional 
capability  is  the  freedom  to  choose  the  plane  of  focus  during  reconstruction.  This  remaining  freedom 
is  removed  if  incoherent  light  is  used  during  reconstruction. 

Introducing  a  diffuser  plate  immediately  before  to  the  test  section  produces  an  infinite  number 
of  point  sources.  This  provides  three-dimensional  information  since  the  point  of  view,  the  plane 
of  focus  and  the  angular  field  of  view  (and  hence  depth  of  focus)  may  be  varied  arbitrarily  during 
reconstruction.  The  information  obtained  by  this  technique  is  averaged  along  the  line  of  sight  and 
also  over  a  conical  volume  either  side  of  the  plane  of  focus.  The  conical  volume  is  determined  by 
the  solid  angle  subtended  by  the  aperture  of  the  reconstruction  imaging  optics.  In  conventional 
holography  of  solid  objects  this  simply  reduces  the  depth  of  field;  however  for  phase  objects  a  finite 
path  length  is  still  required  in  order  to  generate  an  image.  Thus  three-dimensional  information  is 
obtained  at  the  expense  of  additional  angular  mollification  of  the  fine  flow  structure.  This  point 
is  not  well  treated  in  the  literature  and  is  of  considerable  importance  in  obtaining  high  resolution 
images. 


D.3  Principles  of  operation 

The  physical  principles  of  holography  are  described  in  the  literature  (Caulfield  [4],  Merzkirch  [32]). 
Here  we  demonstrate  the  techniques  of  white  light  reconstruction  and  Fourier  plane  conjugate  image 
separation  that  are  features  of  the  current  device. 

Consider  the  interference  of  the  plane  test  section  wavefront,  T,  with  the  reference  beam,  Ro 


D-4 


(Figure  D.2).  The  spatial  variations  of  the  two  fields  in  the  z  =  0  hologram  plane  at  a  given  point 
in  time  are  described  by, 

T  —  Rq  — 

where  <j)  and  (j}o  are  the  phases  at  a  given  instant  in  time.  Both  beams  are  monochromatic  with 
wavenumber,  ko  -  27r/Ao.  The  reference  beam  is  planar  and  forms  an  angle  0  with  respect  to  the 
test  section  beam  in  the  x-z  plane  so  that. 


<f)o(x,  y)  =  koxsin  0. 

The  modulation  of  the  phase  of  test  section  beam  by  the  gas  inhomogeneities  is  described  hy  <i){x,y). 
At  the  hologram  plane  the  media  responds  to  the  intensity  of  the  interfering  wavefronts  yielding  an 
image  of  density, 

lH{x,y)  -  \T^Rof  =  +  ^ikoo:  sin  =  2  [1  +  cos{<P  -  koxsinO)] .  (D.l) 

Holographic  fringes  are  therefore  formed  on  the  film  according  to, 

27rn  —  <j){x,  y)  -  kox  sin  0. 

In  the  case  of  a  uniform  test  section  beam,  =  0,  we  see  that  the  holographic  fringes  are  formed 
with  spacing,  5  2iT/koO,  for  small  crossing  angle,  0.  Phase  variations  are  recorded  as  a  modulation 
of  the  local  fringe  spacing. 

Consider  reconstruction  of  the  hologram  with  a  reference  wave  of  the  same  geometry  but  ar¬ 
bitrary  wavelength,  R  =  ^^k^vsme  ^  image  formed  on  the  film  modulates  the  amplitude  of  the 

reconstruction  beam  to  produce  the  following  wavefront  immediately  downstream  of  the  hologram, 


D{x^  y)  =  2  [1  +  cos((f)  -  koX  sin  ^)]  e 


ikxsind  _  2^'^kxsin$  sin  0)  _|_  gZ(  — +  sin 


The  first  term  corresponds  to  the  transmitted  component  of  the  reconstruction  beam.  The  second 
term  is  the  reconstructed  test  section  beam  and  it  is  deflected  by  an  angle  that  is  proportional  to 
the  shift  in  wavenumber  of  the  reconstruction  beam.  An  additional  wave,  known  as  the  conjugate 
wave,  is  symmetrically  disposed  with  respect  to  the  reference  beam  and  is  also  deflected  because  of 
the  difference  in  reconstruction  wavelength.  Additional  higher  order  diffracted  beams  will  be  present 
since  the  actual  recording  process  will  deviate  from  this  ideal.  Each  of  these  orders  is  brought  to 
focus  by  the  first  lens  shown  in  Figure  D.2.  By  placing  an  aperture  at  the  focus  corresponding  to 
the  reconstructed  test  section  beam,  the  additional  unwanted  images  are  removed.  Each  of  the  foci 
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broaden  as  the  spatial  frequencies  recorded  on  the  hologram  increase.  Overlap  of  the  adjacent  orders 
is  prevented  by  the  spatial  filter  used  during  recording  of  the  hologram.  Maximum  resolution,  for  a 
given  crossing  angle,  is  achieved  with  apertures  whose  diameters  are  equal  to  the  spacing  between  the 
foci.  Note  that  background  luminosity  determines  the  aperture  diameter  and  so  the  small  crossing 
angle  does  not  limit  the  resolution  of  the  interferometer.  Thus  spatial  filtering  recovers  the  second 
term  of  equation  D.2  which  is  the  reconstructed  test  section  beam, 

^  — sin  0) 

Double  pulse  holography  interferes  two  such  beams.  When  the  second  lens  images  the  hologram 
onto  the  interferogram  recording  plane  the  original  phase  distribution  is  recovered  and  we  have, 

^i(4>o+(k-ko)x  sm$) 

Since  we  have  recorded  an  image  plane  hologram  the  phase  differences  caused  by  propagation  between 
the  hologram  and  the  interferogram  planes  vanish.  At  the  interferogram  recording  plane  the  response 
of  the  media  yields  an  image  of  density, 

^  2 

Irix^y)  Ti  +  T2  2  [1  +  cos(((^2  +  -  A:o)xsin^)  -  {(pi  +  (^  -  kQ)xsmO))]  . 

An  additional  consequence  of  image  plane  recording  is  that  the  effect  of  the  wavelength  shift  is 
common  to  both  beams  and  cancels.  Thus  the  hologram  may  be  reconstructed  with  white  light  and 
the  fringe  shift  in  the  final  interferogram  is  independent  of  the  reconstruction  wavelength, 

If{x,  y)  =  2[l-\-  cos((/)2  -  0i)] . 

Fringes  are  observed  in  the  interferogram  at  locations, 

27rm  =  <^2  ■”  • 

Coherence  of  the  two  polychromatic  beams  is  ensured  by  accurate  focusing  and  the  use  of  achromatic 
lenses  for  reconstruction. 
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Figure  D.2;  Geometry  of  hologram  reconstruction  and  conjugate  image  separation. 
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Appendix  E  Surface  junction  thermocouple  sensors 
E.l  Introduction 

The  measurement  of  aerodynamic  heat  transfer  rates  is  important  in  the  design  of  thermal  protection 
systems  for  hypervelocity  aerospace  vehicles.  High  enthalpy  flows  may  be  created  in  wind  tunnels 
for  only  a  fraction  of  a  second  because  of  the  high  specific  energies  involved  and  the  limited  strength 
of  materials  at  elevated  temperatures.  Thus  the  problem  is  transient  in  nature  and  a  common 
measurement  technique  is  to  sense  the  time  history  of  the  surface  temperature  of  the  model.  The 
unsteady  conduction  of  heat  within  a  solid  body  is  described  by  linear  partial  differential  equations 
and  so  the  heat  transfer  rates  at  the  surface  of  a  body  may  be  inferred  from  the  surface  temperature 
measurements.  The  operating  principles,  and  typical  sensor  designs  are  described  by  Schultz  & 
Jones  [40]. 

We  have  conducted  experiments  in  the  T5  free  piston  driver  shock  tunnel  at  GALCIT.  This 
facility  is  capable  of  simulating  the  high  density  hypervelocity  flows  about  space  vehicles  at  speeds 
of  order  5000  m/s.  Under  these  extreme  conditions  the  flow’  in  typical  facilities  is  contaminated  with 
the  products  of  mechanical,  physical  and  chemical  breakdown  of  various  materials.  Some  erosion 
of  the  surface  of  the  gauge  is  probable  under  these  conditions  and  this  impacts  the  design  of  the 
transducer. 

A  commonly  used  technique  for  gauge  fabrication  is  to  deposit  a  thin  metal  film  (typically 
Platinum)  on  a  glass  or  ceramic  substrate.  The  surface  temperature  is  then  sensed  through  the 
temperature  dependence  of  the  resistivity  of  the  metal  film.  This  technique  yields  a  relatively 
sensitive  device  however  the  thin  film  is  prone  to  damage  in  hostile  environments.  Any  erosion  of 
the  metal  film  alters  the  sensitivity  of  the  gauge  and  so  these  devices  are  unsuitable  for  some  wind 
tunnels.  Observe  also  that  since  the  sensitivity  of  the  gauge  depends  on  the  geometry  of  the  device, 
calibration  of  each  individual  sensor  is  required. 

An  alternative  physical  mechanism  for  temperature  measurement  is  the  thermoelectric  EMF 
produced  at  a  junction  of  dissimilar  metals.  If  a  metallic  substrate  is  used,  the  gauge  is  capable  of 
tolerating  a  larger  transient  heat  flux  due  to  the  higher  thermal  diffusivity  and  the  larger  allowable 
temperature.  The  figure  of  merit  is  TmaxVp^-  A  particularly  robust  design  results  when  the 
thermocouple  material  itself  is  used  as  the  substrate.  The  sensitivity  of  the  gauge  is  then  unaffected 
by  surface  erosion  or  gauge  geometry  and  depends  only  on  physical  properties  of  the  materials.  Thus 
calibration  of  individual  gauges  is  not  required. 
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The  design  problem  in  producing  such  a  gauge  is  to  insulate  the  two  conductors  that  form  the 
thermocouple  so  that  only  a  thin  surface  junction  is  formed.  Common  practice  is  to  construct 
the  gauge  from  coaxial  cylindrical  thermocouple  elements  that  are  insulated  by  a  layer  of  epoxy  or 
ceramic  adhesive.  The  surface  junction  is  formed  by  electrolytic  deposition  of  a  metal  film  or  by 
abrasion  of  the  exposed  surface.  In  either  case  very  fine  tolerances  are  required  to  produce  a  small 
gap  between  the  thermocouple  elements  and  the  effective  depth  of  the  junction  scales  with  the  width 
of  the  gap  between  the  conductors. 


E.2  Analysis  of  heat  conduction  in  gauge  substrate 

In  order  to  illustrate  the  principles  of  gauge  operation  and  to  assess  the  effect  of  the  finite  thickness 
of  the  surface  junction,  consider  the  unsteady  conduction  of  heat  in  a  one-dimensional  semi-infinite 
solid.  For  small  temperature  gradients  this  process  is  described  by  the  usual  linear  partial  differential 


equation, 


d-T  __ 

dx^  a  dt  ’ 


where  T  is  the  temperature,  x  is  the  spatial  coordinate  normal  to  the  surface,  t  is  time  and  a  is  the 
thermal  diffusivity.  Note  that  a  =  where  k  is  the  thermal  conductivity,  p  is  the  density  and  c 
is  the  specific  heat  of  the  material.  For  the  following  uniform  initial  condition,  and  instantaneously 
applied,  constant  surface  heat  flux,  qo] 

T{x,0)  =  Ti, 


the  solution  is, 


AT  ^T-Ti^ 


2qo  ^ 

xV  V 


—X  qox 


1  ~  erf 


2V^J  ' 


Expanding  this  solution  asymptotically  in  powers  of  /  as  f  ^  oo  at  constant  x  we  obtain. 


AT= 

k  y  TT  k 


2  k  V  / 


The  first  term  is  asymptotically  dominant  for  — <<  1- 

(at)  2 


and  this  expresses  the  usual  result  that  the  surface  temperature  rises  with  \/t  for  a  constant  heat 
flux.  The  criterion  — ^  «  1  provides  a  limit  on  the  thickness  of  the  surface  sensing  element  for 

(at)2 

practical  designs.  The  next  order  term,  ^x,  gives  the  error  in  the  surface  temperature  measurement 


E-3 


4 

3 


2 


1 


0 

0  2  4  6  8  10 

Time 

Figure  E.l:  Asymptotic  behavior  of  a  surface  temperature  sensor.  For  illustrative  pur¬ 
poses  a=l,k  =  ljqo=^liX  =  l. 


associated  with  this  finite  thickness  and  the  estimate  is  also  asymptotically  valid  provided  — ^  << 

(at)2 

1.  Extracting  the  exponential  prefactor  and  expanding  asymptotically  as  f  0  at  constant  x,  the 
leading  order  behavior  is; 


The  significance  of  the  various  terms  of  the  expansion  is  illustrated  in  figure  E.l.  At  long  times  the 
solution  asymptotes  to  the  true  surface  temperature,  given  by  the  first  term  in  the  expansion  E.3, 
with  a  constant  error  given  by  the  second  term. 


E.3  Description  of  new  gauge  design 

A  new  gauge  design  has  been  developed  and  tested  with  the  aim  of  simplifying  the  construction  of 
surface  junction  thermocouple  sensors.  The  essential  feature  of  this  design  is  the  use  of  a  slightly 
tapered  center  conductor.  The  resulting  interference  between  the  tapered  center  pin  and  the  sharp 
edged  outer  conductor  forms  a  thin  surface  junction  and  provides  reliable  insulation  (Figure  E.2). 
The  protruding  pin  is  severed  and  the  surface  is  polished  once  the  gauge  has  been  assembled.  This 
design  results  in  a  robust  transducer  whose  sensitivity  is  determined  solely  by  the  physical  properties 
of  the  thermocouple  materials.  There  is  a  further  weak  dependence  on  the  taper  angle  however  this 
is  constant  for  all  gauges  and  is  determinable  by  either  calibration  or  analysis. 

The  effective  thickness  of  the  junction  may  be  estimated  from  the  surface  roughness  of  the 
thermocouple  components  and  the  taper  angle  of  the  center  conductor  (Figure  E.3).  The  thickness 
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Figure  E.2:  Layout  of  new  surface  junction  thermocouple  sensor. 


Figure  E.3:  Effective  thickness  of  thermocouple  junction. 


of  the  junction,  averaged  over  the  two  surfaces  may  be  estimated  to  be, 


thickness  =  taper  x  roughness. 


It  is  possible  to  achieve  a  0.3  fJ.  finish  with  conventional  machining  operations  and  with  a  5:1  taper 
the  average  junction  depth  is  1.5  /im.  For  the  gauges  described  here  type  E  (Chromel-Constantan) 
thermocouples  were  used  and  the  average  thermal  diffusivity  for  these  materials  is  5  x  10“®  m^/s.  A 
response  time  of  0.5  us  follows  from  the  condition  — ^  ~  1.  The  test  time  for  typical  hypervelocity 

{at)2 

wind  tunnels  lies  in  the  range  1-10  ms  and  so  this  was  acceptable.  The  length  of  the  gauge  must 
be  sufficient  so  that  the  temperature  rise  at  the  rear  surface  is  negligible.  This  is  determined  from 
the  condition  — >>  1;  where  here  the  scales  are  understood  to  be  the  test  time  and  the  length 

(ai)2 

of  the  gauge. 

A  second  dissimilar  metal  junction  is  formed  between  the  outer  conductor  and  the  model.  This 
junction  extends  along  the  full  length  of  the  gauge  since  these  surfaces  are  not  insulated  and  so  the 
thermoelectric  EMF  is  small.  In  the  current  experiments  this  common  mode  EMF  and  noise  pickup 
were  reduced  by  the  use  of  a  simple  differential  amplifier  (Figure  E.4).  The  cold  junction  was  formed 
at  ambient  temperature  at  the  amplifier  circuit  board.  This  arrangement  was  adequate  given  the 
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Figure  E.4:  Schematic  of  amplifier  circuit. 


constancy  of  the  cold  junction  temperature  during  the  short  run  time  and  the  large  hot  junction 
temperature  rise  experienced  in  hypervelocity  wind  tunnels  (^200K). 

The  center  pin  was  located  in  the  concentric  outer  conductor  and  secured  with  epoxy  or  ceramic 
adhesive.  Sufficient  force  was  applied  to  ensure  adequate  registration  of  the  mating  parts  however 
care  was  taken  not  to  deform  the  sharp  edge  of  the  outer  conductor.  Thermocouple  extension  wires 
were  resistance  welded  to  the  rear  of  the  assembly  by  discharging  a  capacitor  through  the  wire. 
The  thermocouple  was  then  installed  flush  with  the  outer  surface  of  the  model  and  the  protruding 
center  pin  was  removed.  The  excess  material  was  removed  by  filing  and  the  surface  was  honed  to 
produce  the  surface  junction.  Some  degradation  of  the  adhesive  bond  between  the  conductors  was 
anticipated  however  the  thermal  penetration  depth  during  the  test  time  was  small  so  that  the  bond 
at  the  rear  maintained  the  integrity  of  the  device. 


E.4  Deducing  the  surface  heat  transfer  rate 


If  the  unsteady  conduction  of  heat  in  the  thermocouple  sensor  is  modeled  by  the  analysis  presented 
in  section  E.2  then  given  the  time  history  of  the  surface  temperature  and  the  physical  properties 
of  the  gauge  materials  it  is  possible  to  solve  an  inverse  problem  for  the  time  history  of  the  surface 
heat  flux.  The  procedure  is  described  by  Schultz  &  Jones  [40]  whereby  the  calculation  is  performed 
using  a  problem  specific  quadrature  scheme  or  by  an  electronic  analogue. 

Observe  however  that  the  governing  equations  are  linear  and  hence  the  dynamic  behavior  of  the 
system  may  be  described  by  a  convolution  integral, 

AT{x,i)=  f  g(x,t  -  r)q(T)dT.  (E.6) 

Jo 


The  unit  impulse  response  function,  g{x,t),  is  given  by  the  derivative  of  equation  E.2, 

,  ,  dAT(x^t)  I  a  —x^  ^ 
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Note  that  the  integral  E.6  is  improper  in  the  limit  a;  — 0  and  a  finite  value  of  x  is  taken  to  regularize 
the  problem  (as  estimated  in  section  E.2).  Recasting  the  problem  in  this  general  form  allows  the 
use  of  spectral  deconvolution  methods  that  provide  a  mean  square  optimal  estimate  of  the  time 
history  of  the  surface  heat  flux.  This  procedure  provides  optimal  rejection  of  the  errors  associated 
with  the  differentiation  of  noisy  experimental  data.  Theoretical  background  is  given  in  Papoulis  [34] 
and  the  numerical  implementation  with  discrete  Fourier  transform  algorithms  is  described  by  Press 
et  al.  [37].  A  concise  description  of  these  methods  applied  to  an  analogous  problem  is  given  by 
Sanderson  &  Simmons  [39]. 

If  the  noisy  experimental  data  in  discrete  time  is. 


Zi  —  Tli  j 


(E.8) 


then  the  deconvolution  problem  comprises  determining  the  acausal  filter,  <j>i,  that  gives  a  mean 
square  optimal  estimate,  q*,  of  the  time  history  of  the  surface  heat  flux.  Thus, 

J^<l^i-kZk  =  j29i-kql;  f  =  0,l...A-l.  (E.9) 

k=0  k=0 

By  the  convolution  property  of  the  discrete  Fourier  transform,  and  since  jf;  =  0  for  i  <  0, 


^jZj=GjQ*  j  =  0,l...N-l. 


(E.IO) 


Upper  case  symbols  denote  the  discrete  Fourier  transform  of  the  corresponding  lowercase  quantities. 
It  may  be  shown  that  the  optimal  filter  is  determined  by. 


Szz^j  —  Vj 
SzZi 


(E.ll) 


where  SzZj  is  the  discrete  power  spectral  density  of  the  surface  temperature  measurement  and  Vj 
is  the  power  spectral  density  of  the  signal  noise.  By  inverse  discrete  Fourier  transform,  the  optimal 
estimate  of  the  time  resolved  surface  heat  flux  is  then. 


^jZj 


G 


3  J 


i,j  =  0,l...A-l. 


(E.12) 


Note  that  all  time  series  must  be  padded  to  twice  their  length  with  trailing  zeroes  to  enforce  causality 
of  the  right-hand  side  of  equation  E.9  when  discrete  Fourier  transforms  are  used  in  this  manner. 
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Figure  E.5:  Stagnation  point  heat  transfer  measurement  for  sphere  in  hypervelocity 
FLOW  USING  DUMMY  CONSTANTAN  -  CONSTANTAN  THERMOCOUPLE  AT  CONDITION  C.  The  large 
drop  in  the  pitot  pressure  at  t  —  1.5ms  was  caused  by  the  passage  of  an  oblique  shock  wave. 


E.5  Accuracy  of  measurements 

It  has  been  demonstrated  that  the  sensitivity  of  the  transducer  depends  only  on  the  thermo-physical 
properties  of  the  materials.  All  transducers  used  in  the  current  experiments  were  manufactured 
from  the  same  piece  of  raw  material.  Comparison  of  symmetrically  disposed  transducers  on  the 
cylindrical  model  indicated  that  the  repeatability  of  the  measurement  for  individual  realizations  of 
the  design  was  within  ±10%. 

The  flow  environment  in  T5  was  extremely  hostile  and  since  the  thermocouple  junction  was 
electrically  exposed  to  the  flow  it  was  necessary  to  demonstrate  that  erroneous  EMFs  were  not 
generated  by  physical  mechanisms  other  than  the  thermoelectric  effect.  A  dummy  gauge  of  identical 
design  was  manufactured  with  both  conductors  made  of  Constantan  so  that  no  thermoelectric  EMF 
was  generated.  Figure  E.5  shows  the  results  obtained  with  this  dummy  gauge  at  a  high  enthalpy 
condition.  Apart  from  an  initial  transient  that  decays  rapidly,  no  erroneous  EMFs  are  discernible 
on  this  scale.  The  deconvolved  heat  transfer  rate  is  three  orders  of  magnitude  less  than  the  actual 


heat  transfer  rate. 
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Appendix  F  Technical  data 


